Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,634,286	T→C	L164L (TTA→CTA)	ytjC →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,634,286	0	T	C	100.0%	58.5 / NA	19	L164L (TTA→CTA)	ytjC	putative phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/7);  total (12/7)

GTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGATTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGTCGTTGAAACCGCAGGGGACATCTCG  >  NC_000913/4634156‑4634405                                                                                                                                   |                                                                                                                        gTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTa                                                                                                                       >  2:408551/1‑133 (MQ=255)           gAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGATGAGTACGATTTTCGGACTACTAGCATggg                                                                                                             >  1:143450/1‑133 (MQ=255)                cAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCagag                                                                                                        >  2:326297/1‑133 (MQ=255)                  gCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAgcg                                                                                                      <  2:346905/133‑1 (MQ=255)                           tCAATGCGGCACTGGAATCTTGCCGGGACTTACCGTAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACg                                                                                             >  2:59467/1‑133 (MQ=255)                               tGCGGCACTGGAATCTTGCCGGGATTTATCGCTGGGAAGCCTACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTg                                                                                         <  2:187833/133‑1 (MQ=255)                                    cACTGGAATCTTGCCTGGACTTACCGCACGGAGGCCGACCGTTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTaa                                                                                    <  2:154702/133‑1 (MQ=255)                                    cACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTaa                                                                                    <  2:57571/133‑1 (MQ=255)                                    cACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTaa                                                                                    <  1:386197/133‑1 (MQ=255)                                                gCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTc                                                                        >  2:272967/1‑133 (MQ=255)                                                              cAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGTGTAACTGTTCTTTTGCGCGCGTGGATTATc                                                          >  1:149666/1‑133 (MQ=255)                                                                 ggAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAgg                                                       >  1:323884/1‑133 (MQ=255)                                                                                  ggTGAGTCATGGTATGGCCCTGTGATGCCTTGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCg                                      <  1:160199/133‑1 (MQ=255)                                                                                             gTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGt                           >  1:470421/1‑133 (MQ=255)                                                                                                 tGCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGTCGtt                       >  2:291134/1‑133 (MQ=255)                                                                                                  gCACTGGGATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGTAAGCCTGGGGCTGGCGGCTGGATGGGTCGTTg                      >  2:209253/1‑133 (MQ=255)                                                                                                        ggATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTCGGTCGTTGAAACCg                >  2:113797/1‑133 (MQ=255)                                                                                                         gATGCCTGGTGAGTACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGTCGTTGAAACCGc               <  1:59467/133‑1 (MQ=255)                                                                                                                      tACGATTCTCGGACTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGTCGTTGAAACCGCAGGGGACATCTCg  >  2:7151/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                        GTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGGGATGCCTGGTGAGTACGATTCTCGGATTACCAGCATGGGCAGAGCGCCGCTTACGTCTGCGTAACTGTTCTATTTCGCGCGTGGATTATCAGGAAAGCCTGTGGCTGGCGTCTGGCTGGGTCGTTGAAACCGCAGGGGACATCTCG  >  NC_000913/4634156‑4634405

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