Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,931,069	C→T	E9K (GAA→AAA)	ravA ←	regulatory ATPase RavA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,931,069	0	C	T	94.1%	54.0 / NA	17	E9K (GAA→AAA)	ravA	regulatory ATPase RavA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/7);  minor base G (0/1);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACGAAAGAC  >  NC_000913/3930823‑3931309                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gcgcCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCt                                                                                                                                                                                                                                              <  1:231690/251‑1 (MQ=255)                  ggACCAAAAACTTCTTCCGGCGTGGAAAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCTGGCGCTCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCCACGTTCATAACGCCCCTTTTCCAACGAACTGCTCAGGCGGGCAATTCTTTTCGCTAATAAATGAGGGttcgg                                                                                                                                                                                                                             >  2:68769/1‑247 (MQ=255)                            cttcttCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGa                                                                                                                                                                                                                     >  1:236780/1‑249 (MQ=255)                                                               tCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGTCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCGTTCGTCAAAAGTTCTGCGTAAATTg                                                                                                                                                                                 <  1:242891/250‑1 (MQ=255)                                                                         gggCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGAc                                                                                                                                                                      <  1:120942/251‑1 (MQ=255)                                                                                                     ccGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGtctc                                                                                                                                          <  1:6943/251‑1 (MQ=255)                                                                                                              cAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACAAAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTAAAAgggggg                                                                                                                                 >  1:290262/1‑251 (MQ=255)                                                                                                                                                    cacTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCAttt                                                                                            >  1:61968/1‑250 (MQ=255)                                                                                                                                                                    tGCCGCTAATAAACACAAGCGGCTGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCCGGCGGGAAAATCTTTTTGCTAATTAATGAGGGTGAGCCCTAATGAAGTGGCGTCCTTTCGGCAAAAGTTCTGCGTAAATTGCGAATATAGACGTTTCTTGCTGGGGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGCGTACGATTTTGcc                                                                           >  2:213027/1‑251 (MQ=255)                                                                                                                                                                    tGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATTTGGCCGCATTTAGGATTACGATTTTGcc                                                                           >  2:213037/1‑251 (MQ=255)                                                                                                                                                                              aaaCACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTg                                                                 <  2:331378/251‑1 (MQ=255)                                                                                                                                                                                 cacaAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCAt                                                              >  1:51959/1‑251 (MQ=255)                                                                                                                                                                                            ggCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTGCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTGTCGTCAAAAGTTCTGCGGAAATTGCGAGTATAGACGTTTCTGGCTGGTGGCTAAAATAGTCTCAAATGGGGGGTATTTTTCGTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCTTGCATACTGTGCGCtt                                                   <  2:354922/251‑1 (MQ=255)                                                                                                                                                                                                     aCGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTgg                                           <  1:168788/250‑1 (MQ=255)                                                                                                                                                                                                                                       aGGCGGGAAATTCTTTTCGCTAATAAATGAGGATGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGAATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACg        >  1:227930/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACGa       >  1:161796/1‑250 (MQ=255)                                                                                                                                                                                                                                             gAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACGAAAGAc  <  1:62274/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 GCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACGAAAGAC  >  NC_000913/3930823‑3931309

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