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.2% 124.9 / ‑3.6 36coding (262/729 nt)gadWtranscriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (17/18);  total (18/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 = 1.00e+00

TATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGG  >  NC_000913/3664222‑3664483
                                                                                                                                         |                                                                                                                              
tATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATggg                                                                                                                             <  2:556401/139‑1 (MQ=255)
 aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggg                                                                                                                            >  1:69165/1‑139 (MQ=255)
   tGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggg                                                                                                                          <  2:6370/139‑1 (MQ=255)
        aGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggtgac                                                                                                                     <  2:69165/139‑6 (MQ=255)
            aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                                  <  2:133048/139‑1 (MQ=255)
            aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                                  <  1:719566/139‑1 (MQ=255)
            aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGt                                                                                                                  >  1:464068/1‑139 (MQ=255)
                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGt                                                                                                   <  1:639080/139‑1 (MQ=255)
                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATggg                                                                                                                             >  1:152812/1‑112 (MQ=255)
                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATggg                                                                                                                             <  2:152812/112‑1 (MQ=255)
                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                              <  1:550507/139‑1 (MQ=255)
                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                              >  1:423938/1‑139 (MQ=255)
                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                              >  1:232838/1‑139 (MQ=255)
                                  agaAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCt                                                                                            >  1:97914/1‑139 (MQ=255)
                                     aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         <  2:423938/139‑1 (MQ=255)
                                     aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         <  1:635378/139‑1 (MQ=255)
                                                          gTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                    <  2:232838/139‑1 (MQ=255)
                                                                   atTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAat                                                           <  2:341593/139‑1 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                                 >  2:99550/1‑68 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                                 <  1:99550/68‑1 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAata                                                          <  2:645458/139‑1 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAata                                                          >  2:341998/1‑139 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGACATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTgtgt                                                                                        >  1:721746/1‑109 (MQ=255)
                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGACATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTgtgt                                                                                        <  2:721746/109‑1 (MQ=255)
                                                                     tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggt                                                                                                                        <  1:704644/75‑3 (MQ=255)
                                                                     tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATgggggggt                                                                                                                        >  2:704644/1‑73 (MQ=255)
                                                                       cTCAGTATCGGTGAATTTACCTCGGCTTCATTGAAGCTATCAAAAATCGCCGTCACCAGTGGCACGGGGGGGTGACTCGGTGATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGTTAGCTGACGAATAttt                                                       >  2:607317/1‑139 (MQ=255)
                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                        >  2:307375/1‑80 (MQ=255)
                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACaa                                                                                                        <  1:307375/80‑1 (MQ=255)
                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                             >  2:531772/1‑139 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         >  2:147726/1‑90 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         >  1:226849/1‑90 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         <  2:226849/90‑1 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                         <  1:147726/90‑1 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                       >  2:1928/1‑92 (MQ=255)
                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                       <  1:1928/92‑1 (MQ=255)
                                                                                               gATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGtta                                                                                                                <  2:86920/58‑1 (MQ=255)
                                                                                               gATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGtta                                                                                                                >  1:86920/1‑58 (MQ=255)
                                                                                                              tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGt                                                                         >  2:395627/1‑82 (MQ=255)
                                                                                                              tCAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGt                                                                         <  1:395627/82‑1 (MQ=255)
                                                                                                               cAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtc               >  2:585913/1‑139 (MQ=255)
                                                                                                               cAAAAATCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtc               <  1:585913/139‑1 (MQ=255)
                                                                                                                     tCGCCGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTcgcg         >  1:554108/1‑139 (MQ=255)
                                                                                                                        ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCg                                                                          >  1:220440/1‑71 (MQ=255)
                                                                                                                        ccGTCACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCg                                                                          <  2:220440/71‑1 (MQ=255)
                                                                                                                            cACCAGCGGCATGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCgg  >  1:385541/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                              
TATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGG  >  NC_000913/3664222‑3664483

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