Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,295,419	A→G	23.1%	intergenic (‑3/‑347)	opdN ← / → PP_3765	outer‑membrane porin D/H‑NS family protein MvaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,295,419	0	A	G	23.1%	28.6 / 2.8	17	intergenic (‑3/‑347)	opdN/PP_3765	outer‑membrane porin D/H‑NS family protein MvaT
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base G (2/2);  total (8/9)
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

ATAGACATTGCGTGCTTCGATCTTGGCACTGGCGTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCATTGGACCATGGTCT  >  NC_002947/4295273‑4295546                                                                                                                                                   |                                                                                                                                aTAGACATTGCGTGCTTCGATCTTGGCACTGGCGTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGaa                                                                                                                              <  1:462476/150‑1 (MQ=255)        ttGCGTGCTTCGATCTTGGCACTGGCGTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGAc                                                                                                                       <  2:614253/150‑1 (MQ=255)            gTGCTTCGATCTTGGCACTGGCGTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCGGAATGAAGCCGAGt                                                                                                                   >  1:74790/1‑150 (MQ=255)                                  gTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCGGAATGAAGCCGAGTTAACAACCTGTGGCACGTGTTg                                                                                             <  2:74790/150‑1 (MQ=255)                                                        tGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCGGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTcc                                                                       >  1:624314/1‑150 (MQ=255)                                                                             aaTAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTg                                                                                 >  1:792072/1‑119 (MQ=255)                                                                             aaTAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTg                                                                                 <  2:792072/119‑1 (MQ=255)                                                                                aGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGa                                               >  2:350645/1‑150 (MQ=255)                                                                                          cGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATcgcg                                     >  1:540708/1‑150 (MQ=255)                                                                                                  cGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGct                             <  1:227299/150‑1 (MQ=255)                                                                                                             ggCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCGGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGc                  <  2:624314/150‑1 (MQ=255)                                                                                                                gAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCAt               >  1:314751/1‑150 (MQ=255)                                                                                                                    gggaaaggggaaaACACGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCATTGGa           <  1:242362/150‑1 (MQ=255)                                                                                                                     ggaaaggggaaaACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCt                                                 >  2:308509/1‑111 (MQ=255)                                                                                                                     ggaaaggggaaaACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCt                                                 <  1:308509/111‑1 (MQ=255)                                                                                                                       aaaggggaaaACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCATTGGACCa        <  2:314751/150‑1 (MQ=255)                                                                                                                             gaaaACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCATTGGACCATGGTCt  >  2:56316/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                ATAGACATTGCGTGCTTCGATCTTGGCACTGGCGTCTTCAAGAAACCCGGCAGCCTGGGCGTTGACACCGGCACACAATAGAAACAGGCCGGTAACACGGCGCCCGGGGGCGAGGGGGAAAGGGGAAAACATGCGAGGAAACATCCAGAATGAAGACGAGTTAACAACCTGTGGCACGTGTTGGCCACGGTACTGCGGCAGGTCCGAACATGAGTCGGCCATCATCTGAAATGATCGCGCAAAGGCTCTGGCGCGGGCCATTGGACCATGGTCT  >  NC_002947/4295273‑4295546

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