Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,790,585	G→A	S220S (AGC→AGT)	ydiU ←	UPF0061 family protein YdiU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,790,585	0	G	A	100.0%	71.3 / NA	21	S220S (AGC→AGT)	ydiU	UPF0061 family protein YdiU
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/8);  total (13/8)

TATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTA  >  NC_000913/1790512‑1790658                                                                          |                                                                          tATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAAcc                                                                    <  2:307492/81‑1 (MQ=255)       tcaCCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGAc                                                              <  2:307493/81‑1 (MQ=255)        caCCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGACg                                                             >  2:307526/1‑81 (MQ=255)                gAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGt                                                     >  2:307533/1‑81 (MQ=255)                     aaGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCtca                                                >  2:307522/1‑81 (MQ=255)                             cGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTg                                         >  1:307525/1‑80 (MQ=255)                                   ccATTGGGCAATTAACGACGCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGAt                                   <  1:307494/80‑1 (MQ=255)                                                    acgCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGaa                 >  2:307531/1‑81 (MQ=255)                                                    acgCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGaa                 >  1:307535/1‑81 (MQ=255)                                                    acgCGGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGaa                 >  1:307547/1‑81 (MQ=255)                                                     cgcgGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGaa                 >  1:307527/1‑80 (MQ=255)                                                     cgcgGTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAAt                <  1:307495/81‑1 (MQ=255)                                                        ggTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAg              >  1:307536/1‑80 (MQ=255)                                                         gTACGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGc             >  1:307532/1‑80 (MQ=255)                                                            cGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAg         >  1:307537/1‑81 (MQ=255)                                                            cGTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAg         >  2:307534/1‑81 (MQ=255)                                                             gTGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGt        >  1:307539/1‑81 (MQ=255)                                                              tGCGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTc       <  2:307496/81‑1 (MQ=255)                                                                cGACAACATCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCa      <  2:307497/80‑1 (MQ=255)                                                                   caacaTCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTa  <  2:307499/81‑1 (MQ=255)                                                                   caacaTCACTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTa  <  1:307498/81‑1 (MQ=255)                                                                          |                                                                          TATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCAAGATGTGACCAGTAATGACGAATAGCAAAGTCAGCTA  >  NC_000913/1790512‑1790658

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