Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	844,942	G→T	16.6%	P357H (CCC→CAC)	PP_0728 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	844,942	0	G	T	16.6%	75.0 / 3.1	29	P357H (CCC→CAC)	PP_0728	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (2/3);  total (14/15)
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

TGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGCAACCCT  >  NC_002947/844809‑845077                                                                                                                                      |                                                                                                                                        tGAAGTGCCGGAGCATGGTCTTCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGGTTGTCGCCGCCATTAGGGTGTTGAAGCGAGTTTCGCCCACGCCGATAGCTGCGGGa                                                                                                                                       <  1:142504/136‑1 (MQ=255)           gAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGTGAATGGCGAAct                                                                                                                             >  2:174076/1‑136 (MQ=255)           gAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGTGAATGGCGAAct                                                                                                                             >  2:20453/1‑136 (MQ=255)                                        cAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGtta                                                                                                >  2:237177/1‑135 (MQ=255)                                        cAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGtta                                                                                                >  2:18000/1‑135 (MQ=255)                                          ccACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTtctc                                                                                              <  1:160010/135‑1 (MQ=255)                                                           tCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGtga                                                                             >  1:73193/1‑136 (MQ=255)                                                            cGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGtgat                                                                            <  2:5415/136‑1 (MQ=255)                                                                gcAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGaaa                                                                        >  2:130336/1‑136 (MQ=255)                                                                gcAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGaaa                                                                        >  2:474710/1‑136 (MQ=255)                                                                            cagGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCgg                                                            >  2:447830/1‑136 (MQ=255)                                                                                     cgccgcCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTaaaa                                                   >  1:158026/1‑136 (MQ=255)                                                                                         gccAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATcc                                               >  2:117072/1‑136 (MQ=255)                                                                                                     ttGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCg                                   >  1:337936/1‑136 (MQ=255)                                                                                                        aaTCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGAt                                >  2:455295/1‑136 (MQ=255)                                                                                                         aTCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATc                               >  2:181886/1‑136 (MQ=255)                                                                                                             aGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCg                           <  1:447830/136‑1 (MQ=255)                                                                                                             aGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCg                           <  2:82088/136‑1 (MQ=255)                                                                                                                aTCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCAGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGttt                        >  2:427536/1‑136 (MQ=255)                                                                                                                             aTAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCAGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACaa           <  2:256303/136‑1 (MQ=255)                                                                                                                               aGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATg         <  1:474710/136‑1 (MQ=255)                                                                                                                               aGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATg         <  1:455295/136‑1 (MQ=255)                                                                                                                               aGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATg         <  1:8831/136‑1 (MQ=255)                                                                                                                               aGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATg         <  1:130336/136‑1 (MQ=255)                                                                                                                                gCTGCGTGTATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAATCCAGCAATGCGGTGGGTAAAAATCCAATCACAACCCGGATCAGCGTTTTCGATTAAGACAATGc        <  2:444298/136‑1 (MQ=255)                                                                                                                                gCTGCGTGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGc        <  1:174076/136‑1 (MQ=255)                                                                                                                                gCTGCGTGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGc        <  1:20453/136‑1 (MQ=255)                                                                                                                                gCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGc        <  2:361260/136‑1 (MQ=255)                                                                                                                                gCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGc        <  2:88080/136‑1 (MQ=255)                                                                                                                                     gggAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGCAAccc   <  2:225306/136‑1 (MQ=255)                                                                                                                                      ggAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGCAAcccg  >  1:190176/1‑135 (MQ=255)                                                                                                                                      |                                                                                                                                        TGAACTGCCTGAGCATGGTCATCCCAACCAAAATCATGCCAGCACCTGTATACGCCGCTCGGCGCAGTCGAGCAGCAGGCTTGTCGCCGCCAATAGGGTTTTGAATCGAGTATCGCCCACGCCGATAGCTGCGGGAATGGCGAACTCTCGGGCTCGGATCGACATATGCGAGTTCTCGCCTCCGTAGGCGGTGATGAAACCAGCAATGCGGTGGGTAAAAATCCAATCAAAACCCGGATCAGCGTTTTCGATTAAGACAATGCAACCCT  >  NC_002947/844809‑845077

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