Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,129,232	A→G	20.9%	intergenic (+38/+29)	PP_5376 → / ← PP_5377	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,129,232	0	A	G	20.9%	57.0 / 6.6	28	intergenic (+38/+29)	PP_5376/PP_5377	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/15);  new base G (3/3);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCAT  >  NC_002947/6129085‑6129374                                                                                                                                                    |                                                                                                                                               ggggAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAgc                                                                                                                                              >  2:294421/1‑150 (MQ=255)           tgACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCg                                                                                                                                    <  1:210826/150‑1 (MQ=255)             aCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGtt                                                                                                                                  <  2:366020/150‑1 (MQ=255)             aCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGtt                                                                                                                                  >  2:177740/1‑150 (MQ=255)               gTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTgg                                                                                                                                <  2:424007/150‑1 (MQ=255)                                    tATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGcc                                                                                                           >  1:89034/1‑150 (MQ=255)                                     aTGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCt                                                                                                          <  2:280492/150‑1 (MQ=255)                                             tCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATActc                                                                                                  <  1:88830/150‑1 (MQ=255)                                              cGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATActct                                                                                                 <  2:344089/150‑1 (MQ=255)                                              cGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATActct                                                                                                 <  1:22991/150‑1 (MQ=255)                                               gACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtttg                                                                                                                                >  1:114206/1‑116 (MQ=37)                                               gACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtttg                                                                                                                                <  2:114206/118‑3 (MQ=37)                                                      aaGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCAc                                                                                         >  1:346594/1‑150 (MQ=255)                                                            aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATc                                                                                   >  1:337420/1‑150 (MQ=255)                                                            aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATc                                                                                   <  1:294421/150‑1 (MQ=255)                                                            aaCATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATc                                                                                   <  1:100907/150‑1 (MQ=255)                                                                   tcGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGCTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCg                                                                            >  2:5628/1‑150 (MQ=255)                                                                                         cAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCt                                                      <  1:189803/150‑1 (MQ=255)                                                                                            ttGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGTTTGGa                                                                                                                              >  2:277904/1‑75 (MQ=2)                                                                                            ttGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCCAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGTTTGGa                                                                                                                              <  1:277904/75‑1 (MQ=11)                                                                                                      cAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                         <  2:288649/150‑1 (MQ=255)                                                                                                         aaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGac                                      <  2:337420/150‑1 (MQ=255)                                                                                                              aGCTCTGAAGGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCAt                                 <  2:392113/150‑1 (MQ=255)                                                                                                                 tctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCa                              <  2:346594/150‑1 (MQ=255)                                                                                                                  ctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCAc                             <  2:89034/150‑1 (MQ=255)                                                                                                                           gTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAcc                    <  1:5628/150‑1 (MQ=255)                                                                                                                                     aaaCGGCGCTTTTAAGGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCt          >  1:331584/1‑150 (MQ=255)                                                                                                                                             ccTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCAt  >  1:151385/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                               GGGGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCAT  >  NC_002947/6129085‑6129374

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