Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,035,788	G→A	14.8%	A132A (GCC→GCT)	PP_0894 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,035,788	0	G	A	14.8%	55.3 / 4.1	27	A132A (GCC→GCT)	PP_0894	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base A (2/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

ACATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGG  >  NC_002947/1035639‑1035935                                                                                                                                                      |                                                                                                                                                    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)                                                                                                                                                    cAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTgg  <  1:104810/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACATCGCGGTGGACCAGGCGCATGTGGATGGCAACCTGGTGACCGCACCGGCATGGCCCGCGCACCCGGCCTGGCTGGCGGCGTTCCTTAAAGTGCTGGGCACGCGTATTGTCTGACAGGCTTGATTGAGTTGGCAGCGGGTTGCTTCAGGCAACCTGCTCCACCCAGTCATTGAGGTTGTAGTAGTTGGTGACCCGGGCAATCTTGCCACAGTGAATATAGAAGAACGCCCCTGCCGGCAGCACATAGGTCTGGCCGTTTGCTGTCGGCAACCCTTCGTCATCGGCCAGGTATTGG  >  NC_002947/1035639‑1035935

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