Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,517	T→G	6.9%	I292S (ATC→AGC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,517	0	T	G	6.9%	242.2 / 9.9	88	I292S (ATC→AGC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base G (2/4);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GGCCGCCCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGAC  >  NC_002947/1550372‑1550659                                                                                                                                                  |                                                                                                                                               ggctgCCGACTGCTGATCGTTGCCGAAGACGTTGAAGGTGGAGCGCTGGCTCCCATGGTAGTCAAAAACCTCAGATGCATGGTCAAGCTTGCTGCAGTTAAGTGTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCt                                                                                                                                            >  2:1174017/1‑150 (MQ=255)   ccgccCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                          >  1:50170/1‑150 (MQ=255)   ccgccCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                          <  1:1443250/150‑1 (MQ=255)       ccACTGCTGATCGTTGACGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTg                                                                                                                                      >  2:2021133/1‑150 (MQ=255)        cACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGGTGGCTACCCTGGTAGTCAACCCCATCTGCAGCTGCCTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCTgc                                                                                                                                     <  1:2223122/150‑2 (MQ=255)         aCTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCTgca                                                                                                                                    <  2:2250382/150‑3 (MQ=255)             ctgATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCACGGCCATGCTGCAGGACATCGCTGTCCTGACTggc                                                                                                                                <  2:848014/150‑1 (MQ=255)               gATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAATGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTggcgg                                                                                                                              >  1:1375205/1‑150 (MQ=255)                   gTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAg                                                                                                                          <  1:1310375/150‑1 (MQ=255)                      gCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTc                                                                                                                       <  2:2050310/150‑1 (MQ=255)                           aGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCAtctc                                                                                                                  <  2:1376282/150‑1 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTc                                                                                                                       >  2:1626472/1‑141 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTc                                                                                                                       <  1:1626472/141‑1 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                              <  2:1375205/150‑1 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                              >  1:2385834/1‑150 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                              >  1:1985775/1‑150 (MQ=255)                                  gAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaagaa                                                                                                           >  1:787603/1‑150 (MQ=255)                                  gAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaagaa                                                                                                           >  2:330322/1‑150 (MQ=255)                                   aaGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGaaa                                                                                                          >  2:1294353/1‑150 (MQ=255)                                      gTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCg                                                                                                       >  1:538949/1‑150 (MQ=255)                                      gTGAAGCGCTGGCCACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTGCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCg                                                                                                       <  2:804859/150‑1 (MQ=255)                                        gAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGc                                                                                                     <  2:50170/150‑1 (MQ=255)                                          aGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCt                                                                                                   <  2:1985775/150‑1 (MQ=255)                                                 gCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTg                                                                                            <  1:614090/150‑1 (MQ=255)                                                    aCCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCt                                                                                                                                            >  1:1585048/1‑99 (MQ=255)                                                    aCCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCt                                                                                                                                            <  2:1585048/99‑1 (MQ=255)                                                     cccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAg                                                                                                                          >  2:168914/1‑116 (MQ=255)                                                     cccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAg                                                                                                                          <  1:168914/116‑1 (MQ=255)                                                      ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAAc                                                                                       >  2:2265262/1‑150 (MQ=255)                                                      ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAAc                                                                                       <  1:330322/150‑1 (MQ=255)                                                      ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAACGCCATGCTGCAGGGCAACGCTTTCCTGACTGGCGGCCAGATCATCTCCGAATAAATCGGCCTGTCCCTGGacac                                                                                       >  2:2309126/1‑150 (MQ=255)                                                       cTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAAcc                                                                                      >  1:796084/1‑150 (MQ=255)                                                        tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                     >  2:224896/1‑150 (MQ=255)                                                        tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                     >  1:2010437/1‑150 (MQ=255)                                                        tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCCCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                     >  2:1273747/1‑150 (MQ=255)                                                           tAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTa                                                                                  <  2:787603/150‑1 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCg                                                                                                                                              <  1:632756/85‑1 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCg                                                                                                                                              >  2:632756/1‑85 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCt                                                                                                   <  2:998461/128‑1 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCt                                                                                                   >  1:998461/1‑128 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTcc                                                                                               >  1:792448/1‑132 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTcc                                                                                               <  2:792448/132‑1 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACTGCTACCCTg                                                                             >  2:2369310/1‑150 (MQ=255)                                                                          gcgGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTgg                                                                   <  1:651154/150‑1 (MQ=255)                                                                          gcgGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCt                                                                                                                                       >  2:1896345/1‑82 (MQ=255)                                                                          gcgGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCt                                                                                                                                       <  1:1896345/82‑1 (MQ=255)                                                                            ggCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGt                                                                 <  1:224896/150‑1 (MQ=255)                                                                            ggCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGt                                                                 <  2:796084/150‑1 (MQ=255)                                                                               aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                              >  2:1187008/1‑150 (MQ=255)                                                                               aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                              <  1:1174918/150‑1 (MQ=255)                                                                               aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                              <  2:234432/150‑1 (MQ=255)                                                                                   tCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCACAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATATCCGAAGAAATCGGCCTGTCCCTGGCAACCGCTACCCTGGCGCACCTGGGTAACGCCa                                                          >  2:124566/1‑150 (MQ=255)                                                                                             tgcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCAt                                                >  2:1659870/1‑150 (MQ=255)                                                                                              gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                               <  2:1652375/150‑1 (MQ=255)                                                                                              gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                               <  1:1294353/150‑1 (MQ=255)                                                                                                 gTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGCCCAGCTCATCTCCGAAGAACTCGGCCTGTCCCTGGAATCCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTg                                            <  2:579129/150‑1 (MQ=255)                                                                                                  tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGt                                           <  2:2010437/150‑1 (MQ=255)                                                                                                  tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCGCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGt                                           >  1:543404/1‑150 (MQ=255)                                                                                                    aaGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTGCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTcc                                         >  2:425762/1‑150 (MQ=255)                                                                                                    aaGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGc                                                                                                     >  2:1223049/1‑90 (MQ=255)                                                                                                    aaGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGc                                                                                                     <  1:1223049/90‑1 (MQ=255)                                                                                                     aGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCATCTGGGTAACGCCAAGCGCGTCATCCTGTCCa                                        <  2:2012891/150‑1 (MQ=255)                                                                                                      ggCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCaa                                       <  2:133323/150‑1 (MQ=255)                                                                                                       gCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAg                                      >  1:2240099/1‑150 (MQ=255)                                                                                                         tCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                     <  1:757681/101‑1 (MQ=255)                                                                                                         tCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                     >  2:757681/1‑101 (MQ=255)                                                                                                           cGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaa                                  >  2:728331/1‑150 (MQ=255)                                                                                                            gggCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaaa                                 <  1:1273747/150‑1 (MQ=255)                                                                                                               cTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                              <  2:1087096/150‑1 (MQ=255)                                                                                                                  cGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                           <  2:538949/150‑1 (MQ=255)                                                                                                                  cGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                           <  1:2021133/150‑1 (MQ=255)                                                                                                                  cGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTCTCTTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                           <  1:2309126/150‑1 (MQ=255)                                                                                                                    gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                           <  1:1845256/148‑1 (MQ=255)                                                                                                                    gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                           >  2:1845256/1‑148 (MQ=255)                                                                                                                         cgccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCatcatc                    <  2:1453748/150‑1 (MQ=255)                                                                                                                          gccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTggg                                                                  >  1:1407771/1‑103 (MQ=255)                                                                                                                          gccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTggg                                                                  <  2:1407771/103‑1 (MQ=255)                                                                                                                               aaGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACaccacc                          >  1:203816/1‑138 (MQ=255)                                                                                                                               aaGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACaccacc                          <  2:203816/138‑1 (MQ=255)                                                                                                                               aaGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGt              <  1:1187008/150‑1 (MQ=255)                                                                                                                                aGGCCATGCTGCAGGACATCGCTGTCCTCACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTg             >  1:430532/1‑150 (MQ=255)                                                                                                                                 ggCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGc            >  2:456077/1‑150 (MQ=255)                                                                                                                                   ccATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTg          >  1:138707/1‑150 (MQ=255)                                                                                                                                        ctgcAGTACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt     <  1:839164/150‑1 (MQ=255)                                                                                                                                         tgcAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTg    >  1:895405/1‑150 (MQ=255)                                                                                                                                          tcAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                              <  2:2361755/122‑1 (MQ=255)                                                                                                                                          tcAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                              >  1:2361755/2‑123 (MQ=255)                                                                                                                                           cAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGac  >  2:2119270/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                               GGCCGCCCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGAC  >  NC_002947/1550372‑1550659

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