BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,861	T→C	71.6%	intergenic (‑944/‑224)	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,861	0	T	C	71.6%	236.2 / 97.4	157	intergenic (‑944/‑224)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (18/27); new base C (61/51); total (79/78)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

ATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCT  >  CP000730/1929721‑1929994

aTTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgtgt                                                                                                                                             <  2:585006/141‑1 (MQ=255)

     gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                                      <  2:553588/141‑1 (MQ=255)

        tAGACGATTCAATCTGTCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑GAA‑G‑G‑Tg                                                                                                                                 >  1:333468/1‑141 (MQ=255)

        tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑Tg                                                                                                                                 >  2:280709/1‑141 (MQ=255)

         aGACGATGCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGa                                                                                                                                <  2:101703/141‑1 (MQ=255)

                tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAt                                                                                                                         <  1:510197/141‑1 (MQ=255)

                tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAt                                                                                                                         <  1:558996/141‑1 (MQ=255)

                  aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAt‑‑tt                                                                                                                     <  2:478985/141‑1 (MQ=255)

                      tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGt                                                                                                                 <  1:579224/141‑1 (MQ=255)

                        tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTc                                                                                                               <  2:604374/141‑1 (MQ=255)

                             tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCAcatc                                                                                                          <  2:307666/141‑1 (MQ=255)

                                aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTGTCTCTCAGGCGTCAATTTAGACGAACGAAGAAGCTTTG‑AAA‑G‑G‑TGATGCTTAT‑‑GTTGGGTCACATCATa                                                                                                       >  1:345156/1‑141 (MQ=255)

                                 aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAg                                                                                                      >  1:507327/1‑141 (MQ=255)

                                    atatAACTTATCATAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCAc                                                                                                   >  1:97926/1‑141 (MQ=255)

                                                 atcggaagtagatatgCTTGTATATTCACAGATGGTATAATCAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                     <  1:513887/127‑1 (MQ=255)

                                                 atcggaaggagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                     <  1:71181/127‑1 (MQ=255)

                                                 aacggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                     <  2:316223/127‑1 (MQ=255)

                                                  tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                    >  2:632488/14‑141 (MQ=255)

                                                   cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATGGCACCAGTCATCAGTGGCtg                                                                                   >  1:272690/13‑141 (MQ=255)

                                                   cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATGCCACCAGTCATCGGTGGCtg                                                                                   >  2:277398/13‑141 (MQ=255)

                                                   cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGGCA‑TAACT‑T‑TGCTAATCAT‑‑TCTCGTTCACATCATACCACCAGCCATAAGTGGCtg                                                                                   >  1:382966/13‑141 (MQ=255)

                                                       agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                               <  1:259322/133‑2 (MQ=255)

                                                        gtaggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCGGTCATCAGTGGCTGTGctg                                                                              <  1:354262/134‑3 (MQ=255)

                                                        gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                                                               >  2:358101/1‑141 (MQ=255)

                                                         tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTCTCTCTCATTCTTCAATT‑‑GACGAAGAGAGGAGGGGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                             <  1:288791/135‑4 (MQ=255)

                                                         tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                             <  1:26783/135‑4 (MQ=255)

                                                         tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                             <  2:629877/135‑4 (MQ=255)

                                                            cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                           <  1:408338/141‑1 (MQ=255)

                                                            atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGCCTGTGCTAGTGc                                                                          >  2:190326/4‑141 (MQ=255)

                                                             tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                                                                         <  1:205851/139‑8 (MQ=255)

                                                               tgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACACCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGt                                                                        <  1:421486/141‑1 (MQ=255)

                                                                 cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTTTGctgtttcttatt                                                                     >  2:292403/1‑130 (MQ=255)

                                                                 cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                     >  2:463506/1‑141 (MQ=255)

                                                                 cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                     <  2:550349/141‑4 (MQ=255)

                                                                  ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                    <  1:538881/141‑5 (MQ=255)

                                                                  ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                    <  2:356769/141‑5 (MQ=255)

                                                                  ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                    >  1:638863/1‑140 (MQ=255)

                                                                  ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TCCTAATCAT‑‑TTTGGTTCACATCATAGCACCAGTCTTCAGTGGCTGTGCTGTTGCGtatta                                                                    >  2:282518/1‑140 (MQ=255)

                                                                  ttGTATATTCACAGATGGTATAACAATAGAGTCTCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                    <  1:102810/141‑5 (MQ=255)

                                                                   tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAT‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                    <  1:424911/141‑1 (MQ=255)

                                                                   tGTATATTCACAGTTGGTATAACAATAGAGTCGCCTATATCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGGTGCATTGC‑T‑T‑TGCTAATCAT‑‑TTTCGATCACATCATAGCACCAGTCATCAGGGGCTATGctgtggggtattat                                                                   >  2:29800/1‑128 (MQ=255)

                                                                   tGTATATTCACAGCTGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                   <  1:278405/141‑6 (MQ=255)

                                                                    gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                  >  1:242208/12‑138 (MQ=255)

                                                                      atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                >  1:572382/10‑136 (MQ=255)

                                                                       tattcacagatggtATAATAATGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                               <  2:274578/133‑7 (MQ=255)

                                                                       tattatcagatggtATAATAATAGTGTTTGCTATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                               <  1:120571/133‑7 (MQ=255)

                                                                         ttcacagatggtATAATGATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                             >  1:469921/7‑141 (MQ=255)

                                                                          tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                            >  2:91718/6‑141 (MQ=255)

                                                                           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTCTCCTTCTTg                                                           >  1:458604/5‑141 (MQ=255)

                                                                           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                           >  2:609170/5‑141 (MQ=255)

                                                                           cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTg                                                           >  2:534203/5‑141 (MQ=255)

                                                                            acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                          >  2:456364/4‑141 (MQ=255)

                                                                             ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                          <  2:161790/141‑1 (MQ=255)

                                                                             ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                          <  1:448480/141‑1 (MQ=255)

                                                                             caggtggtataATAATAGAGTCGCCTATCTCTCATGCGTCAATT‑‑GACGGCGAGAGGGGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTTCGTATTATACTTATTGGc                                                         <  2:285528/134‑1 (MQ=255)

                                                                              agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                        <  1:437605/140‑1 (MQ=255)

                                                                               tatggtATAATAATGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCTTAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                       <  2:320702/140‑2 (MQ=255)

                                                                               gattgaataacaATATAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTACTCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                       <  2:588350/135‑2 (MQ=255)

                                                                               gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                       >  2:508681/1‑140 (MQ=255)

                                                                                 tgttataATAATAGAGTCGTCTATCTCTCAGGCGTCAA‑‑TGGACGAATAGAGGAGGTCCATTAC‑T‑T‑TGCTAATCCT‑‑TTTCGTTCGCATCATAGCACCAGTCAGCAGTGGCTGTGCTGTGGCGTATTATACTTATTGGCTTAg                                                     <  2:270883/138‑1 (MQ=255)

                                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                     <  2:415240/138‑4 (MQ=255)

                                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                     <  1:612933/138‑4 (MQ=255)

                                                                                 tggtataATAATAGAGTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                     <  1:446155/138‑4 (MQ=255)

                                                                                 gagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                     <  1:638521/138‑4 (MQ=255)

                                                                                  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCttggt                                                    >  1:573848/3‑137 (MQ=255)

                                                                                  gatATAATAGTAGAGCCGCCTATCTCTCAGGCGACAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTAGTTCACACCATAGAACCAGGCAACAGTGGCTGTGGCATTGCGTGTTTTTcatcttcactaat                                                     >  2:625437/1‑129 (MQ=255)

                                                                                   gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                   >  1:496047/2‑141 (MQ=255)

                                                                                   gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCTGTTGCGTATTATACTGATTGGCTTAGTa                                                   >  1:74474/2‑141 (MQ=255)

                                                                                   atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                    >  1:652215/1‑141 (MQ=255)

                                                                                     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                 >  2:493284/1‑139 (MQ=255)

                                                                                     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑‑TGCTAATCATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                >  2:388492/1‑141 (MQ=255)

                                                                                      aaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                <  1:99481/140‑4 (MQ=255)

                                                                                       aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaaag                                               >  2:647690/1‑137 (MQ=255)

                                                                                       aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                               >  2:111394/1‑141 (MQ=255)

                                                                                         taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                             >  1:566730/4‑141 (MQ=255)

                                                                                           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                           >  2:498300/2‑140 (MQ=255)

                                                                                            tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                           <  1:407826/141‑1 (MQ=255)

                                                                                             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                         >  2:297534/1‑138 (MQ=255)

                                                                                              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                        <  1:617180/141‑5 (MQ=255)

                                                                                                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                      <  1:14081/141‑7 (MQ=255)

                                                                                                 tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                     <  1:618995/141‑8 (MQ=255)

                                                                                                 tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                     >  2:500234/1‑134 (MQ=255)

                                                                                                 tCGCCTATCTCCCAGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                     <  1:8777/141‑8 (MQ=255)

                                                                                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                    >  1:261191/1‑133 (MQ=255)

                                                                                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                    >  2:379385/1‑133 (MQ=255)

                                                                                                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                    >  1:358919/1‑133 (MQ=255)

                                                                                                    ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                  >  2:173791/1‑131 (MQ=255)

                                                                                                    ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                  >  1:81893/1‑131 (MQ=255)

                                                                                                    ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                  >  2:641759/1‑131 (MQ=255)

                                                                                                    ccTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                  <  1:602673/141‑11 (MQ=255)

                                                                                                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                 >  1:144924/1‑130 (MQ=255)

                                                                                                      tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatttag                                >  2:164795/1‑129 (MQ=255)

                                                                                                      tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                >  2:13627/1‑129 (MQ=255)

                                                                                                      tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                >  1:93114/1‑129 (MQ=255)

                                                                                                       aTCTCTCAGGCGTCAATT‑‑GCCGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatg                               >  1:172730/1‑128 (MQ=255)

                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                               >  1:168804/1‑128 (MQ=255)

                                                                                                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                               >  2:577827/1‑128 (MQ=255)

                                                                                                                gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                       >  2:526166/1‑141 (MQ=255)

                                                                                                                    cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                   >  1:464158/1‑141 (MQ=255)

                                                                                                                                     ggaggTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt  >  2:26504/1‑141 (MQ=255)

ATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCT  >  CP000730/1929721‑1929994

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: `‑`