Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,648,553	G→A	V16M (GTG→ATG)	ydfA →	Qin prophage; uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,648,553	0	G	A	100.0%	116.2 / NA	30	V16M (GTG→ATG)	ydfA	Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/13);  total (17/13)

ACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCT  >  NC_000913/1648405‑1648694                                                                                                                                                     |                                                                                                                                              aCCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCatg                                                                                                                                             <  2:488378/151‑1 (MQ=255)       atGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCa                                                                                                                                       <  1:170066/151‑1 (MQ=255)           tAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAAc                                                                                                                                    <  1:19484/150‑1 (MQ=255)             gTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAg                                                                                                                                 <  2:48930/151‑1 (MQ=255)                    cTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGc                                                                                                                           >  1:488945/1‑150 (MQ=255)                      aaTCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCAc                                                                                                                         <  2:328566/150‑1 (MQ=255)                       aTCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGt                                                                                                                       >  2:309725/1‑151 (MQ=255)                                gCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                                                                               <  2:547206/150‑1 (MQ=255)                                     aaCACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                                                         >  1:576228/1‑151 (MQ=255)                                     aaCACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                                                         >  2:301924/1‑151 (MQ=255)                                     aaCACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                                                                         >  2:603296/1‑151 (MQ=255)                                               aCGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaaa                                                                                               <  2:520637/151‑1 (MQ=255)                                                cGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                                                                                >  1:507735/1‑149 (MQ=255)                                                  gtgtTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                                                                                             >  2:318857/1‑150 (MQ=255)                                                     ttCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                                                                          >  1:346331/1‑150 (MQ=255)                                                                cAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAcc                                                                              >  1:359661/1‑151 (MQ=255)                                                                     ccGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaa                                                                          <  2:411145/150‑1 (MQ=255)                                                                       gCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATg                                                                       >  1:440489/1‑151 (MQ=255)                                                                        cTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                                                                      >  2:614855/1‑151 (MQ=255)                                                                              cAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                                                                >  2:322234/1‑151 (MQ=255)                                                                              cAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                                                                >  2:122246/1‑151 (MQ=255)                                                                               aGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGt                                                               <  1:171794/151‑1 (MQ=255)                                                                                             aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaag                                                 <  1:330952/151‑1 (MQ=255)                                                                                               gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                >  2:95747/1‑150 (MQ=255)                                                                                                     aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAg                                         <  2:46673/151‑1 (MQ=255)                                                                                                        aTGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtg                                      >  2:390422/1‑151 (MQ=255)                                                                                                                tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                              >  2:393059/1‑151 (MQ=255)                                                                                                                        ttGGCAACAACGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCg                      >  2:596875/1‑151 (MQ=255)                                                                                                                            caacaaCGAATCTCTGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac                  <  1:28964/151‑1 (MQ=255)                                                                                                                                         ctGGTGTACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTAAACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGaa        <  2:151820/148‑1 (MQ=255)                                                                                                                                            gtgtACGGCATGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCt  <  2:427585/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              ACCGATATGTTAGTCATGGCTAATCTTGTTTGCATCAACACACCGCACGGTGTTCTCAGCAAACAGTTCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCT  >  NC_000913/1648405‑1648694

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