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.5%	130.2 / NA	41	D70V (GAT→GTT)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base T (19/20);  minor base G (0/1);  total (20/21)
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

AACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTT  >  NC_002947/335979‑336259                                                                                                                                                   |                                                                                                                                          aaCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tt                                                                                                                                         <  6:193669/148‑2 (MQ=255)         ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGc                                                                                                                                 <  6:416186/148‑1 (MQ=255)         ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCg                                                                                                                                <  4:165006/149‑1 (MQ=255)           tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCAACTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGa                                                                                                                               >  8:301876/1‑148 (MQ=255)                 tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGc                                                                                                                         >  8:28652/1‑148 (MQ=255)                    tACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAg                                                                                                                     >  6:146895/1‑149 (MQ=255)                        cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGgcgc                                                                                                                 <  8:299681/149‑1 (MQ=255)                           agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                              >  4:418980/1‑149 (MQ=255)                           agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                              >  7:77035/1‑149 (MQ=255)                           actacccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCtga                                                                                                              >  8:67634/5‑149 (MQ=255)                             tacccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTTTTCGACATGGCCAAGGCGCtgat‑g                                                                                                           <  7:67634/147‑1 (MQ=255)                               cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑Gaa                                                                                                         <  2:90466/149‑1 (MQ=255)                                ccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAAc                                                                                                        <  6:367210/149‑1 (MQ=255)                                      tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATc                                                                                                  <  2:281301/149‑1 (MQ=255)                                         ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCCCCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑TGACCTGTGCGACATGGCCAAGGCGCGGCG‑GAACGATATc                                                                                                  >  2:328974/1‑146 (MQ=255)                                         ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGACTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCgag                                                                                               <  4:17471/149‑1 (MQ=255)                                           gTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGAt                                                                                             <  3:418980/149‑1 (MQ=255)                                              gACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGc                                                                                          <  2:141763/149‑1 (MQ=255)                                                cgcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGcc                                                                                         >  5:271869/1‑148 (MQ=255)                                                   ggCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAgg                                                                                     >  2:226343/1‑149 (MQ=255)                                                     ccTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGt                                                                                   <  5:443720/149‑1 (MQ=255)                                                      cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tta                                                                                                                                        >  5:227541/1‑94 (MQ=255)                                                      cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑tta                                                                                                                                        <  6:227541/96‑3 (MQ=255)                                                         cGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTTAT‑GAACGATATCGAGATCGCCAAGGGTAAg                                                                                <  7:28652/148‑1 (MQ=255)                                                                   tcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGAAAGGCGCATAAgg                                                                      >  7:116586/1‑148 (MQ=255)                                                                     gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                     <  1:226343/147‑1 (MQ=255)                                                                       tCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                  >  1:131856/1‑148 (MQ=255)                                                                         cTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAc                                                               <  5:146895/149‑1 (MQ=255)                                                                                gCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCa                                                        >  4:364704/1‑149 (MQ=255)                                                                                        aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                 <  7:301876/148‑1 (MQ=255)                                                                                          gCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCa                                                                    <  8:166505/127‑1 (MQ=255)                                                                                          gCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCa                                                                    >  7:166505/1‑127 (MQ=255)                                                                                           cACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGc                                             >  5:292262/1‑149 (MQ=255)                                                                                              cgggTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAGT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATCTCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGc                                          <  1:328974/147‑1 (MQ=255)                                                                                               tGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGcc                                         >  3:63408/1‑149 (MQ=255)                                                                                                 gTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                       <  4:63408/149‑1 (MQ=255)                                                                                                 gTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCt                                        <  8:77035/148‑1 (MQ=255)                                                                                                                           ccATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGATAGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGCCGCGGTCATTTGTGGGAGCGGCATTGCGTCGCGAATGGGCTGCAAagcag              >  1:113759/1‑149 (MQ=255)                                                                                                                           ccATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACTCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagcagc             >  3:359445/1‑149 (MQ=255)                                                                                                                           ccATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagcagc             >  4:221581/1‑149 (MQ=255)                                                                                                                           ccATAAATTCCTGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagcagc             >  2:295921/1‑149 (MQ=255)                                                                                                                                      tGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGc                        <  5:352162/127‑1 (MQ=255)                                                                                                                                      tGGGGGCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGCCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGc                        >  6:352162/1‑127 (MQ=255)                                                                                                                                        ggggCAGCCGATT‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt  <  3:364704/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                          AACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGA‑T‑ACCTGTGCGACATGGCCAAGGCGCTGAT‑GAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTT  >  NC_002947/335979‑336259

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