Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	304,439	A→G	13.3%	V111A (GTG→GCG)	PP_0249 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	304,439	0	A	G	13.3%	69.0 / 3.7	30	V111A (GTG→GCG)	PP_0249	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base G (2/2);  total (18/12)
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.12e-01

TGAACTCCACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGCCACGCCGAGGTCA  >  NC_002947/304302‑304576                                                                                                                                          |                                                                                                                                          tGAACTCCACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCa                                                                                                                               <  1:321392/150‑1 (MQ=255)      tCCACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTgc                                                                                                                          >  2:266077/1‑150 (MQ=255)        cACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAGTTGCGCCACCgc                                                                                                                                    >  2:78671/1‑137 (MQ=255)        cACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAGTTGCGCCACCgc                                                                                                                                    <  1:78671/138‑2 (MQ=255)         aCCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCa                                                                                                                       >  2:550313/1‑150 (MQ=255)          ccTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAg                                                                                                                      >  2:325235/1‑150 (MQ=255)                acacacCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGc                                                                                                                >  2:358073/1‑150 (MQ=255)                     ccATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCGTGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAAcc                                                                                                           >  1:13702/1‑150 (MQ=255)                      cATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAGTTGCGCCACCgc                                                                                                                                    >  1:124539/1‑123 (MQ=255)                      cATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAGTTGCGCCACCgc                                                                                                                                    <  2:124539/124‑2 (MQ=255)                           cccAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGc                                                                                                     >  1:146160/1‑150 (MQ=255)                                  ccAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGt                                                                                              >  1:396059/1‑150 (MQ=255)                                          gTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCg                                                                                      >  1:342644/1‑150 (MQ=255)                                                           ggTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGt                                                                     <  1:358073/150‑1 (MQ=255)                                                             tgtGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCCACAGTGCCACGCAATGGCAGTGGCGCAACTCCACCACCGCGCCCAGGTGAGCAGTGAAGCTAACCCACGGCTGGCCGTCGCGCTCAAAGGGTTCCACTTCGGttg                                                                   <  2:146160/148‑3 (MQ=255)                                                                  gCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGc                                                              <  1:550313/150‑1 (MQ=255)                                                                       gCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCt                                                         >  1:173044/1‑150 (MQ=255)                                                                                  tctTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTcc                                                                             >  1:280047/1‑119 (MQ=255)                                                                                  tctTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTcc                                                                             <  2:280047/119‑1 (MQ=255)                                                                                     tGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCAc                                           <  1:235287/150‑1 (MQ=255)                                                                                      gATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCAcc                                          >  2:490925/1‑150 (MQ=255)                                                                                         ggTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGgt                                       >  1:484989/1‑150 (MQ=255)                                                                                           ttttCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTg                                                                    >  1:313623/1‑119 (MQ=255)                                                                                           ttttCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGTGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTg                                                                    <  2:313623/119‑1 (MQ=255)                                                                                                            ggTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCGGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGc                    >  1:382237/1‑150 (MQ=255)                                                                                                                atctgCAGCGGCGGTTGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCgg                <  2:342644/146‑1 (MQ=255)                                                                                                                  tCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGcc              >  1:408120/1‑150 (MQ=255)                                                                                                                  tCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGcc              >  1:152145/1‑150 (MQ=255)                                                                                                                          gcggTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGCCACGCCGAg      <  1:325235/150‑1 (MQ=255)                                                                                                                              tGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGCCACGCCGAGGTca  <  2:396059/150‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          TGAACTCCACCTCCCACACACCATCACCCAGCGCCAGTGACGTGCGGATGACGTAGCGGGTCTGGGCGTGGCCGTTGGAGCTCTTGATGGTTTTCATCGTCACCAACGGTGCCTCGCAGCGGCGGTGGCGCAACTGCACCACCGAGCCCAGGTGCGCAGTGAAGCGAACCCAGGGCTGGCCGTCGCGCTCGAAGGGTTCCACTTCGGTGGCATGCAGGCTCGACGTGCTGGCACCGGTGTCGATCTTGGCGCGCAGCCCGGCCACGCCGAGGTCA  >  NC_002947/304302‑304576

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