Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	418,295	C→G	8.1%	L26L (CTC→CTG)	PP_0345 →	heptose kinase WapP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	418,295	0	C	G	8.1%	262.0 / 10.1	84	L26L (CTC→CTG)	PP_0345	heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base G (4/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

CTGCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCA  >  NC_002947/418147‑418441                                                                                                                                                     |                                                                                                                                                    ctgctCGCCTGGCTGGAGTGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTc                                                                                                                                                     <  1:216437/149‑1 (MQ=255)    ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGTGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAg                                                                                                                                                  >  1:365860/1‑149 (MQ=255)    ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAg                                                                                                                                                  >  7:57273/1‑149 (MQ=255)      cGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAggg                                                                                                                                                <  2:419751/149‑1 (MQ=255)       gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                 <  1:140189/147‑1 (MQ=255)       gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGc                                                                                                                                               >  1:226663/1‑149 (MQ=255)       gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGATCGCCCTCCAgg                                                                                                                                                 <  3:416748/147‑1 (MQ=255)         ctggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCg                                                                                                                                              >  6:378826/1‑148 (MQ=255)            gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGaa                                                                                                                                            <  6:24743/147‑1 (MQ=255)            gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACg                                                                                                                                          >  5:248873/1‑149 (MQ=255)            gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACg                                                                                                                                          >  1:218574/1‑149 (MQ=255)            gcgggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAAc                                                                                                                                           <  6:335425/145‑1 (MQ=255)             ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                         >  8:70312/1‑149 (MQ=255)              tggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                        >  4:267587/1‑149 (MQ=255)              tggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGAGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                        >  8:404694/1‑149 (MQ=255)                gAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAGCGGCTg                                                                                                                                      >  4:305485/1‑149 (MQ=255)                gAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGCGTTCGCCCTGGAGGGCGAAcaccgc                                                                                                                                      <  3:414972/149‑6 (MQ=255)                 aGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTga                                                                                                                                     <  4:401331/149‑1 (MQ=255)                 aGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTga                                                                                                                                     <  4:394973/149‑1 (MQ=255)                     cAAGGCCCAGAAACTCTACGATAGAAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGACGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccg                                                                                                                                 >  6:210939/1‑149 (MQ=255)                          cccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacc                                                                                                                            >  8:304374/1‑149 (MQ=255)                          cccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacc                                                                                                                            <  2:136916/149‑1 (MQ=255)                           ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgaccc                                                                                                                           >  2:308700/1‑149 (MQ=255)                           ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgaccc                                                                                                                           >  8:351119/1‑149 (MQ=255)                             aGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccgc                                                                                                                         >  5:25900/1‑149 (MQ=255)                              gAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCt                                                                                                                        >  6:343926/1‑149 (MQ=255)                                aaCTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCACGGCAAACGGCTGACCCGCGACCCGCTCa                                                                                                                      <  1:245797/149‑1 (MQ=255)                                   tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGt                                                                                                                    >  6:64037/1‑148 (MQ=255)                                   tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGt                                                                                                                    >  8:303235/1‑148 (MQ=255)                                     tccGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACGTCGGCAGCCTCGACGCGGTGTTCGCGCTCCAGGGCGAACGGCTGACCCGAGACGCGCGCAGTg                                                                                                                   >  6:334866/3‑147 (MQ=255)                                     tACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTg                                                                                                                   >  5:348340/1‑147 (MQ=255)                                       cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                 <  6:248873/147‑1 (MQ=255)                                       cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAg                                                                                                                <  7:70312/148‑1 (MQ=255)                                        gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                              <  4:227760/149‑1 (MQ=255)                                        gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                              >  7:268720/1‑149 (MQ=255)                                            gcaagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTa                                                                                                           <  8:268720/148‑1 (MQ=255)                                             caagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACg                                                                                                         >  8:82205/1‑149 (MQ=255)                                              aagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccg                                                                                                                                 <  2:5065/124‑1 (MQ=255)                                              aagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccg                                                                                                                                 >  1:5065/1‑124 (MQ=255)                                                       tATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGt                                                                                               >  5:29757/1‑149 (MQ=255)                                                         tGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTg                                                                                              <  7:351119/148‑1 (MQ=255)                                                           gggATGCACTCTGATGGCGAGTTGGACACTTGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                           >  8:81335/1‑149 (MQ=255)                                                           gggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                           <  5:343926/149‑1 (MQ=255)                                                               tGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGAt                                                                                       <  6:29757/149‑1 (MQ=255)                                                                   ctctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                   <  7:134040/149‑1 (MQ=255)                                                                   ctctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                   <  8:57273/149‑1 (MQ=255)                                                                           gCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGt                                                                           >  8:127382/1‑149 (MQ=255)                                                                                 tGGACACTGGCGCCGGGGTACTAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGTATCAACTACTACGTCAAg‑c                                                                     <  2:218574/148‑1 (MQ=255)                                                                                   gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCt                                                                  <  7:304374/149‑1 (MQ=255)                                                                                    acacTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTa                                                                 >  3:54246/1‑149 (MQ=255)                                                                                     cacTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTac                                                                <  5:426322/149‑1 (MQ=255)                                                                                           cgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCggg                                                          <  3:96257/149‑1 (MQ=255)                                                                                           cgcAGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGTCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCggg                                                          <  5:210939/149‑1 (MQ=255)                                                                                            gcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCgggg                                                         <  5:64037/149‑1 (MQ=255)                                                                                               ggggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAAGGCGAACGGCTGACCCGCGACCCGCTCAGTTAAGTGGTACGCATCTAGCGTGACGGGATAAACTACTACGTCAAG‑CCCTACACCGGCGCCg                                                      >  8:152102/1‑149 (MQ=255)                                                                                                gggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCgg                                                     <  6:348340/149‑1 (MQ=255)                                                                                                      gAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGca                                               >  6:232251/1‑149 (MQ=255)                                                                                                      gAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGca                                               >  3:212960/1‑149 (MQ=255)                                                                                                       aaCACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGcac                                              <  4:54246/149‑1 (MQ=255)                                                                                                         cacCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACAt                                            <  4:212960/149‑1 (MQ=255)                                                                                                         cacCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGACAAGCACAt                                            <  2:226663/149‑1 (MQ=255)                                                                                                                cggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAATGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCgg                                     <  4:75913/149‑1 (MQ=255)                                                                                                                        ttCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTggg                             >  1:57046/1‑149 (MQ=255)                                                                                                                            gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                         <  5:378826/149‑1 (MQ=255)                                                                                                                            gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                         >  2:163551/1‑149 (MQ=255)                                                                                                                             cAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                        <  7:152102/149‑1 (MQ=255)                                                                                                                              aGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                        <  1:308700/149‑1 (MQ=255)                                                                                                                                ccTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCAc                     <  7:82205/149‑1 (MQ=255)                                                                                                                                 cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                    >  8:46008/1‑149 (MQ=255)                                                                                                                                  tCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                    >  8:10878/1‑148 (MQ=255)                                                                                                                                   cGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCa                  >  7:281595/1‑149 (MQ=255)                                                                                                                                    gACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATAGACCGTGACGGGATCAACTACTACGTCAAC‑CGCTACACCGGGGCCCGCAAGCACAAGCGCCCGTACCTGGGGCGGCCACGCa                  >  6:274143/1‑148 (MQ=255)                                                                                                                                       gcgGGGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCaa              <  5:129235/149‑1 (MQ=255)                                                                                                                                        cgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg             <  2:57046/149‑1 (MQ=255)                                                                                                                                         ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAgg            <  2:365860/149‑1 (MQ=255)                                                                                                                                         ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAgg            >  7:171191/1‑149 (MQ=255)                                                                                                                                         ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGCTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAgg            <  3:305485/149‑1 (MQ=255)                                                                                                                                            gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGcc          >  3:85001/1‑148 (MQ=255)                                                                                                                                            gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGcc          >  4:414875/1‑148 (MQ=255)                                                                                                                                            gtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGcc          >  6:353881/1‑148 (MQ=255)                                                                                                                                             ttCGCCCTACAGGGCGAACGGCTGACCCGCTACCCGCTCAGTGAAGTGGTAAGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACAACGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCGACGCATCAAGGACg         >  8:197446/1‑148 (MQ=255)                                                                                                                                               cGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAAt      >  2:419967/1‑149 (MQ=255)                                                                                                                                                gCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATg     >  4:264147/1‑149 (MQ=255)                                                                                                                                                gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAAt      <  7:127382/148‑1 (MQ=255)                                                                                                                                                   cTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTTACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCa  <  7:404694/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAG‑CGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCA  >  NC_002947/418147‑418441

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