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.G94.6% 125.5 / 0.1 37coding (262/729 nt)gadWtranscriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (18/17);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

TAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CAT‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664232‑3664479
                                                                                                                               |                                                                                                                          
tAAGGGGTTCAGCTCTTGCGTATGAGGAAATATCGATAAGAAGGGAAGGAAGAGCATATTGCTTATTATGGGGGAAGTTGTGCTGATTGGAGTGAAGCTATCAAAAATAGAATTCACCAACGGACGG‑GGGGGGGGACTg                                                                                                                >  2:838207/1‑139 (MQ=255)
         cAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGttaattaa                                                                                                       >  1:896478/1‑139 (MQ=255)
            ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    <  1:184736/139‑1 (MQ=255)
                 tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                               <  2:896478/139‑1 (MQ=255)
                    tATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCtt                                                                                            >  2:740540/1‑139 (MQ=255)
                      tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    <  1:138035/129‑1 (MQ=255)
                      tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACaa                                                                                                    >  2:138035/1‑129 (MQ=255)
                      tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  1:434758/1‑139 (MQ=255)
                      tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  1:464838/139‑1 (MQ=255)
                      tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  2:464838/1‑139 (MQ=255)
                        agaAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCt                                                                                        <  1:861890/139‑1 (MQ=255)
                           aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  2:434758/139‑1 (MQ=255)
                           aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     >  2:579047/1‑139 (MQ=255)
                             aCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   >  2:1053419/1‑139 (MQ=255)
                                         aGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   >  2:163266/1‑127 (MQ=255)
                                         aGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   <  1:163266/127‑1 (MQ=255)
                                           gAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGt                                                                     <  1:579047/139‑1 (MQ=255)
                                                    agCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGttaa                                                                                                           <  2:853904/92‑1 (MQ=255)
                                                    agCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGttaa                                                                                                           >  1:853904/1‑92 (MQ=255)
                                                     gCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CAT‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACg                                                          >  1:358306/1‑139 (MQ=255)
                                                         atTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAat                                                       >  2:780792/1‑139 (MQ=255)
                                                           tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                <  1:187830/128‑1 (MQ=255)
                                                           tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                >  2:187830/1‑128 (MQ=255)
                                                                       gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTg                                                                                                                >  2:143082/1‑68 (MQ=255)
                                                                       gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTg                                                                                                                <  1:143082/68‑1 (MQ=255)
                                                                       gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  2:169689/95‑1 (MQ=255)
                                                                       gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     >  1:169689/1‑95 (MQ=255)
                                                                       gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                         >  1:765296/1‑139 (MQ=255)
                                                                           aTTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCa                                     <  1:1053419/139‑1 (MQ=255)
                                                                            tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACaaa                                                                                                   <  2:237604/76‑1 (MQ=255)
                                                                            tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACaaa                                                                                                   >  1:237604/1‑76 (MQ=255)
                                                                            tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                    <  2:1028649/139‑1 (MQ=255)
                                                                                  ccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTc                              >  2:427471/1‑139 (MQ=255)
                                                                                        tCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCa                        <  1:427471/139‑1 (MQ=255)
                                                                                               aGCTATCAAAAATCGCCGTCACCAGCGG‑CAT‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATc                <  2:358306/139‑1 (MQ=255)
                                                                                                    tCAAAAATCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCCAACTGTCAGCTCATCGATATCCAgt            >  2:978539/1‑139 (MQ=255)
                                                                                                        aaaTCGCCGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgtc        <  2:951470/139‑1 (MQ=255)
                                                                                                              ccGTCACCAGCGG‑CATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAGTATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  >  2:718711/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                          
TAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGG‑CAT‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664232‑3664479

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