Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,056,565 G→A 50.0% G12G (GGC→GGT torS ← hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,056,5650GA50.0% ‑1.5 / 32.8 26G12G (GGC→GGTtorShybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base G (4/9);  new base A (8/5);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAGGTTGAG  >  W3110S.gb/1056508‑1056631
                                                         |                                                                  
cAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGc                                                                   <  1:1156394/59‑1 (MQ=255)
cAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGc                                                                   <  1:926265/59‑1 (MQ=255)
cAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGc                                                                   <  1:3631449/59‑1 (MQ=255)
cAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGc                                                                   <  1:2866485/59‑1 (MQ=255)
  ggTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGccc                                                                 >  1:1084478/1‑59 (MQ=255)
  ggTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGccc                                                                 >  1:3180971/1‑59 (MQ=255)
  ggTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGccc                                                                 >  1:2365626/1‑59 (MQ=255)
     tGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTc                                                       >  1:992873/1‑66 (MQ=255)
     tGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTc                                                       >  1:3648170/1‑66 (MQ=255)
     tGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTc                                                       >  1:2188266/1‑66 (MQ=255)
            tCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTCTTCGGGTCAggg                                           >  1:2982243/1‑71 (MQ=255)
               caccaGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTCTTCGGGTCAGGGTTa                                        <  1:2408256/71‑1 (MQ=255)
                                gggTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTCTTCGGGTc                                               <  1:3122561/47‑1 (MQ=255)
                                  gTTAACAGCGCCATCAGGGCAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACt                     >  1:2577019/1‑71 (MQ=255)
                                   ttAACAGCGCCATCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTa                                        <  1:1575338/51‑1 (MQ=255)
                                   ttAACAGCGCCATCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTa                                        <  1:281270/51‑1 (MQ=255)
                                             cATCAGGGCAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGaaaa          <  1:1902510/71‑1 (MQ=255)
                                               tCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGaaaa          <  1:1467082/69‑1 (MQ=255)
                                               tCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGaaaa          <  1:512975/69‑1 (MQ=255)
                                               tCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGaaaa          <  1:1699109/69‑1 (MQ=255)
                                                 aGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAg        >  1:2405257/1‑69 (MQ=255)
                                                   ggCAAAACCCATCCAGAGTCTTCGGGTCAGGGTTaaa                                      <  1:1996977/37‑1 (MQ=255)
                                                   ggCAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAGGTTg    <  1:1967065/71‑1 (MQ=255)
                                                     cAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTGACGGTCGGTGCACTTTAGGTGAAAAAGGTTGAg  >  1:1937842/1‑71 (MQ=255)
                                                     cAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAGGTTGAg  >  1:386197/1‑71 (MQ=255)
                                                     cAAAACCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAGGTTGAg  >  1:564307/1‑71 (MQ=255)
                                                         |                                                                  
CAGGTTGTACCATCCCACCAGGGTACTGGTCAGGGTTAACAGCGCCATCAGGGCAAAGCCCATCCAGAGTCTTCGGGTCAGGGTTAAATTCACGGTCGGTGCACTTTAGGTGAAAAAGGTTGAG  >  W3110S.gb/1056508‑1056631

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: