Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,570,763	C→T	G428D (GGC→GAC)	gadB ←	glutamate decarboxylase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,570,763	0	C	T	100.0%	85.0 / NA	37	G428D (GGC→GAC)	gadB	glutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/15);  total (22/15)

CCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTT  >  NC_000913/1570626‑1570899                                                                                                                                          |                                                                                                                                         ccgttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTcg                                                                                                                                        >  2:55487/1‑140 (MQ=255)    ttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGcga                                                                                                                                     >  2:890970/1‑140 (MQ=255)    ttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGcga                                                                                                                                     >  2:1097867/1‑140 (MQ=255)       ccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGcgacga                                                                                                                                  >  1:251032/1‑140 (MQ=255)               cGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAAt                                                                                                                          >  1:155326/1‑140 (MQ=37)               cGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAAt                                                                                                                          <  1:484627/140‑1 (MQ=25)                  tATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATgcg                                                                                                                       >  2:990080/1‑140 (MQ=25)                     cAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCAt                                                                                                                    <  2:155326/140‑1 (MQ=21)                           tGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATcaccac                                                                                                              >  2:423212/1‑140 (MQ=18)                           tGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATcaccac                                                                                                              >  1:854089/1‑140 (MQ=17)                            gTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACg                                                                                                             <  2:141776/140‑1 (MQ=17)                               tAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATg                                                                                                          <  1:990080/140‑1 (MQ=17)                                  aGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCg                                                                                                       >  2:561384/1‑140 (MQ=18)                                   gCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCgg                                                                                                      <  2:854089/140‑1 (MQ=17)                                                                                                       cTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                  >  1:622490/1‑140 (MQ=11)                                                                                                         tCCAGCAACAGTTCAGCGAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                >  1:805527/1‑140 (MQ=12)                                                                                                          ccAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATAc                               >  1:504301/1‑140 (MQ=17)                                                                                                           cAGCAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACa                              >  1:1218923/1‑140 (MQ=17)                                                                                                              cAACAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                           <  1:423212/140‑1 (MQ=18)                                                                                                                 cAGTTCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                        <  2:633416/140‑1 (MQ=18)                                                                                                                    ttCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGGGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATc                     >  1:1069915/1‑140 (MQ=12)                                                                                                                    ttCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATAc                               >  2:47131/1‑130 (MQ=17)                                                                                                                    ttCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATAc                               <  1:47131/130‑1 (MQ=18)                                                                                                                    ttCAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATc                     >  2:167550/1‑140 (MQ=17)                                                                                                                      cAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                   <  2:1112256/140‑1 (MQ=18)                                                                                                                      cAGCAAAGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                   <  2:804333/140‑1 (MQ=18)                                                                                                                            aGTCCATTTCGAAGTCGCGACGACACCTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                   <  2:1111412/134‑1 (MQ=18)                                                                                                                            aGTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                   >  1:1111412/1‑134 (MQ=17)                                                                                                                             gTCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCtt            <  2:504301/140‑1 (MQ=18)                                                                                                                              tCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc           >  2:353570/1‑140 (MQ=17)                                                                                                                              tCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc           >  2:453889/1‑140 (MQ=17)                                                                                                                              tCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc           >  2:541151/1‑140 (MQ=17)                                                                                                                              tCCATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGGGTATCCCGGATCTTc           >  2:844555/1‑140 (MQ=17)                                                                                                                               ccATTTCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa          <  2:1059123/140‑1 (MQ=18)                                                                                                                                 tttacaaaGCCGGGACGACACATAATGCGCATCACCACGTTGTCGGTGGTTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAcc        <  2:622490/134‑1 (MQ=11)                                                                                                                                  tttCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa       <  1:1096121/140‑1 (MQ=18)                                                                                                                                  tttCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa       <  2:974479/140‑1 (MQ=18)                                                                                                                                   ttCGAAGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGTCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCAt      >  1:319582/1‑140 (MQ=17)                                                                                                                                       aaGTCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCttt  <  2:596173/140‑1 (MQ=18)                                                                                                                                          |                                                                                                                                         CCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTT  >  NC_000913/1570626‑1570899

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