Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,568,848 C→T A408T (GCC→ACC)  gadB ← glutamate decarboxylase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,8480CT98.4% 145.9 / ‑4.9 64A408T (GCC→ACC) gadBglutamate decarboxylase isozyme
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/34);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAA  >  NC_000913/1568799‑1568937
                                                 |                                                                                         
aTAATGCGCTTCACCACGATGTCGTTGCCTGCCCCGCCGAGAGTGAAGGTCGAAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                                <  2:2879139/93‑1 (MQ=11)
aTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                                <  2:2455123/93‑1 (MQ=14)
 tAATGCGCATCACCACGATGTCGTTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               <  1:4468226/93‑1 (MQ=12)
 tAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               >  1:620215/1‑93 (MQ=17)
 tAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               <  1:1634773/93‑1 (MQ=17)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  1:4371449/93‑1 (MQ=18)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  2:808296/93‑1 (MQ=18)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  2:387711/93‑1 (MQ=18)
  aaTGCGCACCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  2:4226923/93‑1 (MQ=16)
    tGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACa                                            <  2:3615719/93‑1 (MQ=18)
      cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                          <  1:5391/93‑1 (MQ=21)
      cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                          >  1:3605066/1‑93 (MQ=18)
      cgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAgg                                          >  2:2924867/1‑93 (MQ=18)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         >  1:763847/1‑93 (MQ=18)
       ccATCACCCCGATGAGGATGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         <  2:2426606/92‑1 (MQ=11)
       ccATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         <  1:2874488/92‑1 (MQ=21)
        cATCACCCCGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        <  1:1729262/93‑1 (MQ=255)
        cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        <  2:2278939/93‑1 (MQ=18)
        cATCACCACGAGGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        <  1:111683/93‑1 (MQ=255)
          tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                      >  2:1665892/1‑93 (MQ=18)
            accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTAt                                    >  2:2378345/1‑93 (MQ=18)
             ccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATc                                   <  1:873167/93‑1 (MQ=21)
              cacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATcc                                  >  1:245340/1‑93 (MQ=18)
              cacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGGGGATcc                                  >  2:1361814/1‑93 (MQ=255)
               aggATGTTGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                 <  1:1081151/91‑1 (MQ=255)
                ggATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCg                                <  2:3319230/92‑1 (MQ=18)
                 gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                               >  1:1098371/1‑93 (MQ=18)
                 gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                               <  1:3336167/93‑1 (MQ=21)
                  aTGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGa                              <  1:510548/93‑1 (MQ=21)
                  aTGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGa                              >  2:638024/1‑93 (MQ=18)
                  aTGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGa                              <  2:1741698/93‑1 (MQ=21)
                    gTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATc                            <  2:2655990/93‑1 (MQ=21)
                     tCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCt                           <  1:3114946/93‑1 (MQ=21)
                      cGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCtt                          <  1:3111302/93‑1 (MQ=18)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         >  1:58504/1‑93 (MQ=18)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         <  1:616007/93‑1 (MQ=21)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGGGTATCCCGGATCTTc                         >  1:2132061/1‑93 (MQ=16)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        >  2:434342/1‑93 (MQ=18)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        >  2:4451204/1‑93 (MQ=18)
                         tGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAc                       <  2:2063863/93‑1 (MQ=21)
                         tGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAc                       >  2:2420262/1‑93 (MQ=18)
                          ggCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGAGTATCACGGATCTTCAcc                      >  1:3421867/1‑93 (MQ=255)
                           gCTTCACCGCCGAGAGTGGAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGACATACAGGGTGTATCTCGGATCTTCACCa                     >  2:1436529/1‑93 (MQ=11)
                             ttCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATc                   >  1:3901957/1‑93 (MQ=18)
                             ttCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATc                   >  2:444266/1‑93 (MQ=18)
                              tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                  <  1:2609153/93‑1 (MQ=18)
                               ccccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                 <  2:1220901/91‑1 (MQ=16)
                               cACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                 <  1:4519001/93‑1 (MQ=21)
                               cACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                 >  2:4631299/1‑93 (MQ=18)
                               cACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                 <  1:2232226/93‑1 (MQ=18)
                               cACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCtt                 <  1:2275842/93‑1 (MQ=21)
                                 ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTc               >  1:4269317/1‑93 (MQ=18)
                                 ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTc               <  2:1916528/93‑1 (MQ=21)
                                  cgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCa              <  2:3372249/93‑1 (MQ=21)
                                    ccgAGAGTGAAGGTCGGAACCTGCCCGCCGCGGAGACGCAGACGTTCAGAGAGGGCATACAGGGTGTATCCCGGGTCTTCCCCATCTTTCAGt            >  2:2038448/1‑93 (MQ=11)
                                    ccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGt            >  2:2068852/1‑93 (MQ=18)
                                    ccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGt            >  2:1406080/1‑93 (MQ=18)
                                     cgAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           >  1:506018/1‑93 (MQ=33)
                                        gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa        >  2:3657306/1‑93 (MQ=18)
                                        gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa        >  2:702518/1‑93 (MQ=18)
                                          gTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAg      >  2:2886000/1‑93 (MQ=18)
                                           tGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGc     >  1:3745377/1‑93 (MQ=18)
                                            gAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa    <  1:324456/93‑1 (MQ=21)
                                             aaGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCCCCATCTTTCAGTTTGAAGCaa   <  2:2097216/93‑1 (MQ=16)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  >  1:4477155/1‑93 (MQ=18)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  <  1:3553412/93‑1 (MQ=21)
                                                 |                                                                                         
ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAA  >  NC_000913/1568799‑1568937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: