BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,848	C→A	67.4%	intergenic (‑931/‑237)	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,848	0	C	A	67.4%	150.2 / 157.2	166	intergenic (‑931/‑237)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (18/36); new base A (62/49); total (80/86)
Fisher's exact test for biased strand distribution p-value = 7.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-02

Only 100 of 166 total aligned reads displayed.

AATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACC  >  CP000730/1929711‑1929980

aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                                                           <  2:372449/141‑1 (MQ=255)

 aTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAAGCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag                                                                                                                                          <  1:115892/141‑1 (MQ=255)

   ttGTCCAAATAAGTGGAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                                                        <  2:601120/141‑1 (MQ=255)

   ttGGCCAATTAAGTGGAGACGATTCAATCTATCAGCTAGAATTATATAACTTATCAGAAGCAAGTGACTGTTCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTACCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                                                        <  1:321033/141‑1 (MQ=255)

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

                tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g‑g                                                                                                                        <  2:500245/141‑1 (MQ=255)

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

                          tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTT‑At                                                                                                           <  2:100526/141‑1 (MQ=255)

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

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

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

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

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

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

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

                                                           gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTgg                                                                        >  1:139232/1‑141 (MQ=255)

                                                           atcggaagtatatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCGCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑GTTCGTTCACATCATAGCACCAGTCACCAGTGGc                                                                       <  2:317311/127‑1 (MQ=255)

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

                                                           atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                       >  1:260077/15‑141 (MQ=255)

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

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

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

                                                              ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                    >  2:277141/12‑141 (MQ=255)

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

                                                                 agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                 <  1:195437/133‑2 (MQ=255)

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

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

                                                                   tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtgatgt                                                               <  1:478515/135‑5 (MQ=255)

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

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

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

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

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

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

                                                                        atgCTTGTATATTCACAGATGGCATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                          <  2:653165/140‑1 (MQ=255)

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

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

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

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

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

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

                                                                             tgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                     >  2:624099/13‑139 (MQ=255)

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

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

                                                                                 aatgcacagatggtAGAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TGGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                 <  1:594747/133‑7 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                         gatggtATAATAATAGAGTCGCCTATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                         <  2:536543/141‑2 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                      tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                            <  1:292403/141‑3 (MQ=255)

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

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

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

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

                                                                                                       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TGCGAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                           >  2:629138/1‑138 (MQ=255)

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

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

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

                                                                                                           tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                       >  1:162660/1‑134 (MQ=255)

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

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

                                                                                                              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                    <  1:97113/141‑11 (MQ=255)

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

                                                                                                                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                  >  1:508557/1‑129 (MQ=255)

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

                                                                                                                 ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAT‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACGTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                 <  1:450432/140‑14 (MQ=255)

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

                                                                                                                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                 >  2:580444/1‑128 (MQ=255)

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

                                                                                                                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                 >  1:568979/1‑128 (MQ=255)

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

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

                                                                                                                                 tttGGACGCAGACAGGAGGTGTA‑TCA‑G‑G‑T‑GACGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc  <  2:63656/141‑1 (MQ=255)

AATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACC  >  CP000730/1929711‑1929980

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