Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,035,784	T→C	14.3%	intergenic (+30/+1)	PP_0893 → / ← PP_0894	ThiJ/PfpI family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,035,784	0	T	C	14.3%	60.3 / 3.0	28	intergenic (+30/+1)	PP_0893/PP_0894	ThiJ/PfpI family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base C (2/2);  total (14/14)
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

ATCGACATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTAT  >  NC_002947/1035635‑1035932                                                                                                                                                      |                                                                                                                                                     aTCGACATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCt                                                                                                                                                      >  2:746273/1‑150 (MQ=255)     aCATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCCCCCGGATTGTCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAg                                                                                                                                                  <  2:699950/150‑1 (MQ=255)     aCATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAg                                                                                                                                                  <  1:733734/150‑1 (MQ=255)           cgGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGACTGGCCTTAAAGTGCTGGGGACGCGTAGTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAAcc                                                                                                                                            >  1:270682/1‑150 (MQ=255)            ggTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCt                                                                                                                                           <  2:692419/150‑1 (MQ=255)             gTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTg                                                                                                                                          <  1:480488/150‑1 (MQ=255)                 aCCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTcc                                                                                                                                      <  2:137035/150‑1 (MQ=255)                   cAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCTAc                                                                                                                                    <  1:877342/150‑1 (MQ=255)                   cAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCAGGTTGCCTGAAGCAACCCGCt                                                                                                                                        <  1:44585/146‑1 (MQ=255)                   cAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCAGGTTGCCTGAAGCAACCCGCt                                                                                                                                        >  2:44585/1‑146 (MQ=255)                      gcgcATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCa                                                                                                                                 >  1:398958/1‑150 (MQ=255)                         cATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCGGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTc                                                                                                                              >  1:873658/1‑150 (MQ=255)                                            tGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTg                                                                                                           <  2:398958/150‑1 (MQ=255)                                                           ggCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATc                                                                                            <  2:843906/150‑1 (MQ=255)                                                                      ccGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTg                                                                                         >  1:426541/1‑142 (MQ=255)                                                                      ccGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTg                                                                                         <  2:426541/142‑1 (MQ=255)                                                                                      cgTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGcc                                                                 <  1:1133923/150‑1 (MQ=255)                                                                                       gTTCCTTAAAGTGCTGGGCACGCGGATTGTCTGACAGGCTTGAGTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGccc                                                                >  1:1118612/1‑150 (MQ=255)                                                                                        ttCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCCTGAAGCAAcccgc                                                                                                                                         >  2:942666/1‑73 (MQ=255)                                                                                        ttCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGACAGCGGGTTGCCTGAAGCAACCCGc                                                                                                                                         <  1:942666/76‑1 (MQ=38)                                                                                                     tGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCGACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGcaca                                                  <  1:746273/150‑1 (MQ=255)                                                                                                                       gACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGATGATTGTCGAACTGGTGCCACAGTGAATATAGAAGAACGACCCAGACCACATTACGTAGTTCTGGCCGTTTGCt                                >  1:327707/1‑150 (MQ=255)                                                                                                                       gACAGGCTTCATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCAGCTGCACCGAGTCACAACGATCGCAGTAATAAGGCGCCCGCGCTGATGTGCCTCGGTGTGGATAGATTAACGCCCCTGCCGGCAGCACCTAGGTCTGGCCGTTTGCt                                >  2:1041780/1‑150 (MQ=255)                                                                                                                          aGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTc                             <  1:121102/150‑1 (MQ=255)                                                                                                                                 attgaGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCAGTGAGGTTGTAGTAGTTGGTGACCCGGGCAAACTTGCCACAGTGAATATAGAAGACCGCCCCTGCCGGCCGCACATATGTCTGGCCGTTTGCGATGGGCAAcc                      >  1:566413/1‑150 (MQ=255)                                                                                                                                      gTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCg                 >  1:502455/1‑150 (MQ=255)                                                                                                                                       ttGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGt                >  2:573625/1‑150 (MQ=255)                                                                                                                                                     cTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCaatacc  >  1:1016961/1‑145 (MQ=255)                                                                                                                                                     cTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTAt  >  1:460490/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ATCGACATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTAT  >  NC_002947/1035635‑1035932

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