Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	666,572	C→T	8.5%	S1146S (AGC→AGT)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	666,572	0	C	T	8.5%	164.4 / 3.2	59	S1146S (AGC→AGT)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/3);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AGTGGACGCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCAAGGGCGAGCCGGCG  >  NC_002947/666427‑666719                                                                                                                                                  |                                                                                                                                                    aGTGGACGCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAgccg                                                                                                                                                   <  7:147089/148‑1 (MQ=255)   tGGACGCTGGGCGTGCGCCCGGCCTACCCTGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAgccggc                                                                                                                                                 <  5:127074/148‑1 (MQ=255)     gACGCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAgccggcc                                                                                                                                                <  3:94688/147‑1 (MQ=255)     gACGCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAgccggcc                                                                                                                                                >  4:94688/1‑147 (MQ=255)        gCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCcgccgc                                                                                                                                           >  7:52908/1‑149 (MQ=255)                     ccGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGa                                                                                                                              <  4:239328/149‑1 (MQ=255)                          cTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCaaa                                                                                                                         >  1:147539/1‑149 (MQ=255)                               cGGCAATGCTCAAGCGCTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTcgc                                                                                                                    >  3:15601/1‑149 (MQ=255)                                ggCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTcgcg                                                                                                                   >  5:145561/1‑149 (MQ=255)                                 gCAATGCTCAAGCACTACCGTGTGTCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTcgcgc                                                                                                                  <  2:82386/149‑1 (MQ=255)                                   aaTGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCgcgcgc                                                                                                                >  4:2616/1‑149 (MQ=255)                                        tCAAGCACTACCGTGTGGCCTTGAAGGCCCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGCGACGGCCGCCGCTGAACCTTGCCGAGCACATTCGCCCGCTGCaa                                                                                                           >  7:150171/1‑149 (MQ=255)                                               cTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTAcc                                                                                                    >  5:99618/1‑149 (MQ=255)                                               cTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTAcc                                                                                                    <  7:296344/149‑1 (MQ=255)                                                tACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTAccc                                                                                                   >  6:304611/1‑149 (MQ=255)                                                  ccGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTAcccc                                                                                                  >  5:117776/1‑148 (MQ=255)                                                    gtgtGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGACCCTTGCCGAGCAAATTCGCCCGCAGCACCCCTACCCCTAt                                                                                               >  4:224471/1‑149 (MQ=255)                                                       tgGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGc                                                                                            <  6:69214/149‑1 (MQ=255)                                                         gCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCt                                                                                           <  5:93991/148‑1 (MQ=255)                                                         gCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTg                                                                                          >  3:39373/1‑149 (MQ=255)                                                         gCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTg                                                                                          >  4:17519/1‑149 (MQ=255)                                                          ccTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGc                                                                                         >  3:23106/1‑149 (MQ=255)                                                          ccTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGc                                                                                         <  4:24174/149‑1 (MQ=255)                                                              gAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCt                                                                                       <  7:265892/147‑1 (MQ=255)                                                                aGGACCAGCCGTGGTCTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGTCCTGGAATCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGATGGTTGACGAGGAATTTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGaa                                                                                   <  4:142466/149‑1 (MQ=255)                                                                aGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGaa                                                                                   >  4:183507/1‑149 (MQ=255)                                                                 ggACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAAc                                                                                  <  1:138341/149‑1 (MQ=255)                                                                      aGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTCGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGac                                                                              >  6:56299/1‑148 (MQ=255)                                                                          gTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGTCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGaccac                                                                           >  2:148437/1‑147 (MQ=255)                                                                          gTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGTCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGaccac                                                                           >  2:142785/1‑147 (MQ=255)                                                                          gTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGaccac                                                                           >  4:354257/1‑147 (MQ=255)                                                                               ctctGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGTCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGc                                                                    <  4:10957/149‑1 (MQ=255)                                                                                  tGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTg                                                                 >  5:177850/1‑149 (MQ=255)                                                                                    cccGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCt                                                                  <  3:17519/146‑1 (MQ=255)                                                                                    cccGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTa                                                               <  8:52908/149‑1 (MQ=255)                                                                                    cccGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTa                                                               <  6:117776/149‑1 (MQ=255)                                                                                     ccGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGGACAGACCACCAGTGGCCTGTAc                                                              >  7:356166/1‑149 (MQ=255)                                                                                     ccGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTAc                                                              >  1:251025/1‑149 (MQ=255)                                                                                      cGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTAcc                                                             >  6:35532/1‑149 (MQ=255)                                                                                        aaGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCt                                                           <  4:39373/149‑1 (MQ=255)                                                                                         aGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCtt                                                          <  6:177850/149‑1 (MQ=255)                                                                                          gCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTc                                                         <  6:145561/149‑1 (MQ=255)                                                                                            aGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGc                                                       <  5:304611/149‑1 (MQ=255)                                                                                             gACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCt                                                      <  4:194390/149‑1 (MQ=255)                                                                                              aCCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTg                                                     >  2:145632/1‑149 (MQ=255)                                                                                                   gccgccTTCGAGCCCGCTGGCCTGGAAGGCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACg                                                 <  4:382179/148‑1 (MQ=255)                                                                                                     cgccTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGcc                                               >  7:26265/1‑148 (MQ=255)                                                                                                           cGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGTCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGc                                          <  1:142785/147‑1 (MQ=255)                                                                                                           cGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGTCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGc                                          <  1:148437/147‑1 (MQ=255)                                                                                                           cGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGc                                          <  3:354257/147‑1 (MQ=255)                                                                                                              gCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACa                                     <  3:2616/149‑1 (MQ=255)                                                                                                                gcGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAg                                    <  4:15601/147‑1 (MQ=255)                                                                                                                ccGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGGACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGc                                   <  8:356166/149‑1 (MQ=255)                                                                                                                      gCCTGGAAGGCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCaa                              <  3:224471/148‑1 (MQ=255)                                                                                                                      gCCTGGAAGCCATCCTGGCCTTGTCCAGCCGGCCGCCGCAGATGCTTGCCGAGCTAATTGGCGCGCGTCAAGCGTACCCCTATCGCTGCCTCCATCAGACCACCAGTGGCCTGTACCCTTCGCTGTTCGCCGATGCCGACAGCCTCaa                              <  8:150171/148‑1 (MQ=255)                                                                                                                      gCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCCAACAGACCTCCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAg                             <  2:147539/149‑1 (MQ=255)                                                                                                                               ccagccTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTggg                     >  2:78147/1‑148 (MQ=255)                                                                                                                                  gccTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATc                 <  8:26265/149‑1 (MQ=255)                                                                                                                                   ccTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCAGTCCCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATc                 >  7:161436/1‑148 (MQ=255)                                                                                                                                    cgggCCTGGTCGAGCCGGCCGCCGCGGAACCTTTCCGCGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCaa               <  6:99618/147‑1 (MQ=255)                                                                                                                                        ccTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCAAGGGc           >  7:135426/1‑149 (MQ=255)                                                                                                                                          ttGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGCCCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCACGCGCCTGGGCATCAAGGGc           >  4:223231/1‑147 (MQ=255)                                                                                                                                               gAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCAAGGGCGAGCcgg    <  5:35532/149‑1 (MQ=255)                                                                                                                                                 gccggccgCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCAAGGGCGAGCcggcg  <  3:223231/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGTGGACGCTGGGCGTGCGCCCGGCCTACCCGGCAATGCTCAAGCACTACCGTGTGGCCTTGAAGGACCAGCCGTGGACTCTGCCCGAAGCAGACCTGGCCGCCTTCGAGCCCGCTGGCCTGGAAGCCAGCCTGGCCTTGTCGAGCCGGCCGCCGCTGAACCTTGCCGAGCAAATTCGCGCGCTGCAAGCCTACCCCTATGGCTGCCTCGAACAGACCACCAGTGGCCTGTACCCTTCGCTGTACGCCGATGCCGACAGCCTCAAGCGCCTGGGCATCAAGGGCGAGCCGGCG  >  NC_002947/666427‑666719

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