Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,520	C→A	6.6%	A293D (GCT→GAT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,520	0	C	A	6.6%	280.6 / 6.1	90	A293D (GCT→GAT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base A (2/4);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GGCCGCCCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATC  >  NC_002947/1550372‑1550668                                                                                                                                                     |                                                                                                                                                     ggctgCCGACTGCTGATCGTTGCCGAAGACGTTGAAGGTGGAGCGCTGGCTCCCATGGTAGTCAAAAACCTCAGATGCATGGTCAAGCTTGCTGCAGTTAAGTGTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCt                                                                                                                                                     >  2:1174017/1‑150 (MQ=255)   ccgccCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                                   <  1:1443250/150‑1 (MQ=255)   ccgccCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                                   >  1:50170/1‑150 (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                                                                                                                                <  1:1626472/141‑1 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTc                                                                                                                                >  2:1626472/1‑141 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                       <  2:1375205/150‑1 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                       >  1:2385834/1‑150 (MQ=255)                               gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                       >  1:1985775/1‑150 (MQ=255)                                  gAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaagaa                                                                                                                    >  2:330322/1‑150 (MQ=255)                                  gAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaagaa                                                                                                                    >  1:787603/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                                                                                                                                                     <  2:1585048/99‑1 (MQ=255)                                                    aCCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCt                                                                                                                                                     >  1:1585048/1‑99 (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                                                                                              >  1:2010437/1‑150 (MQ=255)                                                        tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                              >  2:224896/1‑150 (MQ=255)                                                        tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCCCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                              >  2:1273747/1‑150 (MQ=255)                                                           tAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTa                                                                                           <  2:787603/150‑1 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCt                                                                                                            >  1:998461/1‑128 (MQ=255)                                                                aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCt                                                                                                            <  2:998461/128‑1 (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:234432/150‑1 (MQ=255)                                                                               aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                                       <  1:1174918/150‑1 (MQ=255)                                                                               aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                                       >  2:1187008/1‑150 (MQ=255)                                                                                   tCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCACAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATATCCGAAGAAATCGGCCTGTCCCTGGCAACCGCTACCCTGGCGCACCTGGGTAACGCCa                                                                   >  2:124566/1‑150 (MQ=255)                                                                                             tgcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCAt                                                         >  2:1659870/1‑150 (MQ=255)                                                                                              gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                                        <  1:1294353/150‑1 (MQ=255)                                                                                              gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                                        <  2:1652375/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                                                                                              >  2:757681/1‑101 (MQ=255)                                                                                                         tCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                              <  1:757681/101‑1 (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                                    >  2:1845256/1‑148 (MQ=255)                                                                                                                    gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcaccac                                    <  1:1845256/148‑1 (MQ=255)                                                                                                                         cgccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCatcatc                             <  2:1453748/150‑1 (MQ=255)                                                                                                                          gccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTggg                                                                           <  2:1407771/103‑1 (MQ=255)                                                                                                                          gccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTggg                                                                           >  1:1407771/1‑103 (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                                       >  1:2361755/2‑123 (MQ=255)                                                                                                                                          tcAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                                       <  2:2361755/122‑1 (MQ=255)                                                                                                                                           cAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGac           >  2:2119270/1‑150 (MQ=255)                                                                                                                                                   cGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGAt   >  1:1414008/1‑150 (MQ=255)                                                                                                                                                    gCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACg                         <  1:1168379/127‑1 (MQ=255)                                                                                                                                                    gCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATc  <  2:1765698/150‑1 (MQ=255)                                                                                                                                                    gCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAAACCACCATCATCGACg                         >  2:1168379/1‑127 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGCCGCCCACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATC  >  NC_002947/1550372‑1550668

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