Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,057,926 G→A 77.3% R151W (CGG→TGG)  torR ← DNA‑binding response regulator in two‑component regulatory system with TorS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,057,9260GA77.3% 35.7 / 9.2 22R151W (CGG→TGG) torRDNA‑binding response regulator in two‑component regulatory system with TorS
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (4/13);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GGTCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGTT  >  W3110S.gb/1057900‑1057980
                          |                                                      
ggTCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGc            <  1:779621/71‑1 (MQ=255)
ggTCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGc            <  1:1352384/71‑1 (MQ=255)
  tCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGa          >  1:701681/1‑71 (MQ=255)
  tCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGa          >  1:471123/1‑71 (MQ=255)
  tCAGTTTAATCGGATCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGa          >  1:91224/1‑71 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGa                                      >  1:472644/1‑35 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGtt  >  1:42083/1‑71 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGtt  >  1:1048079/1‑71 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:290128/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:814650/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:684498/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:493958/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:232734/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:213000/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:1929782/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:1588488/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:1347674/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:123181/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGt   <  1:1213201/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATACCAGt   <  1:1285186/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGCGACACATTCAGACAATAACCGGCAAAGCCATAGCAGt   <  1:1795567/70‑1 (MQ=255)
          aTCGGCTCGCCATCCCACTCCAGCGTATGGCGAGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGtt  >  1:507086/1‑71 (MQ=255)
                          |                                                      
GGTCAGTTTAATCGGCTCGCCATCCCGCTCCAGCGTATGGCGCGACACATTCAGGCAATAACCGGCAAAGCGATAGCAGTT  >  W3110S.gb/1057900‑1057980

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