Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,853,793	G→T	G16G (GGC→GGA)	ysdE ←	protein YsdE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,853,793	0	G	T	100.0%	73.9 / NA	23	G16G (GGC→GGA)	ysdE	protein YsdE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (8/15);  total (8/15)

GCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATTCCAGCTATTGCCGATATGGCGCGCGATC  >  NC_000913/3853549‑3854033                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTa                                                                                                                                                                                                                                            <  1:41163/251‑1 (MQ=255)       gAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAAt                                                                                                                                                                                                                                      >  1:212326/1‑251 (MQ=255)                    atCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTgg                                                                                                                                                                                                                         <  2:188411/251‑1 (MQ=255)                          aTTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGtat                                                                                                                                                                                                                    >  2:34630/1‑250 (MQ=255)                                aTCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGc                                                                                                                                                                                                             <  1:104249/251‑1 (MQ=255)                                           cATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATtcatc                                                                                                                                                                                                  <  1:73393/251‑1 (MQ=255)                                                                               gAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAAt                                                                                                                                                               <  2:113542/250‑1 (MQ=255)                                                                                aaaTACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTGAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCTTATATCTGGCTAACAGTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGc                                                                                                                                                             <  1:156633/251‑1 (MQ=255)                                                                                                   gATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAatt                                                                                                                                          <  2:100655/251‑1 (MQ=255)                                                                                                    aTTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAatta                                                                                                                                         <  2:19193/251‑1 (MQ=255)                                                                                                       cAAGTCGCACCAAAGGGTAGAGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTAttt                                                                                                                                      >  2:77278/1‑251 (MQ=255)                                                                                                                        gAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTGTTATAAAAGCATCCGTg                                                                                                                     <  2:228800/251‑1 (MQ=255)                                                                                                                                    gCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAgg                                                                                                         >  2:223328/1‑251 (MQ=255)                                                                                                                                       ccccTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCaa                                                                                                      <  1:55828/251‑1 (MQ=255)                                                                                                                                         ccTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCaa                                                                                                      >  1:43571/1‑249 (MQ=255)                                                                                                                                          cgtttATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTGAGAAAAGCATCCGTGATAATGAAAAGGCAAAGa                                                                                                   <  1:212414/249‑1 (MQ=255)                                                                                                                                                            gggTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTACTTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                                                 >  1:15931/1‑251 (MQ=255)                                                                                                                                                                                                             ctcCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt                                    <  2:72045/247‑1 (MQ=255)                                                                                                                                                                                                             ctcCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATACGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta                                  <  2:165168/249‑1 (MQ=255)                                                                                                                                                                                                               cccGGCAACATTTATTGCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtat                                 <  1:181660/248‑1 (MQ=255)                                                                                                                                                                                                                               gCCGCTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATTCCAGCTATTGCCGat               >  1:46147/1‑250 (MQ=255)                                                                                                                                                                                                                                   cTTTTGTTTACATATTCTTCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGCTGTTGGGATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATTCCAGCTATTGCCGATATGgc          >  1:12920/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCATATTCTTCCGCGAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATTCCAGCTATTGCCGATATGGCGCGCGAtc  <  1:188017/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATTCCAGCTATTGCCGATATGGCGCGCGATC  >  NC_000913/3853549‑3854033

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