Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,570,869 A→G Y393H (TAT→CAT)  gadB ← glutamate decarboxylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,8690AG95.0% 74.0 / ‑1.6 40Y393H (TAT→CAT) gadBglutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (17/21);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

AGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGCGGCAACCTG  >  NC_000913/1570733‑1571004
                                                                                                                                         |                                                                                                                                       
cagtaCAGTTCCAGCAAAGTCCATTCCGAAGCGTCAAAGAAACATAAGGCGTTTCGCCTGGATGGCGGGGCCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCa                                                                                                                                       <  2:519173/136‑1 (MQ=2)
   aaCAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGgt                                                                                                                                   <  1:212329/140‑1 (MQ=18)
    aCAGTTCAGCAAAGTCCATTTGGAAGCCGCGAGGCCAAAAAAGGCGAATCCCCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGgtg                                                                                                                                  <  1:578505/140‑1 (MQ=11)
        ttCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAgg                                                                                                                                     <  2:92440/133‑1 (MQ=17)
        ttCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAgg                                                                                                                                     >  1:92440/1‑133 (MQ=17)
        ttCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGgtctgtc                                                                                                                              >  1:458671/1‑135 (MQ=12)
        ttCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATc                                                                                                                              >  2:681028/1‑140 (MQ=17)
                   tCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTc                                                                                                                    >  1:519152/1‑140 (MQ=17)
                   tCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTc                                                                                                                    >  1:340984/1‑140 (MQ=17)
                    ccATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCa                                                                                                                   <  2:607863/140‑1 (MQ=18)
                       tttCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCa                                                                                                                <  1:857427/140‑1 (MQ=18)
                           gAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCtt                                                                                                            <  2:519152/140‑1 (MQ=18)
                                   gacgacACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTg                                                                                                    >  2:348411/1‑140 (MQ=17)
                                        cacaTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa                                                                                               <  1:702165/140‑1 (MQ=17)
                                                 cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCgg                                                                                      >  2:501538/1‑140 (MQ=17)
                                                   cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGTCGGGa                                                                                    >  1:382064/1‑140 (MQ=17)
                                                   cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGa                                                                                    >  1:618693/1‑140 (MQ=17)
                                                   cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGa                                                                                    >  1:646677/1‑140 (MQ=17)
                                                    aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGAt                                                                                   >  1:80643/1‑140 (MQ=17)
                                                     tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATg                                                                                  >  2:854948/1‑140 (MQ=17)
                                                       accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGcc                                                                                <  1:854948/140‑1 (MQ=17)
                                                                tCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCgg                                                                       <  1:577589/140‑1 (MQ=17)
                                                                tCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCgg                                                                       <  2:340984/140‑1 (MQ=17)
                                                                tCGGTGGCTTCACCGCCGAGAGTGAACGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCgg                                                                       <  1:604878/140‑1 (MQ=17)
                                                                     ggCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGAc                                                                  <  1:182492/140‑1 (MQ=17)
                                                                     ggCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGAc                                                                  <  1:715503/140‑1 (MQ=17)
                                                                     ggCTCCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGAc                                                                  <  2:176182/140‑1 (MQ=17)
                                                                        ttCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCg                                                               <  2:80643/140‑1 (MQ=17)
                                                                          cACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTa                                                             <  1:348411/140‑1 (MQ=17)
                                                                            ccgccgAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGc                                                                     >  2:173059/1‑130 (MQ=17)
                                                                            ccgccgAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGc                                                                     <  1:173059/130‑1 (MQ=17)
                                                                              gccgAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGa                                                         <  1:681028/140‑1 (MQ=17)
                                                                                 gagaGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGa                                                         <  1:80934/137‑1 (MQ=17)
                                                                                 gagaGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGa                                                         >  2:80934/1‑137 (MQ=17)
                                                                                    agTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACt                                                   >  2:399765/1‑140 (MQ=32)
                                                                                        aaGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATa                                               <  2:646677/140‑1 (MQ=17)
                                                                                                      ccAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTgg                                 >  2:301484/1‑140 (MQ=17)
                                                                                                                           ttCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAg            <  1:399765/140‑1 (MQ=32)
                                                                                                                            tCAGAGAGGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGc           >  2:747868/1‑140 (MQ=18)
                                                                                                                                  agGTCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGCGGCAAc     >  1:450273/1‑140 (MQ=18)
                                                                                                                                     tCATGCAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCTATTTCATCCGCCAGATAAGCGGCAACCTg  <  2:21062/140‑1 (MQ=14)
                                                                                                                                         |                                                                                                                                       
AGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGCGGCAACCTG  >  NC_000913/1570733‑1571004

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