Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	435,813	T→A	D60E (GAT→GAA)	thiL →	thiamine monophosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	435,813	0	T	A	96.0%	78.3 / ‑3.6	25	D60E (GAT→GAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (7/17);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GAAGTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAG  >  NC_000913/435683‑435936                                                                                                                                   |                                                                                                                            gaagTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGAt                                                                                                                    <  1:2017205/140‑1 (MQ=255) gaagTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGAt                                                                                                                    <  1:2268679/140‑1 (MQ=255)  aagTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATc                                                                                                                   <  2:2051938/140‑1 (MQ=255)    gTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTg                                                                                                                 <  2:2021869/140‑1 (MQ=255)      tctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGc                                                                                                               <  2:1720065/140‑1 (MQ=255)      tctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGc                                                                                                               >  1:2124664/1‑140 (MQ=255)      tctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGc                                                                                                               >  2:1672021/1‑140 (MQ=255)       ctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCt                                                                                                              >  2:1867457/1‑140 (MQ=255)       ctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCt                                                                                                              >  1:1793289/1‑140 (MQ=255)       ctcGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGCAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCt                                                                                                              >  1:2066084/1‑140 (MQ=255)            cTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATaa                                                                                                         <  2:1730916/140‑1 (MQ=255)              tGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAg                                                                                                       <  2:1683700/140‑1 (MQ=255)                    cGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTgg                                                                                                 <  1:2044987/140‑1 (MQ=255)                        cTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGt                                                                                             <  1:1743612/140‑1 (MQ=255)                                 ggCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAg                                                                                    <  1:1998946/140‑1 (MQ=255)                                                     cAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAGCCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGACCCTAAGAGATCTGGCAGCGATggggg                                                                <  2:2200617/140‑1 (MQ=255)                                                          tCCCCGAGAAACAGACCCAGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGAt                                                           <  2:1765304/140‑1 (MQ=255)                                                           ccccGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGGTGGGGGCCGATc                                                          >  2:2221610/1‑140 (MQ=255)                                                                          ccTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGc                                           <  2:1793289/140‑1 (MQ=255)                                                                                         tGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTAc                            <  2:1855732/140‑1 (MQ=255)                                                                                         tGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTAc                            <  2:2090801/140‑1 (MQ=255)                                                                                           aTACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCg                          >  1:2254314/1‑140 (MQ=255)                                                                                                cTGGTGGCGGGTAACCATTTCCTCCCTGATATCCAACCTGCTGATCTGGCTTATAAAGCACTGGCGGGGAAACCAAACGATCTGGGAGCGCTGGGGGCCCGACCGGCCTGGCGCGCCCCGGCCTTAACCCTAACGGGCgg                     >  2:2212261/1‑139 (MQ=255)                                                                                                             aCCATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGGGAACCTAAACGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCCGGCGCCGGCATTAACCTTACCCGGCGGAGACGAAGCgggg        >  2:1630391/1‑140 (MQ=255)                                                                                                              ccATTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGc       <  1:2125005/140‑1 (MQ=255)                                                                                                                aTTTCCTCCCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTCCCGGACGTAGACGAAGCGTGGCtt     <  2:2124664/140‑1 (MQ=255)                                                                                                                   tcctccCTGATATCGAACCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAg  >  2:2275628/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                            GAAGTTCTCGTCTTGATGTCGAACTGGGCATCGGCGACGATTGCGCACTTCTCAATATCCCCGAGAAACAGACCCTGGCGATCAGCACTGATACGCTGGTGGCGGGTAACCATTTCCTCCCTGATATCGATCCTGCTGATCTGGCTTATAAAGCACTGGCGGTGAACCTAAGCGATCTGGCAGCGATGGGGGCCGATCCGGCCTGGCTGACGCTGGCATTAACCTTACCGGACGTAGACGAAGCGTGGCTTGAG  >  NC_000913/435683‑435936

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