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%	112.1 / NA	36	F43F (TTC→TTT)	gadB	glutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/20);  total (16/20)

TCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAA  >  NC_000913/1571783‑1572050                                                                                                                                       |                                                                                                                                      tcttTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATg                                                                                                                                   >  1:622508/1‑139 (MQ=255)     tGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGACCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGAc                                                                                                                               >  2:426814/1‑139 (MQ=255)      gTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACa                                                                                                                              >  1:1024451/1‑139 (MQ=255)      gTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACa                                                                                                                              <  2:1024451/139‑1 (MQ=255)      gTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACa                                                                                                                              <  1:1341537/139‑1 (MQ=255)             cAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTcgc                                                                                                                       <  2:622508/139‑1 (MQ=255)                               aaaTCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGaa                                                                                                     <  1:1307273/139‑1 (MQ=18)                                         aTTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGa                                                                                           <  2:633863/139‑1 (MQ=16)                                                  cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGACCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGAt                                                                                  <  1:426814/139‑1 (MQ=255)                                                                  ggTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                <  1:179547/125‑1 (MQ=255)                                                                  ggTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAg                                                                                >  2:179547/1‑125 (MQ=255)                                                                   gTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCg                                                                 >  2:1019761/1‑139 (MQ=255)                                                                        gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCa                                                            >  2:1248482/1‑139 (MQ=255)                                                                                  gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                                      >  2:654012/1‑103 (MQ=255)                                                                                  gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTg                                                                                      <  1:654012/103‑1 (MQ=255)                                                                                            gACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGAt                                                                                  >  2:930724/1‑97 (MQ=255)                                                                                            gACGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGAt                                                                                  <  1:930724/97‑1 (MQ=255)                                                                                             aCGAGCGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTcc                                       <  1:551050/139‑1 (MQ=255)                                                                                                  cGTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCt                                  >  1:999951/1‑139 (MQ=255)                                                                                                   gTTGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                                 <  2:999951/139‑1 (MQ=255)                                                                                                     tGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCg                                                                 <  1:1004938/105‑1 (MQ=255)                                                                                                     tGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCg                                                                 >  2:1004938/1‑105 (MQ=255)                                                                                                     tGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGtt                                         <  2:1337430/129‑1 (MQ=255)                                                                                                     tGCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGtt                                         >  1:1337430/1‑129 (MQ=255)                                                                                                      gCCATCAAGATATAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTaaa                              >  1:367207/1‑139 (MQ=255)                                                                                                               atatAATTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                     <  1:1199811/139‑1 (MQ=255)                                                                                                                    aTTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGaa                                                   >  1:674559/1‑104 (MQ=255)                                                                                                                    aTTCATCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGaa                                                   <  2:674559/104‑1 (MQ=255)                                                                                                                        atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGa                                                                                           >  1:601980/2‑60 (MQ=255)                                                                                                                        atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGa                                                                                           <  2:601980/59‑1 (MQ=255)                                                                                                                        atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                                 <  1:447546/117‑1 (MQ=255)                                                                                                                        atCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCtt                                 >  2:447546/2‑118 (MQ=255)                                                                                                                         tCATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt           <  1:1248482/139‑1 (MQ=255)                                                                                                                          cATTGATAATCTGAAATGCGACATCGTCGCGCATTTCGt                                                                                                              <  2:5154/39‑1 (MQ=255)                                                                                                                             tGATAATCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCAt       >  2:823875/1‑139 (MQ=255)                                                                                                                                  aTCTGAAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTaa  <  1:616338/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAA  >  NC_000913/1571783‑1572050

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