Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,272,344	Δ1 bp	coding (408/885 nt)	garR ←	tartronate semialdehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,272,342	0	C	.	100.0%	146.1 / NA	35	coding (410/885 nt)	garR	tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (23/12);  total (23/12)

CGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAAC  >  NC_000913/3272218‑3272471                                                                                                                             |                                                                                                                                  cgcCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCa                                                                                                                          >  2:87334/1‑133 (MQ=255)    cacaATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTg                                                                                                                       >  2:40520/1‑133 (MQ=255)     acaATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGa                                                                                                                      <  2:268819/133‑1 (MQ=255)     acaATGACCTCATTTTCCAGTTTGCTGACGTTCCCTGCACCCATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGa                                                                                                                      <  2:168608/133‑1 (MQ=255)      caATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGAc                                                                                                                     >  1:475766/1‑133 (MQ=255)                aTTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGt                                                                                                           >  2:405570/1‑133 (MQ=255)                aTTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGt                                                                                                           >  1:407783/1‑133 (MQ=255)                        gTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCt                                                                                                   <  2:364587/133‑1 (MQ=255)                         tttGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCGCCATCACTGACAGCATACCGTCGATGGCtt                                                                                                  >  1:202815/1‑133 (MQ=255)                          ttGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CTACCATCATTGACAGCGCACCGTCGGTGGCttt                                                                                                 >  2:208900/1‑133 (MQ=255)                          ttGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCttt                                                                                                 >  1:339756/1‑133 (MQ=255)                          ttGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCttt                                                                                                 <  1:375848/133‑1 (MQ=255)                                         ccACCGATTCACCCGTTATGCACCCCGGAACCCGCCATCGCATTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCt                                                                                  <  2:202815/132‑1 (MQ=255)                                               aTTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGAGAGCGTACCGTCGGTGGTTTTCGGTTAACCGACGCTCACCgg                                                                            >  2:135254/1‑133 (MQ=255)                                                tttCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGGCGAGGGCGTTCGGTTCACCGCGGCTCACCGGa                                                                           >  2:151410/1‑133 (MQ=255)                                                  tCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGc                                                                         <  2:475766/133‑1 (MQ=255)                                                   ccccGGTACGCACCAGGGTACCCGCCATCGTTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTAGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCa                                                                        <  1:221626/133‑1 (MQ=255)                                                         tATGCCCCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGc                                                                  <  2:301566/133‑1 (MQ=255)                                                           tGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCat                                                                >  2:458831/1‑133 (MQ=255)                                                                  cGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATg                                                         >  1:506763/1‑133 (MQ=255)                                                                  cGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATg                                                         >  2:431608/1‑133 (MQ=255)                                                                    gAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGcc                                                       >  2:302125/1‑133 (MQ=255)                                                                               cGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCa                                            >  1:20806/1‑133 (MQ=255)                                                                                                     ttGTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGTCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCt                      >  1:52559/1‑133 (MQ=255)                                                                                                       gTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTg                    >  1:375466/1‑133 (MQ=255)                                                                                                       gTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTg                    >  2:80703/1‑133 (MQ=255)                                                                                                       gTCGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTACGCGTTCGGCGCTTCGTTGATTTCACGGCTTg                    >  2:169357/1‑133 (MQ=255)                                                                                                         cGAAAATAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGcc                  >  2:370306/1‑133 (MQ=255)                                                                                                           aaaaTAGCTTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAg                <  1:87334/133‑1 (MQ=255)                                                                                                           aaaaTAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAg                <  2:59918/133‑1 (MQ=255)                                                                                                              aTAGCCTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCgg             <  2:483963/133‑1 (MQ=255)                                                                                                                aGCCTTGTCGCCG‑CCTCCCTCATTGACAGCCTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTg           <  1:151410/133‑1 (MQ=255)                                                                                                                  ccTTGTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGACTTCACGACTTGCCAGCGGTGCg         >  1:228643/1‑133 (MQ=255)                                                                                                                      gTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAg     >  1:299590/1‑133 (MQ=255)                                                                                                                      gTCGCCG‑CCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAg     >  1:292909/1‑133 (MQ=255)                                                                                                                         gccgccC‑CCTTCCCTGCCAGCGTACCGTCGATGGTTTTGGTTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAAc  <  1:40520/133‑1 (MQ=255)                                                                                                                             |                                                                                                                                  CGCCACAATGACCTGATTTGCCAGTTTGGTGACGTTACCTGCACCGATTTCCCCGGTATGCACCACGGAACCCGCCATCGCTTTCATCAAATCATAGTATTTGTCGAAAATAGCCTTGTCGCCGCCCACCATCACTGACAGCGTACCGTCGATGGCTTTCGGTTCACCGCCGCTCACCGGAGCATCCAGCATATCAATGCCTTTCGCTTTCAGCGCTTCGCTGATTTCACGGCTTGCCAGCGGTGCGATAGAAC  >  NC_000913/3272218‑3272471

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