BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,784	G→C	7.8%	M152I (ATG→ATC) ‡	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,784	0	G	C	7.8%	574.0 / 30.2	180	M152I (ATG→ATC) ‡	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base G (87/79); new base C (4/10); total (91/89)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 183 total aligned reads displayed.

TCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGC  >  NC_002947/549636‑549930

  cTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGaa                                                                                                                                                 <  2:354782/150‑1 (MQ=39)

  cTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGaa                                                                                                                                                 <  1:4004716/150‑1 (MQ=39)

     gTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGtt                                                                                                                                              <  1:3229547/150‑1 (MQ=39)

      tttGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTc                                                                                                                                             <  1:1767351/150‑1 (MQ=39)

        tGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                                                           <  1:3583811/150‑1 (MQ=39)

        gtttCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTcgc                                                                                                                                           <  1:1604452/147‑1 (MQ=255)

         gCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgcg                                                                                                                                          <  2:3253137/150‑1 (MQ=39)

          cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                                                                                                                                         <  2:2085160/150‑1 (MQ=39)

             ggCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACct                                                                                                                                      >  1:1666801/1‑150 (MQ=39)

               ccgccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgc                                                                                                                                    >  1:2142348/1‑150 (MQ=39)

                   cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                                                                >  1:2884663/1‑150 (MQ=39)

                   cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                                                                >  1:194483/1‑150 (MQ=39)

                     aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa                                                                                                                              <  1:591708/150‑1 (MQ=255)

                      tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                                                             <  2:385246/150‑1 (MQ=39)

                       ggcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTcgaag                                                                                                                            >  2:5028980/4‑146 (MQ=255)

                           cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACg                                                                                                                        >  1:914205/1‑150 (MQ=39)

                            gATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGa                                                                                                                       <  2:2301326/150‑1 (MQ=39)

                             aTGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGAc                                                                                                                      <  1:4097632/150‑1 (MQ=39)

                                  gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                 <  2:2006377/150‑1 (MQ=38)

                                  gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                 >  2:1557656/1‑150 (MQ=38)

                                  acagaCCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTc                                                                                                                                <  2:2384525/131‑1 (MQ=255)

                                   cagaCCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCa                                                                                                                <  1:291745/147‑1 (MQ=255)

                                           tGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacga                                                                                                        <  1:3175683/150‑1 (MQ=2)

                                               cacaTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACt                                                                                                    >  1:432520/1‑150 (MQ=37)

                                                 caTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcc                                                                                                  >  1:1283129/1‑144 (MQ=2)

                                                  aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                 >  1:2340943/1‑150 (MQ=37)

                                                  aTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcca                                                                                                 >  1:738686/1‑143 (MQ=2)

                                                   tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                >  2:3089821/1‑150 (MQ=37)

                                                          gTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcat                                                                                         <  1:2672153/150‑1 (MQ=34)

                                                          cgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatcaa                                                                                         <  2:3112262/146‑2 (MQ=2)

                                                              cGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGt                                                                                     >  1:4881103/1‑150 (MQ=34)

                                                               gTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGtt                                                                                    <  2:3514294/150‑1 (MQ=34)

                                                                tCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTc                                                                                   <  2:1487984/150‑1 (MQ=34)

                                                                    gTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCt                                                                               >  2:518861/1‑150 (MQ=34)

                                                                     tAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTc                                                                              >  2:2058641/1‑150 (MQ=34)

                                                                          gTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATg                                                                         <  1:5053137/150‑1 (MQ=34)

                                                                           ttCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCATATgg                                                                        >  1:2078410/1‑150 (MQ=21)

                                                                               gTACATCGTGGTCGTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                                                                    <  2:2180445/150‑1 (MQ=37)

                                                                                 aCATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTgg                                                                  >  2:196830/1‑150 (MQ=37)

                                                                                  cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                                                                 <  1:1097917/150‑1 (MQ=37)

                                                                                  cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                                                                 <  1:5299957/150‑1 (MQ=37)

                                                                                   aTCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGaa                                                                >  1:652031/1‑150 (MQ=37)

                                                                                    tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                               <  2:2142348/150‑1 (MQ=37)

                                                                                        ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                                                           <  2:1188997/150‑1 (MQ=37)

                                                                                          tCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCCTCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAg                                                         >  1:4487655/1‑150 (MQ=37)

                                                                                           cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                        <  2:3158865/150‑1 (MQ=37)

                                                                                           cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                        <  1:2562115/150‑1 (MQ=37)

                                                                                           cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                        <  2:430204/150‑1 (MQ=37)

                                                                                            ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                >  1:1241140/1‑109 (MQ=35)

                                                                                            ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                <  2:1241140/109‑1 (MQ=35)

                                                                                            ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgac                                                       <  2:786427/150‑1 (MQ=37)

                                                                                              ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     >  2:5429590/1‑150 (MQ=37)

                                                                                              ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                     <  2:2678890/150‑1 (MQ=37)

                                                                                               cTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacgac                                                    <  2:1337406/150‑1 (MQ=37)

                                                                                                tGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACa                                                   >  1:1626727/1‑150 (MQ=37)

                                                                                                 gatcAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                                                           >  2:4277122/4‑61 (MQ=255)

                                                                                                 gAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                                                           <  1:4277122/61‑1 (MQ=255)

                                                                                                  aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc                                                 >  1:4101877/1‑150 (MQ=37)

                                                                                                    cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                               >  1:4601671/1‑150 (MQ=37)

                                                                                                     aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                                                         >  1:2509401/1‑75 (MQ=39)

                                                                                                     aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                                                         <  2:2509401/75‑1 (MQ=39)

                                                                                                     aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGaa                                              >  2:3884077/1‑150 (MQ=37)

                                                                                                       ggCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAAt                                            <  1:4802643/150‑1 (MQ=37)

                                                                                                         ccGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                          <  1:1206205/150‑1 (MQ=35)

                                                                                                          cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                         <  1:1363742/150‑1 (MQ=35)

                                                                                                          cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                         <  2:3577833/150‑1 (MQ=35)

                                                                                                          cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                         >  2:1451477/1‑150 (MQ=35)

                                                                                                          cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                         >  1:1632586/1‑150 (MQ=35)

                                                                                                            aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                               <  1:3421160/126‑1 (MQ=35)

                                                                                                            aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                               >  2:3421201/1‑126 (MQ=35)

                                                                                                                ggTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTAc                                   >  2:4528742/1‑150 (MQ=34)

                                                                                                                ggTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTAc                                   >  1:218997/1‑150 (MQ=34)

                                                                                                                    gacgacGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCt                               <  2:1100095/150‑1 (MQ=34)

                                                                                                                      cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                             <  2:3380190/150‑1 (MQ=34)

                                                                                                                      cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                             <  1:2084572/150‑1 (MQ=34)

                                                                                                                       gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTc                                                                                                   >  2:4633362/1‑79 (MQ=34)

                                                                                                                       gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatc                                                                                           <  2:3683676/87‑1 (MQ=33)

                                                                                                                       gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatc                                                                                           >  1:3683676/1‑87 (MQ=33)

                                                                                                                       gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                          >  2:2991824/1‑136 (MQ=34)

                                                                                                                       gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTc                            >  1:413694/1‑150 (MQ=35)

                                                                                                                          gCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaag                         >  2:4394230/1‑150 (MQ=37)

                                                                                                                            tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCTGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       >  2:1323727/1‑150 (MQ=37)

                                                                                                                            tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       >  1:1122967/1‑150 (MQ=37)

                                                                                                                                ctgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                          >  1:882663/1‑127 (MQ=35)

                                                                                                                                  gctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGt                 >  2:170624/1‑150 (MQ=37)

                                                                                                                                   ctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                         <  1:2280471/125‑1 (MQ=35)

                                                                                                                                    tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       <  2:1013607/142‑1 (MQ=37)

                                                                                                                                    tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       >  1:1013607/1‑142 (MQ=37)

                                                                                                                                    tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       <  1:4968286/142‑1 (MQ=37)

                                                                                                                                    tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                       >  2:4968286/1‑142 (MQ=37)

                                                                                                                                     ggAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGa              >  1:712814/1‑150 (MQ=37)

                                                                                                                                      gAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaa             <  1:3686878/150‑1 (MQ=37)

                                                                                                                                       aaCTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaaa            <  2:2755705/150‑1 (MQ=35)

                                                                                                                                            gTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTg       <  2:652031/150‑1 (MQ=34)

                                                                                                                                             tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg      >  1:794577/1‑150 (MQ=34)

                                                                                                                                             tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg      >  1:1472407/1‑150 (MQ=34)

                                                                                                                                               gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt    >  2:2564757/1‑150 (MQ=35)

                                                                                                                                               gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt    >  1:4742711/1‑150 (MQ=35)

                                                                                                                                                 gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc  >  1:382861/1‑150 (MQ=37)

                                                                                                                                                  aTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                             <  2:851410/122‑1 (MQ=32)

TCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGC  >  NC_002947/549636‑549930

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