Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,787	G→C	A308A (GCG→GCC)	cycA →	serine/alanine/glycine/cycloserine:H(+)symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,430,787	0	G	C	100.0%	166.3 / NA	48	A308A (GCG→GCC)	cycA	serine/alanine/glycine/cycloserine:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (18/30);  total (18/30)

GTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCGAAAG  >  NC_000913/4430709‑4430863                                                                               |                                                                             gTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTc                                                                            <  1:915248/81‑1 (MQ=255) gTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTc                                                                            >  2:915324/1‑81 (MQ=255) gTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTc                                                                            <  1:915246/81‑1 (MQ=255) gTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTc                                                                            <  2:915247/81‑1 (MQ=255)  tttGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCctct                                                                           >  1:915318/1‑81 (MQ=255)    tgttgAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCtt                                                                          <  2:915249/80‑1 (MQ=255)     gttgAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTc                                                                         <  2:915250/80‑1 (MQ=255)        gAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGc                                                                      >  2:915322/1‑80 (MQ=255)        gAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGc                                                                      >  2:915321/1‑80 (MQ=255)        gAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGc                                                                      >  1:915311/1‑80 (MQ=255)         aaCTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTa                                                                    <  2:915251/81‑1 (MQ=255)           cTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAAc                                                                  <  2:915252/81‑1 (MQ=255)            tGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTaa                                                                   >  2:915279/1‑79 (MQ=255)              ttCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAgc                                                               <  1:915253/81‑1 (MQ=255)               tCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAgcg                                                              >  1:915327/1‑81 (MQ=255)                  tgtTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAgcggcg                                                           <  1:915254/81‑1 (MQ=255)                       gTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGtctt                                                       <  1:915257/80‑1 (MQ=255)                       gTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGtctt                                                       <  2:915256/80‑1 (MQ=255)                       gTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGtctt                                                       <  2:915255/80‑1 (MQ=255)                        tAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGtcttct                                                     <  1:915259/81‑1 (MQ=255)                        tAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGtcttct                                                     <  2:915258/81‑1 (MQ=255)                          gggCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTtctc                                                    <  1:915260/80‑1 (MQ=255)                          gggCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTtctc                                                    <  1:915263/80‑1 (MQ=255)                          gggCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTctct                                                   <  2:915261/81‑1 (MQ=255)                          gggCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTctct                                                   <  2:915262/81‑1 (MQ=255)                                cctgcTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCa                                               <  2:915264/79‑1 (MQ=255)                                 ctgctgCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGcc                                            <  1:915265/81‑1 (MQ=255)                                    ctgcCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCg                                           <  2:915266/79‑1 (MQ=255)                                    ctgcCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTa                                         >  1:915331/1‑81 (MQ=255)                                    ctgcCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTa                                         >  2:915317/1‑81 (MQ=255)                                         gCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTg                                    >  2:915343/1‑81 (MQ=255)                                          cAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGt                                   >  2:915332/1‑81 (MQ=255)                                           aaGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGtt                                  >  1:915314/1‑81 (MQ=255)                                           aaGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGtt                                  <  1:915267/81‑1 (MQ=255)                                            aGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGttt                                 <  1:915268/81‑1 (MQ=255)                                            aGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGttt                                 >  1:915269/1‑81 (MQ=255)                                             gCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTg                                >  2:915341/1‑81 (MQ=255)                                                tGATCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTc                             >  1:915330/1‑81 (MQ=255)                                                  aTCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCt                            <  2:915271/80‑1 (MQ=255)                                                  aTCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTg                           >  2:915333/1‑81 (MQ=255)                                                  aTCAACTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTg                           <  2:915270/81‑1 (MQ=255)                                                     aaCTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGgcg                        >  1:915339/1‑81 (MQ=255)                                                     aaCTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGgc                         <  1:915272/80‑1 (MQ=255)                                                      aCTTTGTGGTGCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGgcgc                       <  2:915273/81‑1 (MQ=255)                                                               tgCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGgtg              <  1:915274/81‑1 (MQ=255)                                                                gCTGACCTCTGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGgtgt             >  1:915320/1‑81 (MQ=255)                                                                      ctctGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCAcc       <  2:915275/81‑1 (MQ=255)                                                                        ctGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCg      >  2:915328/1‑80 (MQ=255)                                                                         tGCGGCCTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCg      <  2:915276/79‑1 (MQ=255)                                                                           cggccTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCGaa    <  1:915277/79‑1 (MQ=255)                                                                           cggccTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCGAAAg  >  2:915336/1‑81 (MQ=255)                                                                               |                                                                             GTTTGTTGAACTGTTCGTGTTGGTAGGGCTGCCTGCTGCCGCAAGCGTGATCAACTTTGTGGTGCTGACCTCTGCGGCGTCTTCCGCTAACAGCGGCGTCTTCTCTACCAGCCGTATGCTGTTTGGTCTGGCGCAGGAAGGTGTGGCACCGAAAG  >  NC_000913/4430709‑4430863

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