Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,571,902	G→A	D48D (GAC→GAT)	gadB ←	glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,571,902	0	G	A	100.0%	98.8 / NA	32	D48D (GAC→GAT)	gadB	glutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

TTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGC  >  NC_000913/1571779‑1572034                                                                                                                            |                                                                                                                                     ttcttcTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTga                                                                                                                       <  2:58032/139‑2 (MQ=255)  tcttctTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTgaa                                                                                                                      <  2:1355802/139‑3 (MQ=255)        ttGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGa                                                                                                                >  1:1416177/1‑139 (MQ=255)              aTCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACAtcgt                                                                                                          >  1:476808/1‑139 (MQ=255)                                          tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGttt                                                                              >  1:1390486/1‑139 (MQ=255)                                                  tgGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                      >  2:469229/1‑139 (MQ=255)                                                  tgGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                      <  1:469229/139‑1 (MQ=255)                                                      cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGAt                                                                  <  2:470131/139‑1 (MQ=255)                                                        ttttCGTCGTCCCAGGTCTGGTAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                <  2:265897/139‑1 (MQ=255)                                                        ttttCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                <  2:1416177/139‑1 (MQ=255)                                                         tttCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGt                                                               >  2:751627/1‑139 (MQ=21)                                                                      ggTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTc                                                  >  2:25133/1‑139 (MQ=255)                                                                           ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                             <  1:25133/139‑1 (MQ=255)                                                                           ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGAAcc                                             >  2:1119912/1‑139 (MQ=255)                                                                            gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCa                                            >  1:1530452/1‑139 (MQ=255)                                                                                      gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAAt                                  >  2:77469/1‑139 (MQ=255)                                                                                      gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAAt                                  >  2:852079/1‑139 (MQ=255)                                                                                          ggTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGag                              <  1:865748/139‑1 (MQ=255)                                                                                            ttCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGagta                            >  1:927682/1‑139 (MQ=255)                                                                                                gACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTc                        >  2:1441810/1‑139 (MQ=255)                                                                                                     gCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGAcc                   >  1:477777/1‑139 (MQ=255)                                                                                                      cGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAAt                                                         >  1:929728/1‑100 (MQ=255)                                                                                                      cGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCt                  >  2:1142549/1‑139 (MQ=255)                                                                                                      cGTTGCCATCAAGATATAATGCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAAt                                                         <  2:929728/100‑1 (MQ=255)                                                                                                       gTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                 <  2:477777/139‑1 (MQ=255)                                                                                                       gTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                 <  1:852079/139‑1 (MQ=255)                                                                                                              tCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCg          <  1:1142549/139‑1 (MQ=255)                                                                                                                aaGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTg                                                              <  1:992746/85‑1 (MQ=255)                                                                                                                aaGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTg                                                              >  2:992746/1‑85 (MQ=255)                                                                                                                aaGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                 <  2:446072/130‑1 (MQ=255)                                                                                                                aaGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                 >  1:446072/1‑130 (MQ=255)                                                                                                                      taATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGc  <  1:1441810/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                     TTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGC  >  NC_000913/1571779‑1572034

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