BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	207,486	G→T	15.8%	M787I (ATG→ATT)	dnaE →	DNA polymerase III alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	207,486	0	G	T	15.8%	270.7 / 39.0	108	M787I (ATG→ATT)	dnaE	DNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand): ref base G (48/42); new base T (6/11); total (55/53)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 109 total aligned reads displayed.

GGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTT  >  NC_000913/207388‑207585

ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGa                                                                                                    >  1:2845684/1‑100 (MQ=255)

ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGa                                                                                                    <  2:1905309/100‑1 (MQ=255)

 gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                   <  1:195096/100‑1 (MQ=255)

 gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                   >  1:862418/1‑100 (MQ=255)

 gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                   >  1:2723022/1‑100 (MQ=255)

      aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgccg                                                                                              >  1:477556/1‑100 (MQ=255)

      aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgcc                                                                                               >  1:1255230/1‑99 (MQ=255)

       aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGCCTGAAAGCGCACTATCCTGCGGAGTTTATAACGGCGGTAATGACCACCGa                                                                                             >  1:1232853/1‑100 (MQ=255)

        cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGTGGTCATTAccgccgcc                                                                                            <  2:1041168/100‑3 (MQ=255)

        cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                            <  1:2031048/100‑3 (MQ=255)

        cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                            <  1:2491973/100‑3 (MQ=255)

        cAAATAGCACTCTGCGGCCTATGGTGTAATGTCATAGCAAACGTTATTGCTTAAAGTGCATTACCCTGCGGAGTTTATGGCGGCGGTCATTAccgacgca                                                                                            <  2:2035386/100‑3 (MQ=255)

         aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                           <  1:1177579/100‑4 (MQ=255)

         aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                            <  2:2153007/99‑1 (MQ=255)

          aaTCGCACTCTGCGGCCTATGTTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                            <  1:1809262/98‑3 (MQ=255)

          aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                          <  2:423825/100‑5 (MQ=255)

          aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGata                                                                                           <  2:2524924/99‑1 (MQ=255)

           aTCGTACTCTGCGGCCTATGCTTTTTTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgtcgccata                                                                                         <  2:1637916/100‑6 (MQ=255)

           aTCGCACTCTGGGGTCTATGCTTTGTTGTCATATCAAACGTTATGGCTGAAAGCGCATTATCCTGCTGATTTTATGGCGGCGGTAATGAccgccg                                                                                              <  2:2569921/95‑1 (MQ=255)

           aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                         >  2:2404413/1‑100 (MQ=255)

           aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                         <  1:6145/100‑1 (MQ=255)

           aTCGCACCCTGGGGCCTATGCTTTTGTGTCATATCAAACGTTATGGCTGAAAGCGCACTACCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                          <  2:2945883/99‑5 (MQ=255)

            tCGCACTCTGCTGCTTATTCTTTGGTTGCATATCAAACGTTATGGCTGAAAGCGCACTATCCAGCGGAGTTGATGGCGGCGGTAATGACCGCCGATATgg                                                                                        <  1:1882381/100‑1 (MQ=255)

            tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                        <  1:2601359/100‑1 (MQ=255)

             cGCACTCTGCGGCCTATGCCTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                       <  1:2678367/100‑1 (MQ=255)

              gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                        <  2:1355795/98‑1 (MQ=255)

              gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATAGCGGCGGTAATGACCGCCGATATGGac                                                                                      >  1:1329739/1‑100 (MQ=255)

               cACTCTGCGGCCTATGCTTTGGTGTTATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                     <  1:663434/100‑1 (MQ=255)

                aCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccataaactc                                                                                    <  2:2191897/100‑11 (MQ=38)

                 ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                                   >  1:1530242/1‑100 (MQ=255)

                 ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                                   <  1:2421591/100‑1 (MQ=255)

                   ctGCGGCTTATGCTTTGGTGTCATATCAAACTTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccataaac                                                                                      <  1:1176104/95‑9 (MQ=37)

                   ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                                 >  1:1946425/1‑100 (MQ=255)

                   ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                                 >  1:2166904/1‑100 (MQ=255)

                     gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCg                                                                               <  1:426224/100‑1 (MQ=255)

                        gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCAGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgag                                                                             >  2:2490885/1‑99 (MQ=255)

                         ccTATGCTTTGGTGTCGTATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGaa                                                                           >  1:1006159/1‑100 (MQ=255)

                            aTGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCTGAGTTTATGGCGGCGGAAATGACCGCCGATATGGACAACACCGAGAAgg                                                                         >  1:2694658/1‑99 (MQ=255)

                                 ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTGATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                   <  2:1610412/100‑1 (MQ=255)

                                 ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                   <  2:1277056/100‑1 (MQ=255)

                                  tGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                                  <  2:2166881/100‑1 (MQ=255)

                                   ggTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATAACCGCCGATATGGACAACAACGAGAAGGTGGTGGGt                                                                 >  1:1275797/1‑100 (MQ=255)

                                    gtgtGATATCAAGCTTCATGGCTGAAAGCGCACTATCGTGCTGAGTTTAGGGCGGCGTTAGTGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                                  <  2:2824462/98‑1 (MQ=255)

                                     tgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCt                                                               >  1:2620430/1‑100 (MQ=255)

                                      gtCATATGAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                              >  1:1593817/1‑100 (MQ=255)

                                      gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGAAGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                              >  1:1752399/1‑100 (MQ=255)

                                      gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGAGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                              >  1:1038663/1‑100 (MQ=255)

                                        cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggt                                                            >  2:271561/1‑100 (MQ=255)

                                          tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                          >  1:1544814/1‑100 (MQ=255)

                                          tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                          <  2:103431/100‑1 (MQ=255)

                                          tatCAAACGTTATGGCTGAAAGCGCACTAGCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                   >  2:664631/1‑91 (MQ=255)

                                           atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                         >  1:1298521/1‑100 (MQ=255)

                                           atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCAATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                         >  1:85889/1‑100 (MQ=255)

                                            tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATAGTGACAACACCGAGAAGGTAGTGGGTGTTGTGGAt                                                        >  1:1313975/1‑100 (MQ=255)

                                               aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                     <  1:2138835/100‑1 (MQ=255)

                                                  gTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTg                                                   >  2:2247721/1‑99 (MQ=255)

                                                    cgatcttGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTg                                                <  2:1732717/94‑1 (MQ=255)

                                                     aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                               >  1:1338246/1‑100 (MQ=255)

                                                       ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGTGTCTGGTGGATGAGTGCTGGCg                                             <  2:2967187/100‑1 (MQ=255)

                                                       ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGTTGAGTGCTAGCg                                             >  2:290374/1‑100 (MQ=255)

                                                       ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                             <  2:1329177/100‑1 (MQ=255)

                                                       ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGATGGTGGGTATGGTGGATGAGTGCTGGCg                                             >  1:674205/1‑100 (MQ=255)

                                                          tGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgagatcg                                          >  1:1732737/1‑94 (MQ=255)

                                                           gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATg                                         >  1:401884/1‑100 (MQ=255)

                                                            aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgg                                        >  1:454727/1‑100 (MQ=255)

                                                             aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                       <  1:2496760/100‑1 (MQ=255)

                                                             aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                       <  1:1661202/100‑1 (MQ=255)

                                                               gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                     <  2:3008100/100‑1 (MQ=255)

                                                                ctcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCt                                    <  2:53382/98‑1 (MQ=255)

                                                                 gcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg                                   >  2:1354612/1‑100 (MQ=255)

                                                                 gcACTATCCTGCGGAGATTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGCTGAGAGCTGGCGGACGGGGCTg                                   >  2:1483154/1‑100 (MQ=255)

                                                                  cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAAAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCtag                                  >  1:53382/1‑98 (MQ=255)

                                                                      aTCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                              <  1:2217008/100‑1 (MQ=255)

                                                                       tCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATc                             <  1:1098289/100‑1 (MQ=255)

                                                                            cGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcc                        <  1:2035900/100‑1 (MQ=255)

                                                                               aGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacc                     >  1:1656456/1‑100 (MQ=255)

                                                                                gTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGGTGAAAATTCTGcca                       >  2:3039016/1‑97 (MQ=255)

                                                                                gTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                    >  1:2679084/1‑100 (MQ=255)

                                                                                 tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                   <  1:271564/100‑1 (MQ=255)

                                                                                 tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                   <  2:401881/100‑1 (MQ=255)

                                                                                  ttATTGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa                  <  2:636159/100‑1 (MQ=255)

                                                                                  ttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa                  >  1:173649/1‑100 (MQ=255)

                                                                                  ccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGTTGGCGGATGGGGCTGAAAATCCTGCCACCAGa                  >  1:2027730/7‑100 (MQ=255)

                                                                                   tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat                 <  2:1635002/100‑1 (MQ=255)

                                                                                   tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat                 >  1:1676656/1‑100 (MQ=255)

                                                                                   catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat                 >  1:3048510/6‑100 (MQ=255)

                                                                                   catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat                 >  2:2703353/6‑100 (MQ=255)

                                                                                    atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTTGCGGATGGGGCTGAAAATCCTGCCACCAGata                >  1:2153350/5‑100 (MQ=255)

                                                                                    atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGata                >  1:1946293/5‑100 (MQ=255)

                                                                                     tatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTAGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGatat               >  1:1513139/4‑100 (MQ=255)

                                                                                            ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACt          >  1:2739310/1‑98 (MQ=255)

                                                                                            ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc        <  1:995033/100‑1 (MQ=255)

                                                                                            ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTACTAGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc        >  1:767713/1‑100 (MQ=255)

                                                                                              tAATGACTGCTGATATGGATAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGTCATCAGatat               <  2:2607747/91‑1 (MQ=255)

                                                                                              tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg      <  2:2749885/100‑1 (MQ=255)

                                                                                              tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg      <  2:66568/100‑1 (MQ=255)

                                                                                              tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTAGATGAGTGCTAACGAATGGGACTGAAAATCCTGCCACCAGATATAAACTcc        >  1:851682/1‑98 (MQ=255)

                                                                                               aaTGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGt     <  1:664635/100‑1 (MQ=255)

                                                                                                 tGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAATGTTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCt   >  2:2434186/1‑100 (MQ=255)

                                                                                                  gACCGCCGATATGGACAACACCGAGAAGGTGGTAGGTCTGGTGGATGAGTGCTAGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCtt  >  2:45743/1‑100 (MQ=255)

GGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTT  >  NC_000913/207388‑207585

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