Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,270,523	C→A	L96M (CTG→ATG)	yjbQ →	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,270,523	0	C	A	100.0%	92.1 / NA	27	L96M (CTG→ATG)	yjbQ	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/13);  total (14/13)

GACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGAC  >  NC_000913/4270384‑4270662                                                                                                                                            |                                                                                                                                            gACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCCTTTTTTCCTCAGAACCGTTCCCGAAAACGCAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATAAAATCCTCAATGATGGGaaaa                                                                                                                                     >  7:45797/1‑146 (MQ=255)         aTGAGAACTTCGATCCAAACTTACGCCACTTCATGGAGCGTTTTTTCCTCCGCACAGGTCCCGACAACGGAAATTATTAGCATGACTATGAGGGAGAATACGATATGCCTTCTCATATCAAATACTCAATGATGGTAACATAGCTTGTa                                                                                                                            <  8:45797/149‑1 (MQ=255)                cTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCgg                                                                                                                     <  7:7859/149‑1 (MQ=255)                 tGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCg                                                                                                                      <  8:114907/147‑1 (MQ=255)                 tGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCg                                                                                                                      >  7:114907/1‑147 (MQ=255)                                  ccACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTAtt                                                                                                    >  8:48807/1‑148 (MQ=255)                                   cACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTAtt                                                                                                    >  7:129211/1‑147 (MQ=255)                                            gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCa                                                                                          <  1:162868/148‑1 (MQ=255)                                            gAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCa                                                                                          >  2:162868/1‑148 (MQ=255)                                                         tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGt                                                                                                                    <  8:2580/109‑1 (MQ=255)                                                         tccGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGt                                                                                                                    >  7:2580/1‑109 (MQ=255)                                                            gCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTgg                                                                         <  8:129211/149‑1 (MQ=255)                                                                 gTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTggg                                                                    <  7:48807/149‑1 (MQ=255)                                                                              gAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATc                                                       >  8:104360/1‑149 (MQ=255)                                                                                        gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCa                                                                                                  <  1:170215/96‑1 (MQ=255)                                                                                        gCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCa                                                                                                  >  2:170215/1‑96 (MQ=255)                                                                                                 tGAGGGCGCATACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTAGTCAGACCGGCACCTGGCAAGGCAGTTGGCTGGGGGCACATCGCCTCCACGGCGGATCGCGTCg                                      >  6:101128/1‑147 (MQ=255)                                                                                                    gggAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGca                                    >  4:24926/1‑146 (MQ=255)                                                                                                    gggAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGca                                    <  3:24926/146‑1 (MQ=255)                                                                                                     ggAGCAGACGATATGCCTTCTCATATCAGATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTgg                                                                                     >  6:4342/1‑96 (MQ=255)                                                                                                     ggAGCAGACGATATGCCTTCTCATATCACATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTgg                                                                                     <  5:4342/96‑1 (MQ=255)                                                                                                          aGACGATATGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa                           <  7:104360/149‑1 (MQ=255)                                                                                                                tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCgg                                                                                                                     <  6:67316/53‑1 (MQ=255)                                                                                                                tatGCCTTCTCATATCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCgg                                                                                                                     >  5:67316/1‑53 (MQ=255)                                                                                                                             atCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTaaa         >  3:107346/1‑148 (MQ=255)                                                                                                                              tCAAATCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTaaaa        >  3:45793/1‑148 (MQ=255)                                                                                                                                   tCCTCAATGATGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGAc  <  5:101128/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            GACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTACAAGGGGAGTAAAAAATGAC  >  NC_000913/4270384‑4270662

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