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.G96.6% 206.1 / ‑0.5 59coding (262/729 nt)gadWtranscriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (29/28);  total (30/29)
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.99e-01

TGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAA  >  NC_000913/3664225‑3664487
                                                                                                                                     |                                                                                                                                   
tGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggg                                                                                                                              <  2:813692/139‑1 (MQ=255)
  aaaaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggt                                                                                                                            <  2:104223/139‑3 (MQ=255)
    aaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggtga                                                                                                                          <  1:391846/139‑5 (MQ=255)
                 aGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAAc                                                                                                              >  2:984699/1‑139 (MQ=255)
                   ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            <  1:984699/139‑1 (MQ=255)
                   ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            <  2:198723/139‑1 (MQ=255)
                   ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            <  1:1427772/139‑1 (MQ=255)
                   ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            <  1:1188174/139‑1 (MQ=255)
                        tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                       >  1:1307295/1‑139 (MQ=255)
                           tATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCtt                                                                                                    >  2:1370818/1‑139 (MQ=255)
                           tATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCtt                                                                                                    <  1:1370818/139‑1 (MQ=255)
                             tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                                  >  1:429339/1‑139 (MQ=255)
                             tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                                  >  1:28670/1‑139 (MQ=255)
                             tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                                <  2:535839/108‑1 (MQ=255)
                             tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                                >  1:535839/1‑108 (MQ=255)
                             ggagaAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                                  >  1:1042651/2‑139 (MQ=255)
                                                    aaGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTgg                                                                           >  2:861693/1‑139 (MQ=255)
                                                       gTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                        <  1:1009780/139‑1 (MQ=255)
                                                               tatTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGaa                                                                <  2:266869/139‑1 (MQ=255)
                                                                  tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAatat                                                             >  1:1208564/1‑139 (MQ=255)
                                                                    cTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCt                                                                                                     >  1:659579/1‑97 (MQ=255)
                                                                    cTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCt                                                                                                     <  2:659579/97‑1 (MQ=255)
                                                                         tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                                      >  1:1148767/1‑75 (MQ=255)
                                                                         tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                                      <  2:1148767/75‑1 (MQ=255)
                                                                              gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             <  2:1396809/95‑1 (MQ=255)
                                                                              gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             >  1:1396809/1‑95 (MQ=255)
                                                                              gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                                 >  2:1622013/1‑139 (MQ=255)
                                                                              gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATggg                                                                                                                                 >  2:354060/1‑58 (MQ=255)
                                                                              gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATggg                                                                                                                                 <  1:354060/58‑1 (MQ=255)
                                                                                 aaTTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGc                                              <  1:1622013/139‑1 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                                      <  1:403839/50‑1 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                                      >  2:403839/1‑50 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            <  2:1340953/75‑1 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                            >  1:1340953/1‑75 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                        >  2:225570/1‑127 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                        <  1:225570/127‑1 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                            >  2:1232416/1‑139 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                            >  2:1019896/1‑139 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                            >  1:699244/1‑139 (MQ=255)
                                                                                   tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                            >  1:1224158/1‑139 (MQ=255)
                                                                                        cccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgtaagta                                                   <  2:770390/127‑1 (MQ=255)
                                                                                        cccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgtaagta                                                   >  1:770390/1‑127 (MQ=255)
                                                                                         ccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             >  2:1224900/1‑84 (MQ=255)
                                                                                         ccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             <  1:1224900/84‑1 (MQ=255)
                                                                                           ggATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAg                                    <  2:28670/139‑1 (MQ=255)
                                                                                             aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                                  <  2:1208564/139‑1 (MQ=255)
                                                                                             aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                                  <  1:861693/139‑1 (MQ=255)
                                                                                                 aTTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATc                              >  2:1030028/1‑139 (MQ=255)
                                                                                                       gCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATc                        <  2:1042651/139‑1 (MQ=255)
                                                                                                          aTCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                           <  1:734614/69‑1 (MQ=255)
                                                                                                          aTCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                           >  2:734614/1‑69 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             >  2:293178/1‑66 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                             <  1:293178/66‑1 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCt                                                                      <  1:650856/89‑1 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCt                                                                      >  2:650856/1‑89 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGaa                                                                <  2:722231/95‑1 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGaa                                                                >  1:722231/1‑95 (MQ=255)
                                                                                                           tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgt                    >  1:1445041/1‑139 (MQ=255)
                                                                                                              aaaaTCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt                 <  1:1030028/139‑1 (MQ=255)
                                                                                                                  tCGCCGTCACCAGCGGCATGGAGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCa                                             >  2:1585970/1‑107 (MQ=255)
                                                                                                                  tCGCCGTCACCAGCGGCATGGAGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCa                                             <  1:1585970/107‑1 (MQ=255)
                                                                                                                     ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                       <  2:24757/46‑1 (MQ=255)
                                                                                                                     ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                       <  1:977375/46‑1 (MQ=255)
                                                                                                                     ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                       >  1:24757/1‑46 (MQ=255)
                                                                                                                     ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                       >  2:977375/1‑46 (MQ=255)
                                                                                                                         cACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACg                                                                  <  2:1237202/79‑1 (MQ=255)
                                                                                                                         cACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACg                                                                  >  1:1237202/1‑79 (MQ=255)
                                                                                                                             aGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGaa  <  2:429339/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
TGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAA  >  NC_000913/3664225‑3664487

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