Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,571,917	G→A	F43F (TTC→TTT)	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,917	0	G	A	100.0%	144.9 / NA	48	F43F (TTC→TTT)	gadB	glutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/22);  total (26/22)

TTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAA  >  NC_000913/1571786‑1572051                                                                                                                                    |                                                                                                                                       ttGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGa                                                                                                                                 >  2:471714/1‑139 (MQ=255)                  gTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCg                                                                                                                >  1:647426/1‑139 (MQ=255)                       tGGGTAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGc                                                                                                                                   >  2:965333/1‑115 (MQ=9)                       tGGGTAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGc                                                                                                                                   <  1:965333/115‑1 (MQ=2)                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACAtcgtcg                                                                                                                         <  2:356449/113‑1 (MQ=12)                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACAtcgtcg                                                                                                                         >  1:356449/1‑113 (MQ=12)                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGttt                                                                                               >  2:661000/1‑139 (MQ=25)                                   tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGttt                                                                                               >  1:615361/1‑139 (MQ=255)                                    cAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGtttt                                                                                              >  2:1641309/1‑139 (MQ=255)                                    cAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGtttt                                                                                              >  1:789416/1‑139 (MQ=255)                                        ttGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGAtt                                                                                          >  1:1277008/1‑139 (MQ=255)                                                 ttttCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 <  2:1046040/139‑1 (MQ=255)                                                 ttttCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 <  2:789416/139‑1 (MQ=255)                                                                    ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTc                                                                   >  2:562954/1‑134 (MQ=255)                                                                    ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTc                                                                   <  1:562954/134‑1 (MQ=255)                                                                    ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                              <  2:828836/139‑1 (MQ=255)                                                                    ggCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                              >  2:680946/1‑139 (MQ=255)                                                                     gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 >  2:957442/1‑119 (MQ=255)                                                                     gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 >  2:851530/1‑119 (MQ=255)                                                                     gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 <  1:851530/119‑1 (MQ=255)                                                                     gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 <  1:957442/119‑1 (MQ=255)                                                                         aaaGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAAc                                                         <  2:615361/139‑1 (MQ=255)                                                                           aGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGt                                                       <  2:647426/139‑1 (MQ=255)                                                                                         gACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTc                                         >  1:1021343/1‑139 (MQ=255)                                                                                                    ccATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAAt                              <  1:680946/139‑1 (MQ=255)                                                                                                            atatAATTCATCCTTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTTAATCTAGTAGTTCCGACCTTAAATCCGTTACt                      <  2:1031408/139‑1 (MQ=18)                                                                                                              ataATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTg                    >  2:1117490/1‑139 (MQ=255)                                                                                                                 aTTCATCATTGATAGTCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                      <  1:1415034/134‑1 (MQ=255)                                                                                                                 aTTCATCATTGATAGTCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                      >  2:1415034/1‑134 (MQ=255)                                                                                                                   tcatCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                                       >  2:1229122/4‑67 (MQ=14)                                                                                                                   tcatCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                                       <  1:1229122/64‑1 (MQ=255)                                                                                                                     atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                              <  1:52455/89‑1 (MQ=255)                                                                                                                     atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCAcc                                                              >  2:52455/2‑90 (MQ=255)                                                                                                                     atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTGCGACCTTAAATCCGTTACTTGCTTCTTa             >  1:1040791/2‑139 (MQ=255)                                                                                                                     atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGt                          >  2:334751/2‑126 (MQ=255)                                                                                                                     atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGt                          <  1:334751/125‑1 (MQ=255)                                                                                                                      tCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt            >  1:1266974/1‑139 (MQ=255)                                                                                                                      tCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt            <  1:471714/139‑1 (MQ=255)                                                                                                                       cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCg                                                                                     >  2:649330/1‑65 (MQ=255)                                                                                                                       cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCg                                                                                     <  1:649330/65‑1 (MQ=255)                                                                                                                       cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 <  1:919941/69‑1 (MQ=255)                                                                                                                       cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                 >  2:919941/1‑69 (MQ=255)                                                                                                                       cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATc           >  2:856003/1‑139 (MQ=255)                                                                                                                        aTTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATcc          <  1:1490136/139‑1 (MQ=255)                                                                                                                         ttGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGt                                                       <  2:1698911/93‑1 (MQ=255)                                                                                                                         ttGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGt                                                       >  1:1698911/1‑93 (MQ=255)                                                                                                                                 cTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTaaa  >  1:1527376/1‑138 (MQ=255)                                                                                                                                 cTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTaaa  <  2:1527376/138‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       TTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAA  >  NC_000913/1571786‑1572051

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