Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,792,107 Δ1 bp 100% coding (315/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.76.8% 93.7 / 58.5 56coding (315/1320 nt)uhpCmembrane protein regulates uhpT expression
Reads supporting (aligned to +/- strand):  ref base G (4/9);  new base . (10/33);  total (14/42)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGTCTCCGGCATTGTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGCTCTG  >  W3110S.gb/3792041‑3792166
                                                                  |                                                           
tGTCTCCGGCATTGTCAGCGATCGCACAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACggg                                                           >  1:2386357/1‑69 (MQ=255)
  tctcCGGCATTGTCAGCGATCGCTCAAATGCCCGTTATTTTATTGGGATAGGGCTTATCGCCACggca                                                          <  1:3557115/68‑3 (MQ=255)
  tctcCGGCATTGTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACggca                                                          <  1:41235/68‑3 (MQ=255)
         cATTGTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAAc                                                 <  1:834080/69‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                    <  1:2659924/62‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                    <  1:3065638/62‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                    <  1:3077617/62‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                    <  1:3114571/62‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATc                                                    <  1:2236692/62‑1 (MQ=255)
             gTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATAACATTCt                                            >  1:1166246/1‑69 (MQ=255)
                aGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                           <  1:988514/68‑1 (MQ=255)
                aGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTg                                           <  1:2231336/68‑1 (MQ=255)
                aGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGt                                          >  1:1051953/1‑69 (MQ=255)
                  cGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGt                                          <  1:2102723/68‑1 (MQ=255)
                  cGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGt                                          <  1:2107170/68‑1 (MQ=255)
                    aTCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTgg                                      <  1:626791/69‑1 (MQ=255)
                          cAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:1515139/69‑1 (MQ=255)
                          cAAATGCCCGTTATTTTATGGGGATAGGGCTTATAGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:1017676/69‑1 (MQ=255)
                             aTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:2585334/67‑1 (MQ=255)
                             aTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:1772886/67‑1 (MQ=255)
                             aTGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:2993441/69‑1 (MQ=255)
                              tGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             >  1:3588727/1‑68 (MQ=255)
                              tGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACg                            >  1:831722/1‑69 (MQ=255)
                              tGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACg                            >  1:2206964/1‑69 (MQ=255)
                                    ttATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTa                      <  1:1104888/69‑1 (MQ=255)
                                    ttATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTa                      <  1:2385935/69‑1 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCa                                                   <  1:551120/39‑1 (MQ=38)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                >  1:1651436/1‑58 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                >  1:2291213/1‑58 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                >  1:2184331/1‑58 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:2163998/58‑1 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:1753723/58‑1 (MQ=255)
                                     tATTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtct                                 >  1:653297/1‑57 (MQ=255)
                                     tATTTTATGGGGATAGGGATTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                >  1:1920837/1‑58 (MQ=255)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTAt                     >  1:2255316/1‑69 (MQ=255)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGTCATTATCaa                                                  <  1:3470354/39‑1 (MQ=37)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGCATTATCaa                                                  <  1:1018050/39‑1 (MQ=38)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGCATTATCaa                                                  <  1:1285571/39‑1 (MQ=38)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:2479085/57‑1 (MQ=255)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTtctc                                <  1:723255/57‑1 (MQ=255)
                                      aTTTTATGGGGATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTAt                     <  1:62691/68‑1 (MQ=255)
                                              gggATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGttt             <  1:2830031/69‑1 (MQ=255)
                                              gggATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGttt             <  1:2834260/69‑1 (MQ=255)
                                              gggATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGttt             <  1:3240770/69‑1 (MQ=255)
                                              gggATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGttt             <  1:35216/69‑1 (MQ=255)
                                              gggATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGttt             <  1:1611020/69‑1 (MQ=255)
                                                gATAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGc           <  1:3431123/68‑1 (MQ=255)
                                                 aTAGGGCTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCg         <  1:1426867/69‑1 (MQ=255)
                                                  tAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGcc          >  1:1124954/1‑68 (MQ=255)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:111723/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:1196412/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:1696262/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:1700786/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:1743633/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:2350423/43‑1 (MQ=38)
                                                       cTTATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGAc                             <  1:2486826/43‑1 (MQ=38)
                                                        ttATCGCCACGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGCTCTg  <  1:2600699/69‑1 (MQ=255)
                                                            cGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTa                      >  1:147073/1‑46 (MQ=255)
                                                            agCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTa                      >  1:2378145/2‑46 (MQ=255)
                                                                  |                                                           
TGTCTCCGGCATTGTCAGCGATCGCTCAAATGCCCGTTATTTTATGGGGATAGGGCTTATCGCCACGGGCATTATCAACATTCTGTTTGGCTTCTCGACGTCGCTATGGGCGTTTGCCGTGCTCTG  >  W3110S.gb/3792041‑3792166

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: