Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,571,332 G→T I238I (ATC→ATA gadB ← glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,571,3320GT100.0% 47.8 / NA 21I238I (ATC→ATAgadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/7);  total (14/7)

TGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATC  >  NC_000913/1571227‑1571464
                                                                                                         |                                                                                                                                    
tGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATcc                                                                                                     >  1:330590/1‑139 (MQ=18)
             gATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCa                                                                                        >  1:1382846/1‑139 (MQ=18)
                    cacGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTAt                                                                                                                      <  1:136905/102‑1 (MQ=14)
                    cacGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTAt                                                                                                                      >  2:136905/1‑102 (MQ=14)
                         ggCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGa                                                                            >  1:1534199/1‑139 (MQ=17)
                           cAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAAc                                                                          <  1:702131/139‑1 (MQ=18)
                                   aGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt                                                                                   >  1:1451902/1‑122 (MQ=17)
                                   aGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt                                                                                   <  2:1451902/122‑1 (MQ=17)
                                   aGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGtt                                                                  >  2:1506880/1‑139 (MQ=18)
                                                         gacgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                >  1:680809/1‑119 (MQ=18)
                                                         gacgaACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                <  2:680809/119‑1 (MQ=18)
                                                                  tGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGt                                                                                   >  2:937859/1‑91 (MQ=18)
                                                                  tGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAAATTATCCAGCGCATCGTGCAGCGGt                                                                                   <  1:937859/91‑1 (MQ=255)
                                                                            ccACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCGAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCAAGCCGATgg                         >  1:910618/1‑139 (MQ=255)
                                                                            ccACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                >  1:1336784/1‑100 (MQ=18)
                                                                            ccACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                <  2:1336784/100‑1 (MQ=18)
                                                                            ccACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATgg                         >  1:1415670/1‑139 (MQ=25)
                                                                            ccACCGCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATgg                         >  1:408978/1‑139 (MQ=37)
                                                                                 gCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                <  1:1379253/95‑1 (MQ=18)
                                                                                 gCTGGCAGCGTCGATGTGCATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTAc                                                                >  2:1379253/1‑95 (MQ=18)
                                                                                       aGCGTCGATGTGGATGTCTATGTCGATACCGGTATCGGCGTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGACCTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCTCGCCGATGGTGTTTTCGTca              <  2:1534199/139‑1 (MQ=11)
                                                                                                   cATGTCTATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCGTAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATc  >  1:1391533/1‑139 (MQ=255)
                                                                                                         |                                                                                                                                    
TGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCAGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATC  >  NC_000913/1571227‑1571464

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