BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,695	T→G	15.2%	S123A (TCC→GCC) ‡	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,695	0	T	G	15.2%	452.3 / 63.8	164	S123A (TCC→GCC) ‡	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base T (60/79); new base G (8/17); total (68/96)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATC  >  NC_002947/549549‑549844

   gTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTc                                                                                                                                                 >  1:315105/1‑150 (MQ=34)

     cACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAg                                                                                                                                               <  1:3567055/150‑1 (MQ=34)

      aCTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAgg                                                                                                                                              >  2:1206205/1‑150 (MQ=34)

       cTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGt                                                                                                                                             >  2:2672153/1‑150 (MQ=34)

        tACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTa                                                                                                                                            <  1:2369154/150‑1 (MQ=34)

             tCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGt                                                                                                                                       >  1:4962985/1‑150 (MQ=34)

                 gTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                                                                                                   <  2:2906263/150‑1 (MQ=34)

                     aCTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACa                                                                                                                               >  2:1363742/1‑150 (MQ=34)

                       tGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATc                                                                                                                             >  1:3312396/1‑150 (MQ=34)

                        gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                                            >  1:1188918/1‑150 (MQ=34)

                        gCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCg                                                                                                                            <  1:2539827/150‑1 (MQ=34)

                         cccAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGt                                                                                                                           >  2:2562115/1‑150 (MQ=34)

                          ccAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCTCGTCAGGTAGGCGTTCCGTACATCGTg                                                                                                                          <  2:574357/150‑1 (MQ=21)

                            aGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGt                                                                                                                        <  2:3990747/150‑1 (MQ=34)

                                cACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTc                                                                                                                    <  2:2827187/150‑1 (MQ=34)

                                 aCGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTcc                                                                                                                   <  1:2268983/150‑1 (MQ=34)

                                  cGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCt                                                                                                                  <  1:327763/150‑1 (MQ=34)

                                   gCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTg                                                                                                                 >  2:3960075/1‑150 (MQ=34)

                                          aCGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAgg                                                                                                          <  2:5101161/150‑1 (MQ=34)

                                           cGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGc                                                                                                         >  1:3380190/1‑150 (MQ=34)

                                            gTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGcc                                                                                                        >  2:5194281/1‑150 (MQ=35)

                                              gaagaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGa                                                                                                      >  2:2084506/1‑150 (MQ=37)

                                                    cATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                                                <  2:1448833/150‑1 (MQ=37)

                                                     aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                                                                                               <  2:3050406/150‑1 (MQ=37)

                                                         tCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgacg                                                                                           >  1:4907852/1‑150 (MQ=37)

                                                                tgctgcCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGa                                                                                    <  1:3279571/150‑1 (MQ=37)

                                                                 gctgcCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAg                                                                                   >  1:2685745/1‑150 (MQ=37)

                                                                   tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact                                                                                 >  1:570552/1‑150 (MQ=2)

                                                                     cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                                               <  1:3924103/150‑1 (MQ=37)

                                                                     cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                                               >  2:4097632/1‑150 (MQ=37)

                                                                     cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                                               <  1:995926/150‑1 (MQ=37)

                                                                     cacAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttc                                                                               >  2:480945/3‑148 (MQ=2)

                                                                      ccAGATGGACGGCGCTATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                                              >  1:430204/1‑150 (MQ=37)

                                                                      ccAGATGGACGGCGCTATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                                              >  1:3158865/1‑150 (MQ=37)

                                                                      ccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                                              >  1:747200/1‑150 (MQ=37)

                                                                       cAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttctc                                                                             <  1:2741506/150‑2 (MQ=2)

                                                                       cAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttctc                                                                             <  1:2741159/150‑2 (MQ=2)

                                                                        aGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTgg                                                                            <  1:2561002/150‑1 (MQ=37)

                                                                          aTGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGaa                                                                          <  2:1133813/150‑1 (MQ=38)

                                                                          aTGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGaa                                                                          <  1:2237055/150‑1 (MQ=38)

                                                                           tGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAAc                                                                         <  1:2766772/150‑1 (MQ=38)

                                                                           tGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGagc                                                                         >  1:850276/1‑148 (MQ=2)

                                                                               cGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGt                                                                     >  1:5249573/1‑150 (MQ=39)

                                                                               cGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGt                                                                     >  2:3686878/1‑150 (MQ=39)

                                                                               cGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGt                                                                     <  1:2930077/150‑1 (MQ=255)

                                                                                ggCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGCCGCTGAGCTGCTGGAACTGGTc                                                                    <  1:2335151/150‑1 (MQ=255)

                                                                                  cgcgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCga                                                                  <  1:1348759/150‑1 (MQ=39)

                                                                                   gcgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgag                                                                 >  2:4004716/1‑150 (MQ=39)

                                                                                   gcgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgag                                                                 <  2:2298385/150‑1 (MQ=39)

                                                                                    cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                <  2:4579730/150‑1 (MQ=39)

                                                                                    cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                <  2:5103261/150‑1 (MQ=39)

                                                                                     gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                                               <  2:2230723/150‑1 (MQ=39)

                                                                                     gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                                               <  2:1052779/150‑1 (MQ=39)

                                                                                       tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                                             <  2:1998407/150‑1 (MQ=39)

                                                                                       tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                                             <  2:3312396/150‑1 (MQ=39)

                                                                                       tCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAgaactcg                                                             <  1:5369862/150‑6 (MQ=255)

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

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

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

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

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

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

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

                                                                                                 cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                                                                                                   >  1:5212988/1‑70 (MQ=255)

                                                                                                 cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCg                                                                                                                                   <  2:5212988/70‑1 (MQ=255)

                                                                                                 cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                                                   >  2:1282437/1‑150 (MQ=39)

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

                                                                                                   cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                                                 >  1:1626241/1‑150 (MQ=39)

                                                                                                   cGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGaac                                                 <  1:5060636/150‑3 (MQ=255)

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

                                                                                                    ggCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTgacgatg                                                                                        <  2:2965963/110‑3 (MQ=255)

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

                                                                                                        gccgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgctg                                            <  1:1638368/150‑1 (MQ=39)

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

                                                                                                           gATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTcc                                         <  1:2452379/150‑1 (MQ=39)

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

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

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

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

                                                                                                               gTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCt                                     >  1:322605/1‑150 (MQ=39)

                                                                                                                tCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTa                                    <  1:2300129/150‑1 (MQ=39)

                                                                                                                  cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                                                >  1:4953369/1‑88 (MQ=255)

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

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

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

                                                                                                                         acagaCCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTccc                           >  2:2124084/5‑150 (MQ=255)

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

                                                                                                                          cAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCa                          <  2:3218112/150‑1 (MQ=37)

                                                                                                                               ccGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGcga                     >  1:585500/1‑150 (MQ=37)

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

                                                                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgct                                             <  1:78147/120‑1 (MQ=39)

                                                                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgct                                             >  2:78147/1‑120 (MQ=39)

                                                                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGacac               >  1:1100095/1‑150 (MQ=37)

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

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

                                                                                                                                          tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa          <  2:4314122/150‑1 (MQ=37)

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

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

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

                                                                                                                                                  tCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcatc  >  2:501842/1‑150 (MQ=34)

TTCGTCACTACGCTCACGTTGACTGCCCAGGTCACGCTGACTACGTGAAGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATC  >  NC_002947/549549‑549844

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