Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,822	T→C	D363D (GAT→GAC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,822	0	T	C	95.8%	68.5 / ‑6.4	24	D363D (GAT→GAC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/15);  minor base G (0/1);  total (8/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTA  >  NC_000913/4609‑5059                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               gactcaggggcggTTATCAAACGCGATGGACGTGAATCAGCCGAACAACTGGCGGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGGCAATGCGTGAGTGAAAAGAGCGGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGGATCGGGCGCTGCGTGAGCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCa                                                                                                                                                                                                           <  2:183573/239‑1 (MQ=255)  aCTCAGGCGACGTTATCCAACGAGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCAc                                                                                                                                                                                                          >  1:213614/1‑251 (MQ=255)                aTCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGCAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGACTATGGCTTGTTCCTCGGCACCGCGCATCCGGCgc                                                                                                                                                                                            >  1:99283/1‑250 (MQ=255)                            ggACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAagag                                                                                                                                                                                 <  1:69552/250‑1 (MQ=255)                                              caacgggcggcggggggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTGGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCg                                                                                                                                                              <  2:133867/236‑1 (MQ=255)                                              caacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCg                                                                                                                                                              <  1:213681/251‑1 (MQ=255)                                               aacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                             >  1:205149/1‑251 (MQ=255)                                                   ggCCGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                          <  1:68874/250‑1 (MQ=255)                                                             ggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAATAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                                                <  2:240696/250‑1 (MQ=255)                                                              gAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAACGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCTGCACCGCGCATCCGGCGACATTTAAAGAGAGCGTGGAAGCGTTTCTCGGTGAAACGTTGAAT‑Ct                                                                                                                                             >  1:193923/1‑251 (MQ=255)                                                              gAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑Ct                                                                                                                                             <  2:205149/251‑1 (MQ=255)                                                               aagatgtgttCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAAGCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTg                                                                                                                                            <  2:125943/245‑1 (MQ=255)                                                                  agTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGcc                                                                                                                                          >  2:215428/1‑250 (MQ=255)                                                                        cccgccgCAAAATCTGGCAACGGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAagag                                                                                                                                   <  2:40314/250‑1 (MQ=255)                                                                                               aaaGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTAcc                                                                                                            <  1:858/251‑1 (MQ=255)                                                                                               aaaGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATACGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑ATGCCAAAAGAGCTGGCAGAACGTGCTGATTTAcc                                                                                                            >  2:214634/1‑251 (MQ=255)                                                                                                aaGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGAGCTGATTTAccc                                                                                                           >  2:59013/1‑251 (MQ=255)                                                                                                                                aCCACGCAACAGACGATGCGTGCGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGCTTCTCGATGACACGTTGGAT‑CTTCCAACAGAGCTGGCAGAAAGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTtgc                                                                           >  1:46569/1‑251 (MQ=255)                                                                                                                                   ccGCAACAGACAATGCGTGAGTTAAAATAACGGGTCTACCCCTCGGAGCCGCACGCTGCTGTGGCGTATCGTTCGTTGCGTGACCAGTTGAATCCAGGCGAATATGGCCTTTTCCTCGGCACCGCGCTTCCGGGGAAATTTAAAGAGAGCTGGGAAGGGATTCCCGGTGAAACTTTGGAT‑CTGCCAAAAGAGCTGTCAGAACTGGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTtgctgc                                                                        <  1:109145/250‑1 (MQ=255)                                                                                                                                                 gCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCGGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTTCGTTGCGTAAATtg                                                           <  1:215428/250‑1 (MQ=255)                                                                                                                                                            aaGAACTGGGCTACACCTCGGAGCCTCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGTGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcat                                               <  1:107474/251‑1 (MQ=255)                                                                                                                                                                                     cccctcttccgtggtgtttcttGCGCTGCTTGACCAGTAGAATCCAGGCGAATATGGCTTTTTCCTCGGCTCCGGGCAGCGTACGAAAGTTAAGGGGGCCGTGGAAGCCATTCTCGGTGAAGC‑TATGATACTGCCAAAGAAGCTGGAAGATCGTGCTTATTTACCTTTGCTTTCACATAATTCGCCCGCCGATTTTGCTGCGTTGCGTAAATGGATGATGAATCATCAGTAAAATCTATTCATTATCTCaa                      <  1:127321/230‑1 (MQ=255)                                                                                                                                                                                             cggTAGCGTATCGGGCGCTGCGTGACCAGTGGAAACCAGGAGAATAGGGCTTGTTCCTGGGCACCGCGAATCCGGAGAAATTAAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTGGCTTTCACATAATCGGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCg              <  2:140074/249‑1 (MQ=255)                                                                                                                                                                                                         gTGCGCTGTGTGACCAGTTGAATCCAGGCAAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTa  <  1:134683/251‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               GACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAT‑CTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTA  >  NC_000913/4609‑5059

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