Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,376,226	T→C	L453P (CTA→CCA)	dcuD →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,376,226	0	T	C	100.0%	71.5 / NA	21	L453P (CTA→CCA)	dcuD	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/9);  total (12/9)

ACGCTAGCGCGCGCGGTTTCGCCGATTACTGCGGTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCT  >  NC_000913/3376081‑3376367                                                                                                                                                  |                                                                                                                                              acgCTAGCGCGCGCGGTTTCGCCGATTACTGCGGTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCAtttt                                                                                                                                          >  1:184989/1‑151 (MQ=255)                         aTTACTGCGGTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAAc                                                                                                                  <  2:119398/151‑1 (MQ=255)                                  gtggtggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGgccgc                                                                                                           >  1:209546/1‑149 (MQ=255)                                               cGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGt                                                                                            >  2:33358/1‑151 (MQ=255)                                                     gAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTa                                                                                        >  2:178668/1‑149 (MQ=255)                                                         tGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTcac                                                                                  <  2:287508/151‑1 (MQ=255)                                                                                 ggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGaccac                                                           >  1:173118/1‑150 (MQ=255)                                                                                 ggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGaccac                                                           >  2:47980/1‑150 (MQ=255)                                                                                 ggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGaccac                                                           >  2:289106/1‑150 (MQ=255)                                                                                               cGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCg                                            <  2:4468/151‑1 (MQ=255)                                                                                                   ccccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                        >  1:285374/1‑151 (MQ=255)                                                                                                   ccccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                        >  2:41559/1‑151 (MQ=255)                                                                                                   ccccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                        >  2:315315/1‑151 (MQ=255)                                                                                                         ggCAGTCGGTTTCGTGGCGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGTCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGGGGCTGCGGAGGTTCAt                                   >  2:66072/1‑150 (MQ=255)                                                                                                                  tttCGTGGTGAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATg                         <  2:117887/151‑1 (MQ=255)                                                                                                                       tggtgAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgtgt                    <  1:151601/151‑1 (MQ=255)                                                                                                                         gtgAATATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtca                  >  1:94032/1‑151 (MQ=255)                                                                                                                            aaTATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtcatca               <  1:210755/151‑1 (MQ=255)                                                                                                                            aaTATGATTGCCACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtcatc                <  1:205274/150‑1 (MQ=255)                                                                                                                                      ccACAATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCt      <  2:136868/150‑1 (MQ=255)                                                                                                                                         caATCACGCCATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCt  <  2:222873/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              ACGCTAGCGCGCGCGGTTTCGCCGATTACTGCGGTGGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCT  >  NC_000913/3376081‑3376367

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