Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,792,109 Δ1 bp coding (317/1320 nt) uhpC → membrane protein regulates uhpT expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,792,1070G.95.1% 148.0 / 2.9 41coding (315/1320 nt)uhpCmembrane protein regulates uhpT expression
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (11/28);  total (12/29)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGCT  >  W3110S.gb/3792054‑3792163
                                                     |                                                        
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:935621/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:3048642/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:2883515/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:2590977/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:376951/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:485514/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:209927/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:2079032/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:1879600/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:917385/1‑71 (MQ=255)
gTCAGCGATCGCTCAAATGCCCGTTATTATATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                        >  1:711840/1‑71 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:604253/47‑1 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:623118/47‑1 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:633457/47‑1 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:2810122/47‑1 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:1467843/47‑1 (MQ=255)
               aaTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:3168950/47‑1 (MQ=255)
                aTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:1884012/46‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:3630699/43‑1 (MQ=39)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:988098/43‑1 (MQ=39)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:1621104/43‑1 (MQ=39)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                 <  1:1625787/43‑1 (MQ=39)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:413015/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:1131114/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:1254737/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:672186/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:1815382/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:3103177/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:3819512/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:273809/71‑1 (MQ=255)
                   cccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:3253490/71‑1 (MQ=255)
                    ccGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:4063095/71‑1 (MQ=255)
                    ccGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGc                     <  1:2102070/70‑1 (MQ=255)
                                     tAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGt     >  1:4000759/1‑70 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTGGCTATGGGCGTTTGCCGTGct  <  1:151795/70‑1 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGct  <  1:2460160/70‑1 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGct  <  1:3727598/70‑1 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGct  <  1:3628700/70‑1 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGct  <  1:731510/70‑1 (MQ=255)
                                       gggCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGct  <  1:35678/70‑1 (MQ=255)
                                           ttATCGCCACGGCATTATCAACATTCTGTTTGGCt                                 <  1:1694026/35‑1 (MQ=37)
                                                     |                                                        
GTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGCT  >  W3110S.gb/3792054‑3792163

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