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	96.2%	55.2 / NA	27	E9K (GAA→AAA)	ravA	regulatory ATPase RavA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (11/14);  minor base G (0/1);  total (12/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 = 8.95e-01

GAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATT  >  NC_000913/3930829‑3931297                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      gAATATAAAGGGGCCCAAAAACGTCTTCCGGCGTGGAGAAGCGGGTCATCAGATACTCAAACGCGCGGGCATTCTTAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACAGTTCAGGCGGGAAATTCTTTTCGCTAATaaa                                                                                                                                                                                                                             <  1:17315/251‑1 (MQ=255)             gACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGcc                                                                                                                                                                                                                 >  1:4822/1‑251 (MQ=255)                   aaacttttcccgggggggAGAAGCGGGTCACAAGATATTAAAACGCCCGGGCATTCTGAAAGGCAAATTTTAAGCGCCGGGCGAACAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACTCTTCCACCACTTAATGCCGCTAATAACCACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCGTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTGCGCTAATAAATGAGGGCGAGCCATAATg                                                                                                                                                                                                           <  2:294729/235‑1 (MQ=255)                                   ggAGAAGCGGGTCATCAGATATTCAACCGAGCGGGCATTCTGAAAGGCGAATTTTAAGAGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCACCTTTTCCAACGAACTGATCAGGCGGGAAATTCATTTCGATAATAAATGAGGGAGAGTCATAATGAAG‑TGGCGTCCTTTc                                                                                                                                                                                           >  1:284951/1‑250 (MQ=255)                                      gaAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGGGAGCCATAATGAAG‑TGGCGTCCTTTCGTc                                                                                                                                                                                        >  1:455522/1‑250 (MQ=255)                                                                       aTTCTGAAAGGTGAATTTTAAGCGCCGGGCGATCACACTTTTGGCAATACCTCGCGGGCCAAGGAGGAACATACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGAGTGGCTACGATCATAAAACCCCATTTCCAACGAACTGATCAGGCGGGACATTCTTTTTGCTAATAAATGAGGGGGTGCCACAATGAAG‑TAGCGTCCTTACGTCAAAAGTTCGGCGTAAATTGCTAGTATAGAAGttt                                                                                                                                                      >  1:197385/1‑251 (MQ=255)                                                                       aTTCTGAAAGGTGAATGTTAAACGACGTGCGATTAAAATTTTGGCAATACCGGGCGGGCCAAGGAGGAACACTCTTTCACCACCTAATGCCGCTAATAAACACAAGCGGATGGCGTGTATACGTTCATAAAGCCTCTTTTCCAACTAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATGGCGAGTATAGTCGttt                                                                                                                                                      <  2:9684/251‑1 (MQ=255)                                                                              aaGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAAACGTTTCTTGCtg                                                                                                                                               >  1:328570/1‑251 (MQ=255)                                                                               aGGCGAATTTTAAGCGCCTGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCtgg                                                                                                                                              <  2:455522/251‑1 (MQ=255)                                                                               aGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCtgg                                                                                                                                              <  1:245326/251‑1 (MQ=255)                                                                                ggCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCATTCCCAGACGGGCCAAGGAGGAACACACTTACACAACTTAATGCAGCTAATAAACACACGCGGATAGCGTAGCTAGCATCATAAAGGCCCTTTTCCAACGAAATGCTCAGGCGGGAAATTCTTTTCTCTAATAATAGAGGGGCGGCCATAATGAAGGTAGGGTC‑TTTTGTCAAAAGTTCTGCTTAAAATGcttagtgtgaggttttttgctggg                                                                                                                                             >  2:54250/1‑228 (MQ=255)                                                                                           agggcccggCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGACGGCGTGGCTACGTTCATAAAGCCCCTTATCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGGGAGCCATAATGAAG‑GGGCGTCCTTTCGTCAAAAGTTTTGCGTAAATTGCGAGTATAGACGTTTCTTGCGGGTGGCTGAAATAg                                                                                                                                  <  2:39869/244‑1 (MQ=255)                                                                                                 gggCGAACAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGAGTGGCAACGTTCATAAAGACCCTTTTCCAACGAACTGATCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGGAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGGCGATTCTTGATGGGTGCTAAAATAgactaa                                                                                                                             >  1:193129/1‑245 (MQ=255)                                                                                                                      tACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCt                                                                                                                                                    <  2:300474/206‑1 (MQ=255)                                                                                                                      tACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCt                                                                                                                                                    >  1:300474/1‑206 (MQ=255)                                                                                                                         cTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTaaa                                                                                                                                      >  1:78268/1‑217 (MQ=255)                                                                                                                         cTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTaaa                                                                                                                                      <  2:78268/217‑1 (MQ=255)                                                                                                                                                        acTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGAt                                                                     <  1:367441/251‑1 (MQ=255)                                                                                                                                                          ttAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGACTCAAAGGGGGG‑TATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTTGGAGTACGAtttt                                                                  >  2:110629/1‑251 (MQ=255)                                                                                                                                                            aaTGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGG‑TATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGc                                                                <  1:132364/251‑1 (MQ=255)                                                                                                                                                                              aaGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACt                                               >  1:237656/1‑251 (MQ=255)                                                                                                                                                                                               aCGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATGATGAAG‑TGGCGTCATTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGGTAATGGGGCAGCATTTAGGAGTACGATTTTGCCGTTAATCGGGCGTGCTGTGAGCTTTTTTGTggg                              >  1:83031/1‑251 (MQ=255)                                                                                                                                                                                                    cATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAg                         <  2:130306/251‑1 (MQ=255)                                                                                                                                                                                                                  tCCAACGAACTTCTCAGGCGGGAAATTCTTGTCGCTAATTAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGTCGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGcacac            <  1:190721/250‑1 (MQ=255)                                                                                                                                                                                                                     aaCGAACTGCTAAGGCGGGAAATTCTTTTCGCGAATAAAGGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACAttt        <  2:126464/251‑1 (MQ=255)                                                                                                                                                                                                                          acttcTCAGGAGGGAACGTCTTTTCGCTAATAAATTAGGGTGGGCCATAATGAAG‑TGGCGTCCTTTCGTTAAAAGATCTGTGGAAATTGCGAGTTGAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGTGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATGTGGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCatat   <  2:273525/247‑1 (MQ=255)                                                                                                                                                                                                                             gCTCAGGCGGGAAATTCTTTTCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATAtt  <  1:220205/249‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      GAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCACCACTTAATGCCGCTAATAAACACAAGCGGATGGCGTGGCTACGTTCATAAAGCCCCTTTTCCAACGAACTGCTCAGGCGGGAAATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAG‑TGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATT  >  NC_000913/3930829‑3931297

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