Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,570,650 T→C T466A (ACC→GCC)  gadB ← glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,6500TC94.9% 128.7 / ‑0.8 39T466A (ACC→GCC) gadBglutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (20/17);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

AACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTC  >  NC_000913/1570503‑1570797
                                                                                                                                                   |                                                                                                                                                   
aaCGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATg                                                                                                                                                  <  2:120461/151‑1 (MQ=255)
aaCGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATg                                                                                                                                                  <  1:66070/151‑1 (MQ=255)
   ggTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGttt                                                                                                                                               >  1:467655/1‑151 (MQ=255)
          cTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGctg                                                                                                                                        >  1:264206/1‑151 (MQ=255)
           tGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGctgt                                                                                                                                       >  2:657687/1‑151 (MQ=255)
               ggAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGctgttctg                                                                                                                                   >  2:627532/1‑151 (MQ=255)
                        ccTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAAt                                                                                                                          <  1:501930/151‑1 (MQ=255)
                                 ttCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGc                                                                                                                   >  1:262101/1‑149 (MQ=255)
                                     gAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTc                                                                                                             <  2:71418/151‑1 (MQ=255)
                                     gAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTc                                                                                                             >  1:553341/1‑151 (MQ=255)
                                        gCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCggg                                                                                                          <  1:161750/151‑1 (MQ=255)
                                           tATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTg                                                                                                        >  2:79576/1‑150 (MQ=255)
                                            aTTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGa                                                                                                       <  1:132701/150‑1 (MQ=255)
                                            aTTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGAt                                                                                                      <  2:320748/151‑1 (MQ=255)
                                                aGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCt                                                                                                  <  1:137693/151‑1 (MQ=255)
                                                 gAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTg                                                                                                 >  1:153762/1‑151 (MQ=255)
                                                   aCAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTgag                                                                                               <  2:554404/151‑1 (MQ=255)
                                                           ttCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCACGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTc                                                                                        >  2:450113/1‑150 (MQ=255)
                                                             ccGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAgg                                                                                     >  2:589623/1‑151 (MQ=255)
                                                             ccGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAgg                                                                                     <  2:581610/151‑1 (MQ=255)
                                                                        gCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGt                                                                           <  2:465079/150‑1 (MQ=255)
                                                                        gCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTa                                                                          <  1:244989/151‑1 (MQ=255)
                                                                                aCAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCa                                                                  <  2:304901/151‑1 (MQ=255)
                                                                                     cTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAAc                                                             >  1:241734/1‑151 (MQ=255)
                                                                                      tttGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAAc                                                             >  1:360262/1‑150 (MQ=255)
                                                                                      tttGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATAACCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAAc                                                             >  2:283939/1‑150 (MQ=255)
                                                                                        tGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAg                                                           <  2:390354/150‑1 (MQ=255)
                                                                                          gAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTc                                                        <  2:524183/151‑1 (MQ=255)
                                                                                              ttGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGc                                                     >  2:205721/1‑150 (MQ=255)
                                                                                                      aaaCGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCAt                                            <  1:405852/151‑1 (MQ=255)
                                                                                                       aaCGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCAtt                                           <  2:4373/151‑1 (MQ=255)
                                                                                                                       gACACCGTTACCGTTAAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCgacgac                           <  1:625737/151‑1 (MQ=255)
                                                                                                                            cgttacctttaCACGTTATCAGGCATGTTTCAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCACTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAAGAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGacacct                       >  2:72815/1‑148 (MQ=14)
                                                                                                                             gttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATaa                     <  2:394753/151‑1 (MQ=255)
                                                                                                                               taccgttaAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATg                   >  2:34188/1‑151 (MQ=255)
                                                                                                                                 ccgttaAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATgcg                 >  2:226847/1‑151 (MQ=255)
                                                                                                                                     taAACGTTATCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATc             >  2:497075/1‑151 (MQ=37)
                                                                                                                                             aTCAGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGAt     >  2:500842/1‑151 (MQ=18)
                                                                                                                                                aGGCATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTc  >  1:281168/1‑151 (MQ=18)
                                                                                                                                                   |                                                                                                                                                   
AACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTC  >  NC_000913/1570503‑1570797

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