Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,664,363 (G)6→7 coding (256/729 nt) gadW ← transcriptional activator of gadA and gadBC, repressor of gadX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,664,3571.G97.7% 147.4 / ‑3.4 43coding (262/729 nt)gadWtranscriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (21/21);  total (22/21)
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.73e-01

CTATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664221‑3664479
                                                                                                                                          |                                                                                                                          
cTATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgg                                                                                                                          >  1:267003/1‑139 (MQ=255)
  aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                        >  1:730266/1‑139 (MQ=255)
  aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                        >  1:340982/1‑139 (MQ=255)
      aaaaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggt                                                                                                                    <  1:196702/139‑3 (MQ=255)
      aaaaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGAATATTGCTAGGTATCGTTGAATTTACCCCGGATTCATTAAAGCTATCAAAAACCGCCGTCACCAGCGGCATgggggggt                                                                                                                    >  1:504631/1‑137 (MQ=255)
        aaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggtga                                                                                                                  <  1:132656/139‑5 (MQ=255)
           tAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTg                                                                                                                >  1:682940/1‑139 (MQ=255)
            aaGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTgg                                                                                                               <  2:164812/139‑1 (MQ=255)
             aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATTGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                              >  2:150898/1‑139 (MQ=255)
             aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAGTTTACCCCGGATGCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAGGGGGGGGTGACTgga                                                                                                              >  1:347989/1‑138 (MQ=255)
                   tCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGttaatta                                                                                                        >  1:344042/1‑139 (MQ=255)
                     aGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAAc                                                                                                      <  1:744174/139‑1 (MQ=255)
                       ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    <  1:320493/139‑1 (MQ=255)
                         ctttcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAAc                                                                                                  <  2:168157/139‑1 (MQ=255)
                            tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                               <  2:340982/139‑1 (MQ=255)
                              ttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCt                                                                                             <  2:730266/139‑1 (MQ=255)
                               tATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCtt                                                                                            >  2:435248/1‑139 (MQ=255)
                                 tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  2:124093/1‑139 (MQ=255)
                                 tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  2:100967/139‑1 (MQ=255)
                                 tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  2:87906/139‑1 (MQ=255)
                                 tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  1:174565/1‑139 (MQ=255)
                                      aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  1:702143/139‑1 (MQ=255)
                                      aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  1:435248/139‑1 (MQ=255)
                                           tCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAAcc                                                                                <  2:267003/139‑1 (MQ=255)
                                                               agCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGa                                                            >  1:694175/1‑139 (MQ=255)
                                                                   tatTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGaa                                                        <  2:174565/139‑1 (MQ=255)
                                                                     ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAAc                                                                                 <  2:366073/112‑1 (MQ=255)
                                                                     ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAAc                                                                                 >  1:366073/1‑112 (MQ=255)
                                                                      tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCt                                                                              >  2:257559/1‑114 (MQ=255)
                                                                      tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCt                                                                              <  1:257559/114‑1 (MQ=255)
                                                                           aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                >  1:429171/1‑139 (MQ=255)
                                                                           aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                >  1:256661/1‑139 (MQ=255)
                                                                             tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCACAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGAGGTAGCTGACGAATATTTTgtaa                                              >  1:420435/1‑139 (MQ=255)
                                                                             tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgtaa                                              >  1:543383/1‑139 (MQ=255)
                                                                                  gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                               <  2:535448/85‑1 (MQ=255)
                                                                                  gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                               >  1:535448/1‑85 (MQ=255)
                                                                                  gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCg                                                                      <  1:738563/110‑1 (MQ=255)
                                                                                  gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCg                                                                      >  2:738563/1‑110 (MQ=255)
                                                                                  gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                         >  2:56910/1‑139 (MQ=255)
                                                                                       tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    <  1:333061/75‑1 (MQ=255)
                                                                                       tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    >  2:333061/1‑75 (MQ=255)
                                                                                            cccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  2:238633/1‑80 (MQ=255)
                                                                                            cccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  1:238633/80‑1 (MQ=255)
                                                                                                gATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGc                           >  2:413373/1‑139 (MQ=255)
                                                                                                 aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTg                                                 >  1:751583/1‑116 (MQ=255)
                                                                                                 aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTg                                                 <  2:751583/116‑1 (MQ=255)
                                                                                                 aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                          <  2:472323/139‑1 (MQ=255)
                                                                                                    cATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCAt                       <  2:256661/139‑1 (MQ=255)
                                                                                                                         ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  >  1:508661/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                          
CTATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664221‑3664479

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