Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,649,867	T→G	V25G (GTT→GGT)	ydfD →	Qin prophage; lysis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,649,867	0	T	G	100.0%	94.5 / NA	26	V25G (GTT→GGT)	ydfD	Qin prophage; lysis protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/17);  total (9/17)

AAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAA  >  NC_000913/1649736‑1650002                                                                                                                                    |                                                                                                                                        aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  2:1591549/140‑1 (MQ=255) aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  1:1419067/140‑1 (MQ=255) aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  1:1440893/140‑1 (MQ=255) aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  1:1623756/140‑1 (MQ=255) aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  1:1446014/140‑1 (MQ=255) aaataTGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGa                                                                                                                                 <  1:1521107/140‑1 (MQ=255)           tGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGa                                                                                                                       <  1:1608995/140‑1 (MQ=255)             tttCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAgt                                                                                                                     >  1:1477600/1‑140 (MQ=255)                cAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTa                                                                                                                  <  1:1702784/140‑1 (MQ=255)                    gCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGAc                                                                                                              <  1:1427045/140‑1 (MQ=255)                     cTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACt                                                                                                             >  2:1594091/1‑140 (MQ=255)                       ggCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTgc                                                                                                           >  1:1334770/1‑140 (MQ=255)                         cTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTgcag                                                                                                         <  2:1349744/140‑1 (MQ=255)                                        tGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGaaaa                                                                                          >  2:1496339/1‑140 (MQ=255)                                                ggATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGt                                                                                  <  1:1442481/140‑1 (MQ=255)                                                            tGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCg                                                                      <  2:1514198/140‑1 (MQ=255)                                                                 tCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGa                                                                 <  2:1489886/140‑1 (MQ=255)                                                                      aTTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCGGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAg                                                            >  2:1373286/1‑140 (MQ=255)                                                                        ttGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGtt                                                          >  1:1320066/1‑140 (MQ=255)                                                                                     gACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGata                                             <  1:1484996/140‑1 (MQ=255)                                                                                           ttttCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCg                                       >  1:1502469/1‑140 (MQ=255)                                                                                                       cGGAGAGCCAGTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAg                           <  2:1334770/140‑1 (MQ=255)                                                                                                                 gTTGATATTGCAGTCAGTGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTaaaa                 <  2:1643664/140‑1 (MQ=255)                                                                                                                      tatTGCAGTCAGTGGTCACAGGACAATGCGGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGtt            <  1:1373286/140‑1 (MQ=255)                                                                                                                            agtcagtGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGtaa      >  1:1439171/1‑140 (MQ=255)                                                                                                                               cagtGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATACCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGtaataa   >  2:1656991/1‑140 (MQ=255)                                                                                                                                agtGGTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTGCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTaataaa  >  2:1638938/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                        AAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCAGTGTTCACAGGACAATGCAGGAGTGTATGACTGCAGCAACCGAACAGAAAATTCCCGGTAACTGTTACCCGGTCGATAAAGTTATTCACCAGGATAATATCGAAATCCCGGCAGGTCTTTAAAACAGTTCCGTAATAAA  >  NC_000913/1649736‑1650002

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