Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,666,292 T→G E430A (GAA→GCA)  gadA ← glutamate decarboxylase A, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,666,2920TG98.0% 109.0 / ‑4.5 50E430A (GAA→GCA) gadAglutamate decarboxylase A, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (24/25);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AAGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCACCATCTTTC  >  NC_000913/3666145‑3666435
                                                                                                                                                   |                                                                                                                                               
aaGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGa                                                                                                                                              <  2:709884/151‑1 (MQ=255)
 aGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGaa                                                                                                                                             <  2:251985/151‑1 (MQ=255)
  ggCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAg                                                                                                                                            <  1:180684/151‑1 (MQ=255)
  ggCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAg                                                                                                                                            <  1:706907/151‑1 (MQ=255)
   gCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGc                                                                                                                                           <  2:836836/151‑1 (MQ=255)
    cTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGcc                                                                                                                                          >  2:837333/1‑151 (MQ=255)
            tCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGcgacga                                                                                                                                   <  2:95629/150‑1 (MQ=255)
            tCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCgacgac                                                                                                                                  >  1:248055/1‑151 (MQ=255)
                     aCAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATgc                                                                                                                         <  1:190566/151‑1 (MQ=255)
                      cAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATgcg                                                                                                                        <  1:532338/151‑1 (MQ=255)
                         cGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCAt                                                                                                                     >  2:776538/1‑151 (MQ=38)
                           ttatgttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATc                                                                                                                    >  2:167702/1‑150 (MQ=37)
                            tatgttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATcac                                                                                                                  >  1:94926/1‑151 (MQ=25)
                            tatgttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATca                                                                                                                   <  1:17804/150‑1 (MQ=25)
                               gttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATcaccac                                                                                                               <  1:59504/151‑1 (MQ=18)
                               gttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTCCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATcaccac                                                                                                               <  1:519738/151‑1 (MQ=18)
                                ttatCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACg                                                                                                              >  1:848463/1‑151 (MQ=18)
                                     aGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTc                                                                                                         <  1:466343/151‑1 (MQ=18)
                                     aGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTc                                                                                                         >  2:384006/1‑151 (MQ=18)
                                      ggTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCg                                                                                                        <  1:614326/151‑1 (MQ=18)
                                       gtgtgtTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCgg                                                                                                       <  2:9861/151‑1 (MQ=18)
                                          tgtTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTgg                                                                                                    >  2:100129/1‑151 (MQ=14)
                                               aaaGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCa                                                                                               <  1:689461/151‑1 (MQ=14)
                                                   cTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAccgc                                                                                           >  2:398061/1‑151 (MQ=14)
                                                   cTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAccgc                                                                                           >  2:450102/1‑151 (MQ=14)
                                                    tGTTCTACTAGGCAATAACCTGCAGTTTCGGGTGATCGCTGTGCTATTTCAGGGAGGCTTTGTAGACATCCAGCAACACCTCAGCAAAGTCCATGGCGAAGCCGCGACGCCCCATAATACGCATCAACACGAATTCGGCCGCtt                                                                                                 >  2:738826/1‑144 (MQ=2)
                                                       tCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAccgccg                                                                                         >  2:735611/1‑149 (MQ=14)
                                                          gctgGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAGGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCgag                                                                                       <  2:762659/148‑1 (MQ=14)
                                                          gctgGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGagag                                                                                     <  2:93098/150‑1 (MQ=14)
                                                                                                             cTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGt                                  >  1:710302/1‑150 (MQ=14)
                                                                                                               ttGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTa                                 <  1:856021/149‑1 (MQ=14)
                                                                                                                   aGTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGgt                           >  2:193813/1‑151 (MQ=14)
                                                                                                                    gTCTTCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGgtg                          <  1:644575/151‑1 (MQ=14)
                                                                                                                       ttCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTAt                       >  1:18262/1‑151 (MQ=14)
                                                                                                                        tCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATc                      >  2:95992/1‑151 (MQ=14)
                                                                                                                        tCCAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATc                      >  1:609366/1‑151 (MQ=14)
                                                                                                                          cAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATccc                    >  2:560065/1‑151 (MQ=14)
                                                                                                                          cAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATccc                    >  1:806190/1‑151 (MQ=14)
                                                                                                                          cAGCAACAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATc                      <  1:719024/149‑1 (MQ=14)
                                                                                                                              aaCAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCCCCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGAt                <  2:188880/151‑1 (MQ=14)
                                                                                                                               aCATTTCAGCAAAGTCCATTGCGAAGCCGCGACGTACCATAATGCGCATCACCACGATTTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCgg                  <  1:859870/148‑1 (MQ=11)
                                                                                                                               aCAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATc               <  1:900100/151‑1 (MQ=14)
                                                                                                                                cAGTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCt              <  1:703370/151‑1 (MQ=14)
                                                                                                                                 aattCAGCAAAAACCATTGCGAGGCCGCAACGACACATAAGGCGCAACACCAGGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCtt             <  1:734928/149‑1 (MQ=11)
                                                                                                                                  gTTCAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTc            >  2:853202/1‑151 (MQ=14)
                                                                                                                                     cAGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCAcc         >  2:879351/1‑151 (MQ=14)
                                                                                                                                      aGCAAAGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCAcc         >  1:871867/1‑150 (MQ=14)
                                                                                                                                         aaaGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCACCATCt     >  1:695666/1‑151 (MQ=14)
                                                                                                                                          aaGTCCATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCACCATCtt    <  2:543204/151‑1 (MQ=14)
                                                                                                                                              ccATTGCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCACCATCTTTc  >  2:523235/1‑149 (MQ=14)
                                                                                                                                                   |                                                                                                                                               
AAGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCGTACAGGGTGTATCCCGGATCTTCACCATCTTTC  >  NC_000913/3666145‑3666435

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