Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	435,874	G→A	A81T (GCC→ACC)	thiL →	thiamine monophosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	435,874	0	G	A	91.7%	19.8 / NA	12	A81T (GCC→ACC)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (6/5);  minor base C (1/0);  total (7/5)
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

GTAAGAAGTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGC  >  NC_000913/435679‑436069                                                                                                                                                                                                    |                                                                                                                                                                                                    gtAAGAAGTTCTCGTCTTGATGTCGAACTGTGCAGCGGCGATGATTGCGCACTTCTTAATATCCCCGAGAAACAGTCCCTGGCGATCAGCACCGATACGCTGGTGGCGGGTAGCCAGTGCCTCTCTGATATCGATCTTGCTGTGCTGGGTTATAAATCTCTGGCGGTGAACCTAGGCGTTCTGGCAGCTATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTAACCTTan                                                                                                                                                                 <  1:59668/232‑2 (MQ=255)                                       cgacgaTTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTTAACCTAAGCGATCTGGCCGCGATGGGGACCCATCCGGCCTGGCTGACGCTTGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTTAGTCCTTccacgcccattttttttatcttctcc                                                                                                        >  1:402765/1‑226 (MQ=255)                                           gatgtcgcACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATTCGCTGGTGGCGGGTAACCCTTTCCTCCCTGATATCGATCCTTCTGATCTGGCTTATAAAGCCCTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGCCGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAatta                                                                                                    <  1:20820/246‑1 (MQ=255)                                                                                          cACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTccg                                                                                                                               >  1:179485/1‑175 (MQ=255)                                                                                                                 aCCATTTACTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGCCGAAGCGTGGCTTGAGt                                                                                                                                      >  1:387738/1‑147 (MQ=255)                                                                                                                        nctccCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCATGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCCATTATTACGATATGCAACTCATTGGCGGCGCTACCCCGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTc                       >  1:59670/2‑251 (MQ=255)                                                                                                                                   tCGATCCTGCTGAGCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCCTGGCTGACGCTggn                                                                                                                                                                            <  1:47583/91‑2 (MQ=255)                                                                                                                                       tCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACGTGGGGTATCCACGGCTTTGTTCCGATGGGACGAGcc         <  1:135087/250‑1 (MQ=255)                                                                                                                                             tgATCTGGCTTATAAAGCACTGGCGGTGAAACTAAGCCCTCTGGCAGCGCCGGGGACCGATACGGCCTCGCTAACGCTGGCATTACCCTTACCGGACGTAAACGAAGCGTGGCTTGAGTCCTTCAGCGACCGTTTGTTTGATCTTCTCAATTATTACGATATGCCAATCATTGGCGGCGATAACCCGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAAcgc  >  1:235723/1‑251 (MQ=255)                                                                                                                                                                cTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGCCCGCTCCGGCCTGGCTGACGCTGGAATTAACATTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTTAGCGGCAgtttgttt                                                                                                                  >  1:324886/1‑120 (MQ=255)                                                                                                                                                                       tGAACCTAAGCGATCTGGCNGCGATGGGGACCGAGCNGNCNTNGCTGAc                                                                                                                                                                                  <  2:478370/49‑1 (MQ=25)                                                                                                                                                                                         cAGCGATGGGGACCGATCCGGCCTGGCTGACGCTGGCATTACCCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATATTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGc                                                           >  1:434890/1‑150 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                    GTAAGAAGTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAGTCCTTCAGCGACAGTTTGTTTGATCTTCTCAATTATTACGATATGCAACTCATTGGCGGCGATACCACGCGTGGGCCATTATCAATGACGTTGGGTATCCACGGCTTTGTTCCGATGGGACGAGCCTTAACGC  >  NC_000913/435679‑436069

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