Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779	3,374,780	G→A	L19L (CTG→CTA)	dcuD →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,374,780	0	G	A	100.0%	77.7 / NA	19	L19L (CTG→CTA)	dcuD	transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/6);  total (13/6)

TTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCC  >  NC_007779/3374683‑3374872                                                                                                  |                                                                                             tttATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTaaaa                                                                                           <  3:1008233/101‑1 (MQ=255)     tGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACt                                                                                       >  1:1091444/1‑101 (MQ=255)         aTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACaa                                                                                   >  3:1447772/1‑101 (MQ=255)              ccGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTc                                                                              >  1:908492/1‑101 (MQ=255)                     cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAggtggt                                                                       >  3:882156/1‑101 (MQ=255)                     cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAggtggt                                                                       <  3:569854/101‑1 (MQ=255)                                gAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCg                                                            >  3:1334620/1‑101 (MQ=255)                                      aacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTa                                                      >  1:166904/1‑101 (MQ=255)                                             ttCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCt                                               <  3:529635/101‑1 (MQ=255)                                               cGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtga                                             >  3:357022/1‑101 (MQ=255)                                                   aTAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTgatgat                                         >  1:539548/1‑101 (MQ=255)                                                    tAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGatgatg                                        <  1:1046377/101‑1 (MQ=255)                                                             cTGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                               >  3:499596/1‑101 (MQ=255)                                                              tGTCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCt                              <  3:1199193/101‑1 (MQ=255)                                                                tCATCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTgg                            >  1:1407710/1‑101 (MQ=255)                                                                   tCGTATTAATTACGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTa                         >  3:1471225/1‑101 (MQ=255)                                                                              aCGATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtgg              <  3:317193/101‑1 (MQ=255)                                                                                gATGGGCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGgtggtg            >  3:1354153/1‑101 (MQ=255)                                                                                     gCTATTTGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTc       <  1:873196/101‑1 (MQ=255)                                                                                          ttGATCCTAAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGGATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATcc  >  3:1278229/1‑101 (MQ=255)                                                                                                  |                                                                                             TTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCC  >  NC_007779/3374683‑3374872

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