BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,292,012	T→C	5.5%	*855R (TGA→CGA)	USA300HOU_RS11675 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,292,012	0	T	C	5.5%	418.3 / 8.8	146	*855R (TGA→CGA)	USA300HOU_RS11675	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (58/78); new base C (5/3); total (63/83)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 150 total aligned reads displayed.

AAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGTTG  >  USA300TCH1516_ALE/2291873‑2292152

aaaGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTg                                                                                                                                             <  1:1016766/141‑1 (MQ=34)

 aaGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCGATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGa                                                                                                                                            >  1:802396/1‑141 (MQ=21)

  aGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGaa                                                                                                                                           <  1:290984/141‑1 (MQ=34)

    tCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAAtt                                                                                                                                         <  1:42539/141‑1 (MQ=35)

     cTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAAttt                                                                                                                                        >  1:617325/1‑141 (MQ=34)

      tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc                                                                                                                                       >  2:138654/1‑141 (MQ=35)

      tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc                                                                                                                                       >  1:345215/1‑141 (MQ=35)

      tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc                                                                                                                                       <  2:137369/141‑1 (MQ=35)

      tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc                                                                                                                                       >  1:596515/1‑141 (MQ=35)

        tgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCtt                                                                                                                                     >  2:635745/1‑141 (MQ=35)

         gAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCttt                                                                                                                                    >  1:397960/1‑141 (MQ=35)

         gAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCttt                                                                                                                                    >  4:298086/1‑141 (MQ=35)

           gTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTg                                                                                                                                  <  2:4998/141‑1 (MQ=35)

                  tgtATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAg                                                                                                                           <  2:740487/141‑1 (MQ=35)

                   gtATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGa                                                                                                                          >  1:98541/1‑141 (MQ=35)

                   gtATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGa                                                                                                                          >  2:64233/1‑141 (MQ=35)

                     aTGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGaaa                                                                                                                        <  2:323238/141‑1 (MQ=35)

                            gTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTcaac                                                                                                                 <  1:459661/141‑1 (MQ=34)

                            gTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTcaac                                                                                                                 <  4:244101/141‑1 (MQ=34)

                             tGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTcaaca                                                                                                                >  3:300387/1‑141 (MQ=33)

                             tGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTcaaca                                                                                                                <  1:596859/141‑1 (MQ=33)

                              ggTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTcaacaa                                                                                                               <  1:1066873/141‑1 (MQ=32)

                               gTTAGATAATATCGAAGACATCAAAGGATGTTAAGAAATACAATATCTTACCCAGCATTTAATAATGGTGCAATAGAAGGAATTAATAATAAGATAAAATTAATAAAGTGAATTTCTTTCGGTTACAGAAATTTCaacaac                                                                                                              >  4:227316/1‑141 (MQ=11)

                               gTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCaacaac                                                                                                              >  3:451132/1‑141 (MQ=32)

                                ttAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACt                                                                                                             <  4:399155/141‑1 (MQ=32)

                                  aaaTATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACttt                                                                                                           <  2:75147/141‑1 (MQ=32)

                                      attatAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAg                                                                                                       <  1:881727/141‑1 (MQ=32)

                                      attatAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAg                                                                                                       <  1:77798/141‑1 (MQ=32)

                                         ataGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCAc                                                                                                    >  1:460349/1‑141 (MQ=32)

                                          taGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACg                                                                                                   <  2:636447/141‑1 (MQ=32)

                                            gAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGta                                                                                                 <  2:1026550/141‑1 (MQ=32)

                                              aaaCATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGtata                                                                                               <  4:9414/141‑1 (MQ=32)

                                                aCATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAat                                                                                             <  4:155391/141‑1 (MQ=32)

                                                 cATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAatg                                                                                            <  3:293012/141‑1 (MQ=32)

                                                   tCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAatgat                                                                                          >  2:747598/1‑141 (MQ=32)

                                                    cAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAatgatg                                                                                         >  1:383398/1‑141 (MQ=32)

                                                    cAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAAGGAATCAATTTAATATCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAatgatg                                                                                         >  1:102402/1‑141 (MQ=11)

                                                     aaaGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAAtgatga                                                                                        >  1:888076/1‑141 (MQ=32)

                                                      aaGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATgatgat                                                                                       <  3:97343/141‑1 (MQ=32)

                                                       aGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAGATTTCAACAACTTTAAAGCACGTATAATGATGAtt                                                                                      <  2:350724/141‑1 (MQ=17)

                                                            gTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAg                                                                                 >  4:258720/1‑141 (MQ=32)

                                                            gTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAg                                                                                 >  2:247491/1‑141 (MQ=32)

                                                             ttAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGc                                                                                <  4:233964/141‑1 (MQ=32)

                                                              tAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCt                                                                               <  2:460349/141‑1 (MQ=32)

                                                              tAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCt                                                                               <  3:3347/141‑1 (MQ=32)

                                                                aGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTg                                                                             <  4:203091/141‑1 (MQ=32)

                                                                  aaaTACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTa                                                                           <  3:452784/141‑1 (MQ=32)

                                                                    aTACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACa                                                                         <  2:486747/141‑1 (MQ=32)

                                                                         aTTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAatga                                                                                           <  1:461471/118‑1 (MQ=33)

                                                                         aTTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGga                                                                    <  1:987930/141‑1 (MQ=32)

                                                                          tttATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGgag                                                                   <  2:1082000/141‑1 (MQ=32)

                                                                          tttATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGgag                                                                   >  4:101423/1‑141 (MQ=32)

                                                                           ttattaCCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGgaga                                                                  >  1:57637/1‑141 (MQ=32)

                                                                           ttattaCCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGgaga                                                                  >  1:711/1‑141 (MQ=32)

                                                                               taCCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCATCTTGTACAAAGGAGaaaaa                                                              >  4:200827/1‑141 (MQ=17)

                                                                                    aGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAaagaa                                                         <  1:447131/141‑1 (MQ=35)

                                                                                          taataaTGGAGCAATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAAc                                                   <  2:856775/141‑1 (MQ=255)

                                                                                           aataatGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAAcc                                                  >  2:735314/1‑141 (MQ=38)

                                                                                           aataatGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAAcc                                                  >  1:509326/1‑141 (MQ=38)

                                                                                            ataatGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCa                                                 >  1:739304/1‑141 (MQ=38)

                                                                                             taatGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCaa                                                <  2:400721/141‑1 (MQ=37)

                                                                                              aatGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAg                                               <  2:419514/141‑1 (MQ=37)

                                                                                                 ggTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGccc                                            >  2:896474/1‑141 (MQ=38)

                                                                                                    gCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCaat                                         >  4:83338/1‑141 (MQ=38)

                                                                                                         aGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATgg                                    <  1:834350/141‑1 (MQ=38)

                                                                                                          gAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGa                                   <  4:50665/141‑1 (MQ=38)

                                                                                                          gAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGa                                   >  2:514150/1‑141 (MQ=38)

                                                                                                          gAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGa                                   >  1:690739/1‑141 (MQ=38)

                                                                                                           aaGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGAc                                  <  4:277087/141‑1 (MQ=37)

                                                                                                           aaGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGAc                                  >  1:185379/1‑141 (MQ=38)

                                                                                                            aGGAATTAATAATAAGATAAAAGTAATCAAAAGAATGGAATTTGGTTACAGAAATTTCAACAACTTTAAAGAAAGTATAATGATGATTTTCAGCTTGAAAAAAAAAAAAAAAAAGAAGACAACCAAGCCCAATAATGGACt                                 <  4:300387/141‑1 (MQ=11)

                                                                                                                aTTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTAGTACAAAGGAGAAAAAAAGAAGACAACCAAGCACAATACTGGACTGgca                             >  2:1024194/1‑140 (MQ=255)

                                                                                                                  taataataaGATAAAATTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGgccgc                           >  1:480932/1‑141 (MQ=9)

                                                                                                                   aataataaGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAAGGATGATTTGCAGCTTGTACAAATGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGgccgcc                          >  3:130463/1‑141 (MQ=255)

                                                                                                                    ataataaGATAAAATTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGATTGGCCGCCt                         <  1:816025/141‑1 (MQ=12)

                                                                                                                     taataaGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCta                        <  3:105352/141‑1 (MQ=38)

                                                                                                                      aataaGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaa                       <  1:64233/141‑1 (MQ=37)

                                                                                                                       ataaGATAAAATTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaat                      >  2:419402/1‑141 (MQ=9)

                                                                                                                        taaGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGAACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaata                     <  2:441217/141‑1 (MQ=38)

                                                                                                                          aGATAAAATTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaataat                   >  1:1005605/1‑141 (MQ=17)

                                                                                                                              aaaaTTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATaaaa               <  2:201443/141‑1 (MQ=37)

                                                                                                                               aaaTTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATGATAAAAg              <  1:858881/141‑1 (MQ=38)

                                                                                                                               aaaTTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaataa                    >  2:177376/1‑135 (MQ=9)

                                                                                                                               aaaTTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCtaataa                    <  1:177376/135‑1 (MQ=12)

                                                                                                                               aaaTTAATCAAGCGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTATACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATaaaaa              >  1:674568/1‑140 (MQ=2)

                                                                                                                                aaTTAATCCAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTTTCCGCCTAATAATAAAAGc             <  4:217741/141‑1 (MQ=255)

                                                                                                                                 atttaccAAGTGAATTTCTTTTGGTTACCCACATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGCGAAAAACAGAAGCCACCCCACCCCAATAATGGACTGGCCGCCTAATAATAAAAGct            <  2:270318/135‑1 (MQ=11)

                                                                                                                                 aTTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGct            <  3:78434/141‑1 (MQ=39)

                                                                                                                                 aTTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGct            >  3:255552/1‑141 (MQ=40)

                                                                                                                                 aTTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGct            <  2:21701/141‑1 (MQ=40)

                                                                                                                                    aaTCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTa         <  2:557344/141‑1 (MQ=40)

                                                                                                                                    aaTCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACCTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTa         <  1:301327/141‑1 (MQ=255)

                                                                                                                                     aTCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTaa        <  2:1070811/141‑1 (MQ=40)

                                                                                                                                     aTCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTaa        >  1:111452/1‑141 (MQ=40)

                                                                                                                                       cAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTaaaa      >  2:68875/1‑141 (MQ=40)

                                                                                                                                          gTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGtt   >  1:870118/1‑141 (MQ=39)

                                                                                                                                          gTGAAGTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGtt   <  1:513/141‑1 (MQ=255)

                                                                                                                                           tGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGTTg  >  1:35081/1‑141 (MQ=38)

                                                                                                                                           tGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGTTg  <  1:1056287/141‑1 (MQ=38)

                                                                                                                                           tGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGTTg  <  1:554746/141‑1 (MQ=38)

AAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTGGTTACAGAAATTTCAACAACTTTAAAGCACGTATAATGATGATTTTCAGCTTGTACAAAGGAGAAAAAAAGAAGACAACCAAGCCCAATAATGGACTGGCCGCCTAATAATAAAAGCTCTAAAAGTTG  >  USA300TCH1516_ALE/2291873‑2292152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`