BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,864	A→T	14.8%	intergenic (‑947/‑221)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,864	0	A	T	14.8%	552.8 / 66.2	212	intergenic (‑947/‑221)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (88/90); new base T (14/17); total (102/110)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 215 total aligned reads displayed.

GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929728‑1930004

gtAGACGATTCAATCTATCAGCGAGAATCATATAACTTACAAGAAGCAAGTGACTGTGCTTGTAAATTTGCCGATGGTATAATAGTAGAGGCGGCTAACTATCAGGCGTCAATT‑TAGACGGAGAGAGGAGGTGtg‑taggat                                                                                                                                              >  1:102090/1‑137 (MQ=255)

 tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑Tg                                                                                                                                          >  1:638417/1‑141 (MQ=255)

    aCGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑Tg                                                                                                                                      <  2:1003221/141‑1 (MQ=255)

       aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGAGTA‑TAA‑G‑G‑TGA‑TGCtt                                                                                                                                   >  1:188102/1‑141 (MQ=255)

         tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTG‑TAA‑G‑G‑TGA‑TGCTTAt                                                                                                                                 >  2:774601/1‑141 (MQ=255)

          cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAtt                                                                                                                                >  1:851058/1‑141 (MQ=255)

                 tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTc                                                                                                                         <  2:29043/141‑1 (MQ=255)

                 tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTc                                                                                                                         <  2:442753/141‑1 (MQ=255)

                    gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTcaca                                                                                                                      <  1:207966/141‑1 (MQ=255)

                     cTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCAcat                                                                                                                     <  1:106887/141‑1 (MQ=255)

                         aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATa                                                                                                                 <  2:447528/141‑1 (MQ=255)

                         aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATa                                                                                                                 >  2:1042748/1‑141 (MQ=255)

                                 aaCTTATGAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGt                                                                                                         <  2:405118/141‑1 (MQ=255)

                                 aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCGGGGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGt                                                                                                         <  1:201705/141‑1 (MQ=255)

                                 aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGt                                                                                                         <  2:236841/141‑1 (MQ=255)

                                   cTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTTTA‑TCA‑G‑G‑TGA‑TGCTTCTCTTCGTTCACATCATAGCACCAGtca                                                                                                       <  2:1023348/141‑1 (MQ=255)

                                       tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                                                                                   >  1:905344/1‑141 (MQ=255)

                                          atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                               >  2:419563/15‑141 (MQ=255)

                                           tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                              >  1:530107/14‑141 (MQ=255)

                                             ttaactagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                            <  2:558935/130‑1 (MQ=255)

                                             ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                            <  1:985514/130‑1 (MQ=255)

                                             gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                             >  1:1023864/1‑141 (MQ=255)

                                             gCAAGTGACGGGTCTTGGAGATTTGCCGATGATATAAGAGTAGAGGCGCCTATCTCTCAGGCGTCAACT‑TAGAAGCAGAGAGGAGGGGTA‑TAA‑G‑G‑GGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGcgg                                                                                             <  2:308771/141‑3 (MQ=255)

                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:169492/131‑1 (MQ=255)

                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:336630/131‑1 (MQ=255)

                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:603360/131‑1 (MQ=255)

                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:26177/131‑1 (MQ=255)

                                               aagtagatatgattGTATATTCACAGAGGATATAAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGAAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                          >  2:970568/13‑141 (MQ=255)

                                                agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                         <  1:806288/133‑2 (MQ=255)

                                                 gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                        >  2:877495/8‑139 (MQ=255)

                                                  gagatatgCTTGTATATTCACAGATGGTATAATCATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                       <  1:26320/135‑4 (MQ=255)

                                                   gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       <  1:278438/141‑1 (MQ=255)

                                                   gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       >  2:399780/1‑141 (MQ=255)

                                                   gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑CAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       <  1:714266/141‑1 (MQ=255)

                                                   agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                      >  1:767931/6‑137 (MQ=255)

                                                   agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                      >  1:84295/6‑137 (MQ=255)

                                                     atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                    <  1:344168/138‑7 (MQ=255)

                                                      tatgCTTGTATATTCACAGATGGTATAATAATAGAGAGGCCGATGTCTCGGTCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCg                                                                                   >  2:175497/3‑141 (MQ=255)

                                                          cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtgtt                                                                               >  1:322164/1‑141 (MQ=255)

                                                           ttgttgatTTGCCGATGAGATAAGAGTAGAGTCTCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                               <  2:954764/135‑1 (MQ=255)

                                                           ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCATGCGTCAATT‑TAGCCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                               <  2:866725/141‑1 (MQ=255)

                                                           ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttt                                                                               <  2:974418/141‑1 (MQ=255)

                                                            cgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCCTCAGTGGCTGTGCTGTTGCGtattat                                                                             >  2:536012/13‑139 (MQ=255)

                                                             gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                            <  2:782078/130‑4 (MQ=255)

                                                             gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                            <  1:486249/141‑7 (MQ=255)

                                                               atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATAAGTGGCTGTGCTGTTGCGtattatact                                                                          >  1:214158/10‑136 (MQ=255)

                                                                 attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                        >  2:513247/8‑140 (MQ=255)

                                                                  ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                       >  1:495560/7‑141 (MQ=255)

                                                                     caagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                    <  1:863092/138‑1 (MQ=255)

                                                                     acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATTGAAGAGGCATGGACG‑GGCAT‑‑TAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                    <  1:823026/138‑1 (MQ=255)

                                                                       agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGCGGAGCTTCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                  <  1:1040905/140‑1 (MQ=255)

                                                                         atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                >  1:64641/5‑139 (MQ=255)

                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                               <  1:175163/138‑4 (MQ=255)

                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑T‑TGC‑TAATCATTTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                               >  2:583068/4‑141 (MQ=255)

                                                                           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                              <  2:186077/139‑1 (MQ=255)

                                                                           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                              <  1:655840/139‑1 (MQ=255)

                                                                           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                              >  1:1049589/3‑141 (MQ=255)

                                                                           ggtataAAAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                              <  2:614834/139‑1 (MQ=255)

                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTa                                                             >  1:59171/2‑141 (MQ=255)

                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                             >  2:18345/2‑141 (MQ=255)

                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                             <  2:539936/140‑1 (MQ=255)

                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGAAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                             <  1:200143/140‑1 (MQ=255)

                                                                            atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCCGTGAGGGGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                              <  1:767544/141‑1 (MQ=255)

                                                                            agataataGTAGAGTCGCCTATCTCGCAGGCGACAATT‑TAGACGCAGAGAGGAGGGGTA‑TTA‑G‑G‑TGA‑TGCTTATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                              <  1:395164/139‑1 (MQ=255)

                                                                              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                           >  2:144909/1‑139 (MQ=255)

                                                                              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTATTATACTTATTGGCTTAGTaaa                                                           <  2:911961/141‑3 (MQ=255)

                                                                               taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                          <  2:810339/141‑4 (MQ=255)

                                                                                 ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                        <  2:944018/141‑1 (MQ=255)

                                                                                 ataagagagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                        <  2:178958/136‑1 (MQ=255)

                                                                                     tagaggcctcTATCTCTCAGGTGTCAATT‑‑‑GACGAAGTGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTCTTGGCTTAGTAAACGCaaca                                                    <  1:746053/132‑3 (MQ=255)

                                                                                     tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                    >  2:737965/1‑139 (MQ=255)

                                                                                      aaagTCGCCTATCTCTCAGGCAGCAA‑‑‑TGGACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                   <  2:32700/139‑4 (MQ=255)

                                                                                       gagTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                  >  1:218811/1‑137 (MQ=255)

                                                                                        agTCGCCTATCTCTCAGGCGGCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATCGTAGCACCAGGCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAAcacat                                                 >  1:610664/1‑140 (MQ=255)

                                                                                         gTCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                >  1:40418/1‑135 (MQ=255)

                                                                                          tCGCCTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                               <  1:901024/141‑8 (MQ=255)

                                                                                           cGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                               >  1:138422/1‑141 (MQ=255)

                                                                                            gCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                              >  1:598935/1‑141 (MQ=255)

                                                                                             ccTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                            <  2:985247/141‑11 (MQ=255)

                                                                                              cTATCTCTCAGGCGTCAATT‑‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                           >  1:551886/1‑130 (MQ=255)

                                                                                              cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATCGCACCAGTCCTCAGTGGCTGTGCTGTTTCGTATTAGACTTATTGGCTTAGTAAACGCacgaagtaattat                                           >  2:298574/1‑129 (MQ=255)

                                                                                              cTATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                           >  1:975657/1‑130 (MQ=255)

                                                                                               tATCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGGATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                          >  1:438285/1‑129 (MQ=255)

                                                                                                aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                          >  2:846315/1‑141 (MQ=255)

                                                                                                aTCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGTTAAGTAGACGCACTATAAAATAggg                                          >  2:958987/1‑140 (MQ=255)

                                                                                                aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                         >  2:695494/1‑128 (MQ=255)

                                                                                                aTCTCTCAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTAAACGCaacaaataattaata                                         >  2:89379/1‑128 (MQ=255)

                                                                                                 tctctcAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGTCTATTATACTTATTGGCTTAGTAAACGCaacaaataattaagag                                        <  2:437983/141‑15 (MQ=255)

                                                                                                 tctctcAGGCGTCAATT‑‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                        <  2:1062047/141‑15 (MQ=255)

                                                                                                  ctttcAGGCGTCAATT‑TAGACGCAGAGAGTAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGa                                        <  1:286633/138‑1 (MQ=255)

                                                                                                   tctcAGGCGTCAATT‑TACACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGCGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGAc                                       <  2:709455/141‑1 (MQ=255)

                                                                                                       aGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATCCAAAAAAGGTGCCatat                                   >  1:76182/1‑141 (MQ=255)

                                                                                                              aaTT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                            >  2:925655/1‑141 (MQ=255)

                                                                                                                tt‑tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCCTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                          <  2:836007/141‑1 (MQ=255)

                                                                                                                    aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                       <  1:794324/141‑1 (MQ=255)

                                                                                                                     gatgCAGAGTGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATa                      <  2:672060/138‑1 (MQ=255)

                                                                                                                        gCAGGGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                   <  1:345224/141‑1 (MQ=255)

                                                                                                                          agagagGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGCCATATAGCAGCCCCAATAAAAAt                 >  1:132380/1‑141 (MQ=255)

                                                                                                                             gaggagGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc              <  1:699398/141‑1 (MQ=255)

                                                                                                                                          aa‑G‑G‑TGG‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa  >  2:350452/1‑141 (MQ=255)

GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATT‑TAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929728‑1930004

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