Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,570,678	A→G	G456G (GGT→GGC)	gadB ←	glutamate decarboxylase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,570,678	0	A	G	100.0%	153.5 / NA	47	G456G (GGT→GGC)	gadB	glutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/23);  total (24/23)

TATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAC  >  NC_000913/1570546‑1570808                                                                                                                                     |                                                                                                                                   tATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCa                                                                                                                             >  2:1020182/1‑140 (MQ=255) gaTTCAGAACAATATTTTCCGTTGCTAGTGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGCTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCa                                                                                                                             >  2:452707/2‑140 (MQ=255)   ttCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGt                                                                                                                           <  1:790967/140‑1 (MQ=255)      aGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTc                                                                                                                        >  2:767032/1‑140 (MQ=255)        aaCAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCAGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCgg                                                                                                                      >  2:754315/1‑140 (MQ=255)          cAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGt                                                                                                                    <  1:8752/140‑1 (MQ=255)           aaTATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTg                                                                                                                   >  2:556745/1‑140 (MQ=255)                tttCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGc                                                                                                              >  2:6165/1‑140 (MQ=255)                    cGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTgag                                                                                                          <  1:360165/140‑1 (MQ=255)                     gTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTgaga                                                                                                         <  2:504686/140‑1 (MQ=255)                                cAGTGAACAGACTTTGGAAATTGTCCCGAAACGTGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGa                                                                                              <  1:63751/140‑1 (MQ=255)                                cAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGa                                                                                              <  2:752465/140‑1 (MQ=255)                                 aGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAg                                                                                             <  2:392441/140‑1 (MQ=255)                                   tGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGc                                                                                           >  2:335420/1‑140 (MQ=255)                                   tGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGc                                                                                           >  1:356657/1‑140 (MQ=255)                                              tGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCtt                                                                                >  1:577646/1‑140 (MQ=255)                                               ggAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTc                                                                               >  1:747978/1‑140 (MQ=255)                                                     tGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCaa                                                                         <  2:572400/140‑1 (MQ=255)                                                     tGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCaa                                                                         <  1:543902/140‑1 (MQ=255)                                                                 ggTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCaaa                                                             >  2:484008/1‑140 (MQ=255)                                                                 ggTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCaaa                                                             >  2:401138/1‑140 (MQ=255)                                                                 ggTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCa                                                               <  2:947753/138‑1 (MQ=255)                                                                 ggTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCa                                                               >  1:947753/1‑138 (MQ=255)                                                                                 ccgttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcg                                             <  1:49190/140‑1 (MQ=255)                                                                                 ccgttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcg                                             <  1:174940/140‑1 (MQ=255)                                                                                  cgttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcgc                                            >  1:753449/1‑140 (MQ=255)                                                                                  cgttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcgc                                            >  1:976122/1‑140 (MQ=255)                                                                                   gttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAGAGTCCATTTCGAAGCcgcg                                           <  1:754315/140‑1 (MQ=255)                                                                                   gttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcgcg                                           <  1:335420/140‑1 (MQ=255)                                                                                   gttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCcgcg                                           <  1:6165/140‑1 (MQ=255)                                                                                    ttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGcga                                          <  2:577646/140‑1 (MQ=255)                                                                                    ttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGcga                                          >  1:595273/1‑140 (MQ=255)                                                                                     taccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGcgac                                         <  1:398802/140‑1 (MQ=255)                                                                                     taccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGcgac                                         <  1:556745/140‑1 (MQ=255)                                                                                            aaaCGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACAt                                  >  1:469059/1‑140 (MQ=255)                                                                                                  tATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATgcg                            <  2:976122/140‑1 (MQ=21)                                                                                                    tCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCa                          <  1:484008/140‑1 (MQ=21)                                                                                                     cAGGTATGTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCAt                         <  2:595273/140‑1 (MQ=21)                                                                                                            gTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACg                  >  2:397734/1‑140 (MQ=18)                                                                                                            gTTTAAAGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACg                  >  1:359927/1‑140 (MQ=17)                                                                                                                 aaGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTc             >  1:950470/1‑140 (MQ=18)                                                                                                                  aGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCg            >  1:628636/1‑140 (MQ=18)                                                                                                                  aGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCg            >  1:737515/1‑140 (MQ=18)                                                                                                                  aGCTGTTCTGTTGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCg            >  2:179390/1‑140 (MQ=17)                                                                                                                    ctgttctgttGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGt          <  1:397734/140‑1 (MQ=17)                                                                                                                       ttctgttGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGc       >  1:78726/1‑140 (MQ=18)                                                                                                                       ttctgttGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGc       >  2:299162/1‑140 (MQ=18)                                                                                                                            ttGGGCAATGCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAc  <  2:137930/140‑1 (MQ=11)                                                                                                                                     |                                                                                                                                   TATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAC  >  NC_000913/1570546‑1570808

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