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.7% 97.5 / ‑3.1 30coding (262/729 nt)gadWtranscriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (13/16);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

ATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664223‑3664479
                                                                                                                                       |                                                                                                                          
aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                        >  2:367541/1‑139 (MQ=255)
           aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                              >  2:173674/1‑139 (MQ=255)
           aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                              >  2:696341/1‑139 (MQ=255)
                       ctttcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAAc                                                                                                  >  1:115816/1‑139 (MQ=255)
                               tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  2:393791/139‑1 (MQ=255)
                               tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  2:337210/1‑139 (MQ=255)
                                 agaAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCt                                                                                        >  1:452433/1‑139 (MQ=255)
                                    aaaCATCGATAAGCTGGTTAGGTTGGGCATATTGCTCAGTATCGTTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCGCCAGCGATATGGGGGGGGGACTGGTTAATTAACAAACGGTCTTGGCtgtg                                                                                     <  1:173674/139‑1 (MQ=255)
                                    aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  1:653886/139‑1 (MQ=255)
                                      aCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   <  1:337210/139‑1 (MQ=255)
                                      aCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACAGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   >  2:619885/1‑139 (MQ=255)
                                                         gTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                <  2:503814/139‑1 (MQ=255)
                                                            gagCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCATCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTg                                                             >  2:295512/1‑139 (MQ=255)
                                                               cATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACg                                                          <  2:115816/139‑1 (MQ=255)
                                                                  atTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAat                                                       <  2:639036/139‑1 (MQ=255)
                                                                  atTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAata                                                      >  2:549246/1‑139 (MQ=255)
                                                                      cTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATAttt                                                   >  1:640271/1‑139 (MQ=255)
                                                                              cGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAAcc                                                                                <  1:120050/102‑1 (MQ=255)
                                                                              cGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAAcc                                                                                >  2:120050/1‑102 (MQ=255)
                                                                                gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                         >  2:315893/1‑139 (MQ=255)
                                                                                gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                         >  1:653483/1‑139 (MQ=255)
                                                                                   aaTTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACTGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGc                                      <  1:545022/139‑1 (MQ=255)
                                                                                        aCCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACt                                 <  1:367541/139‑1 (MQ=255)
                                                                                               aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                          <  2:452433/139‑1 (MQ=255)
                                                                                                 tCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCa                        <  2:793308/139‑1 (MQ=255)
                                                                                                         gCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGc                                      <  1:63270/117‑1 (MQ=255)
                                                                                                         gCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGc                                      >  2:63270/1‑117 (MQ=255)
                                                                                                             tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTATGTAATGGCAAACTGTCAGCTCATCGATATCCAgt            >  1:796208/1‑139 (MQ=255)
                                                                                                                aaaaTCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt         <  1:315893/139‑1 (MQ=255)
                                                                                                                       ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  <  1:350464/139‑1 (MQ=255)
                                                                                                                       ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  <  1:179076/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                          
ATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664223‑3664479

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