Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,571,850	T→C	T66A (ACC→GCC)	gadB ←	glutamate decarboxylase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,571,850	0	T	C	100.0%	86.3 / NA	27	T66A (ACC→GCC)	gadB	glutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/19);  total (8/19)

ATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGAT  >  NC_000913/1571612‑1572065                                                                                                                                                                                                                                               |                                                                                                                                                                                                                        aTACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAgg                                                                                                                                                                                                             >  2:177547/1‑250 (MQ=255)                 gCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACg                                                                                                                                                                                              >  2:316701/1‑250 (MQ=255)                           cGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCAt                                                                                                                                                                                   <  1:103410/251‑1 (MQ=255)                           cGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATGTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCAt                                                                                                                                                                                   <  1:332355/251‑1 (MQ=255)                                                ggCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTTACGAGCGTTGCCATCAAGATATAATTCGTCATTGa                                                                                                                                                              <  1:221214/251‑1 (MQ=255)                                                           aaCCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAAt                                                                                                                                                   <  1:298557/251‑1 (MQ=255)                                                                  ggTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACAt                                                                                                                                            >  2:37698/1‑251 (MQ=255)                                                                                               tCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACg                                                                                                                <  1:19868/250‑1 (MQ=255)                                                                                                     caggcgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGAt                                                                                                         >  2:209412/1‑251 (MQ=255)                                                                                                      aggcgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGAtt                                                                                                        >  1:310910/1‑251 (MQ=255)                                                                                                       ggcgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTc                                                                                                       <  1:178299/251‑1 (MQ=255)                                                                                                                                 atTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAc                                                                             <  2:293775/251‑1 (MQ=255)                                                                                                                                     acgcaacgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCaaa                                                                         <  1:316701/251‑1 (MQ=255)                                                                                                                                                      ggCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGagtagt                                                         <  2:198413/250‑1 (MQ=255)                                                                                                                                                      ggCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGtt                                                        <  1:37698/251‑1 (MQ=255)                                                                                                                                                            ggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGAcc                                                  <  1:344623/251‑1 (MQ=255)                                                                                                                                                            ggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGAcc                                                  >  2:247606/1‑251 (MQ=255)                                                                                                                                                                   ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATc                                           <  2:105931/251‑1 (MQ=255)                                                                                                                                                                                    gATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCGGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                               >  2:184877/1‑182 (MQ=255)                                                                                                                                                                                    gATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                               <  1:184877/182‑1 (MQ=255)                                                                                                                                                                                       ccAGTTTTTGGTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTACTTATCCa                       <  2:44645/251‑1 (MQ=255)                                                                                                                                                                                        cAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCAt                      <  1:81280/251‑1 (MQ=255)                                                                                                                                                                                             tttGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTa                  <  1:133478/250‑1 (MQ=255)                                                                                                                                                                                                gTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAAc               <  2:310910/250‑1 (MQ=255)                                                                                                                                                                                                  tAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGCAAGTGGACAGGTTCTGACGAGCATTGCAATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCTGACCTTAAATCTGTTACTTGCTTCTTATCCATTTTAAACTcc            >  2:191522/1‑251 (MQ=255)                                                                                                                                                                                                        aTAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTaaaa      <  2:66620/251‑1 (MQ=255)                                                                                                                                                                                                            aTCCATCAATTTGTGGACATTTTCGTCGTCCCAGGCCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGAt  <  2:133975/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                        ATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGAT  >  NC_000913/1571612‑1572065

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