Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	336,125	+T	coding (210/237 nt)	PP_0278 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	336,124	1	.	T	97.1%	114.4 / NA	36	D70V (GAT→GTT)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base T (11/23);  minor base G (0/1);  total (11/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335989‑336263                                                                                                                                          |                                                                                                                                             tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCg                                                                                                                                    >  1:93215/1‑147 (MQ=255)     tctcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATg                                                                                                                               <  3:71079/148‑1 (MQ=255)       tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATgg                                                                                                                              >  1:111074/1‑147 (MQ=255)          tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCaa                                                                                                                          >  2:16801/1‑148 (MQ=255)                   agCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                 <  3:122148/148‑1 (MQ=255)                       cGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAAc                                                                                                             <  6:120668/148‑1 (MQ=255)                           gTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGata                                                                                                         <  5:155991/148‑1 (MQ=255)                              cccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCtg                                                                                                                                        >  4:38075/1‑114 (MQ=255)                               ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGTCCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGTCAAGGCGCTGATGAACGatat                                                                                                        >  1:52387/1‑145 (MQ=255)                               ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCtgt                                                                                                                                       <  3:38075/114‑1 (MQ=255)                                   cGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATc                                                                                                 <  3:56964/148‑1 (MQ=255)                                             tACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAggg                                                                                         >  1:156239/1‑146 (MQ=255)                                                    aTGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCCGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCa                                                                                >  3:14762/1‑148 (MQ=255)                                                        ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGTCCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGTCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                            <  2:52387/148‑1 (MQ=255)                                                        ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                            <  3:129212/148‑1 (MQ=255)                                                          cGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                          >  4:3835/1‑148 (MQ=255)                                                                           aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                          <  2:93215/147‑1 (MQ=255)                                                                           aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                         <  3:159882/148‑1 (MQ=255)                                                                           aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                         <  3:157401/148‑1 (MQ=255)                                                                             cAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGTGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                       <  5:26543/148‑1 (MQ=255)                                                                              aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                      >  7:95336/1‑148 (MQ=255)                                                                               aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTggg                                                     <  7:37423/148‑1 (MQ=255)                                                                                        tGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                            <  3:89830/148‑1 (MQ=255)                                                                                        tGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                            <  5:63152/148‑1 (MQ=255)                                                                                        ggcgcGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                            <  7:140068/147‑1 (MQ=255)                                                                                               gggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑cgc                                    <  8:55325/148‑1 (MQ=255)                                                                                               gggCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑cgc                                    <  5:123878/148‑1 (MQ=255)                                                                                                    atatGGACGCCGACCATTAATTCCTGGGGGCAGCCGCT‑TACCTGTTCGAAATGGCCAAGGCGCTGATGAACGATATCGAGCTCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATg                               <  2:70354/148‑1 (MQ=255)                                                                                                     tatGGAGGACGACCATACATTCCGGGGGGCAGCCG‑AGTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATgg                              <  3:161851/148‑1 (MQ=255)                                                                                                           ggaCGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCa                        <  1:104536/148‑1 (MQ=255)                                                                                                                  cATAAATTCCTCGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGTGTAAGTAGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCTT‑CGCGAATGGGCTGCAAagcagc                 <  4:30582/148‑1 (MQ=255)                                                                                                                       aTTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGGCATGCGTGGGCGCATGGG‑TGCAAATCAGcccccg            >  6:137301/1‑146 (MQ=255)                                                                                                                        ttCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATTTCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCAAAGCAGCCCCGGc           <  4:75892/148‑1 (MQ=255)                                                                                                                        ttCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTTCGG‑AGCGACAGGGGTGGAAAGCACCCCCGGc           >  2:104939/1‑148 (MQ=255)                                                                                                                                ggCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCag   <  3:123085/148‑1 (MQ=255)                                                                                                                                 gCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCaga  <  2:111074/148‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             TCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGT‑CGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335989‑336263

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