BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,870	A→C	72.9%	intergenic (‑953/‑215)	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,870	0	A	C	72.9%	269.0 / 90.4	154	intergenic (‑953/‑215)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (16/26); new base C (59/51); total (77/77)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

AGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGT  >  CP000730/1929730‑1930010

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)

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

                        aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCCTAg                                                                                                                      >  2:349720/1‑141 (MQ=255)

                          cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTACGTTGGCATTATAGCa                                                                                                                    <  2:146991/141‑1 (MQ=255)

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

                           atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCAc                                                                                                                   >  1:400338/1‑141 (MQ=255)

                                    aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAC‑‑TTTCGCTCACAGCATAGCACCAGtcatca                                                                                                          <  2:472377/141‑1 (MQ=255)

                                        atcggaagtatatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCGCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑GTTCGTTCACATCATAGCACCAGTCACCAGTGGc                                                                                                     <  2:317311/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)

                                           ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                  <  1:408990/130‑1 (MQ=255)

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

                                               gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                              >  2:495662/8‑139 (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‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtgatgt                                                                                             <  1:478515/135‑5 (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)

                                                 gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTGT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                             >  1:553498/1‑141 (MQ=255)

                                                  aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                            <  1:87817/141‑1 (MQ=255)

                                                   cTGTGCTTGTATTTTTGCCGATGAAATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                                           <  1:626350/141‑1 (MQ=255)

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

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

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

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

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

                                                          tgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                   >  2:624099/13‑139 (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)

                                                           gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGttttata                                                                                  >  1:575104/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)

                                                                ttcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAAAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                             >  1:429667/7‑141 (MQ=255)

                                                                 tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                            >  2:91718/6‑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                                                                          <  1:448480/141‑1 (MQ=255)

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

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

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

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

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

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

                                                                        tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                     >  1:649279/4‑141 (MQ=255)

                                                                        gagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                     <  1:638521/138‑4 (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:352569/1‑141 (MQ=255)

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

                                                                           tataATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                  >  2:280738/1‑140 (MQ=255)

                                                                            ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                 <  1:445057/141‑3 (MQ=255)

                                                                            ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                 >  1:354903/1‑139 (MQ=255)

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

                                                                             taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                <  2:29742/141‑4 (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)

                                                                                taatagAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                             <  2:376500/138‑1 (MQ=255)

                                                                                 agtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                             <  1:206125/141‑1 (MQ=255)

                                                                                 aatagGGTCGCCTATATCTCAGGCGTCAATT‑‑GACGAACAGAGGCGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAATGCaa                                                            <  1:469100/139‑1 (MQ=255)

                                                                                   tagAGTCGCCTATCCCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                           <  1:53886/141‑1 (MQ=255)

                                                                                    agctgctccTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                         <  1:77642/134‑4 (MQ=255)

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

                                                                                    agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                         <  1:277141/141‑4 (MQ=255)

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

                                                                                       gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                       <  2:400338/141‑1 (MQ=255)

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

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

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

                                                                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGCGATGAGGGTGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaacaatt                                                   <  2:300611/141‑10 (MQ=255)

                                                                                           ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                  >  2:173791/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‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                >  1:508557/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)

                                                                                              ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAT‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACGTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                               <  1:450432/140‑14 (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‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                               >  2:580444/1‑128 (MQ=255)

                                                                                              aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                               >  1:568979/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)

                                                                                                            aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                  <  2:454536/141‑1 (MQ=255)

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

                                                                                                                                     ttaA‑G‑G‑TGATGCGTCT‑‑TTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg          <  1:86504/140‑1 (MQ=255)

                                                                                                                                          g‑g‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa      <  1:280709/141‑1 (MQ=255)

                                                                                                                                                aTGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt  >  2:190265/1‑141 (MQ=255)

AGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGT  >  CP000730/1929730‑1930010

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