Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,877,897	A→G	14.2%	M7V (ATG→GTG)	PP_1684 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,877,897	0	A	G	14.2%	64.4 / 3.2	28	M7V (ATG→GTG)	PP_1684	transporter
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (2/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AGAAAGAGCTGGTTGAACGGCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGATTTTGGCTGGGGGCGTGGGGTGTTT  >  NC_002947/1877761‑1878037                                                                                                                                         |                                                                                                                                             aGAAAGAGCTGGTTGAACGGCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGgtgt                                                                                                                                  <  4:249766/149‑1 (MQ=255)            gTTGAACGGCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGc                                                                                                                       <  6:129714/149‑1 (MQ=255)                    gCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGc                                                                                                                       <  8:281364/141‑1 (MQ=255)                    gCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGc                                                                                                                       >  7:281364/1‑141 (MQ=255)                      tCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGgtgt                                                                                                             >  6:13054/1‑149 (MQ=255)                                      ggATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCtt                                                                                             <  5:13054/149‑1 (MQ=255)                                      ggATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCt                                                                                              <  8:277229/148‑1 (MQ=255)                                       gATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCtt                                                                                             >  7:168179/1‑148 (MQ=255)                                       gATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTg                                                                                            <  6:147952/149‑1 (MQ=255)                                                                aacTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTac                                                                   >  2:147343/1‑149 (MQ=255)                                                                       cTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTacggcacgg                                                            >  6:279967/1‑149 (MQ=255)                                                                               aTACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGcc                                                    >  1:11983/1‑149 (MQ=255)                                                                               aTACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGc                                                     <  1:147343/148‑1 (MQ=255)                                                                                      tataAGCGGGTATATACACGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCt                                             >  6:36639/1‑149 (MQ=255)                                                                                            cGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGcc                                       <  5:279967/149‑1 (MQ=255)                                                                                             gggTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGcc                                       >  1:280431/1‑148 (MQ=255)                                                                                             gggTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCg                                      <  2:280431/149‑1 (MQ=255)                                                                                               gTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGAt                                    <  2:11983/149‑1 (MQ=255)                                                                                                 atatcCCCGCCAACATTGCGATGTGCTGGAGATAACGACAGTGACGTCGGGGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCGGGGCGTACGGCACGGCGTTGGCCTGTTTCTGGCGCCGag                                    >  4:96815/1‑146 (MQ=255)                                                                                                 atataCCCGCAAACCTTGCGATGTGCTGGAGATAACGACAGTGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGAt                                    >  6:158119/1‑147 (MQ=255)                                                                                                  tataCCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGGGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGa                                  <  1:174478/148‑1 (MQ=255)                                                                                                    taCCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGt                                >  3:139418/1‑148 (MQ=255)                                                                                                      cccGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTAc                              <  4:139418/148‑1 (MQ=255)                                                                                                            aaCCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGAtttt                       <  5:36639/149‑1 (MQ=255)                                                                                                                   cgctgtgCTGGAGATAACGACCGTGACTTCGGTGTGGCGAACCAGCGGATGGGTGGTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGCGCCGATGAGTACTGATTTTGGCTgg                 <  3:96815/145‑1 (MQ=255)                                                                                                                   cGATGTGCTGGAGATAACGACAGTGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGATTTTGGCTgg                 <  5:158119/148‑1 (MQ=255)                                                                                                                            ggAGATAACGACAATGACTTCGGTGTGGTGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGATTTTGGCTGGGGGCGTgggg       >  1:179475/1‑149 (MQ=255)                                                                                                                                 tAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGATTTTGGCTGGGGGCGTGGGGTGttt  <  2:52839/149‑1 (MQ=255)                                                                                                                                         |                                                                                                                                             AGAAAGAGCTGGTTGAACGGCTCGGGGTGGGGCAGCGGGATGAACTGGTGCGACTGCTGGAACAACTGGCCTGATTTGATACGACTATAAGCGGGTATATACCCGCAAACCTTGCGATGTGCTGGAGATAACGACAATGACTTCGGTGTGGCGAACCAGCGGATGGGTGTTGGTGGGGGCGGCCTTGATACTGGCGCTGTCTCTGGGCGTACGGCACGGCTTTGGCCTGTTTCTGGAGCCGATGAGTACTGATTTTGGCTGGGGGCGTGGGGTGTTT  >  NC_002947/1877761‑1878037

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