Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,101,111 G→A 30.0% S272L (TCG→TTG)  csgG ← outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,101,1110GA30.0% 56.6 / 29.8 40S272L (TCG→TTG) csgGouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base A (7/5);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAC  >  W3110S.gb/1101049‑1101177
                                                              |                                                                  
tAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACatat                                                              >  1:3152183/1‑69 (MQ=255)
tAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACatat                                                              >  1:2048670/1‑69 (MQ=255)
tAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACatat                                                              >  1:370122/1‑69 (MQ=255)
     aGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCgg                                                         <  1:709657/69‑1 (MQ=255)
     aGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCgg                                                         <  1:2337214/69‑1 (MQ=255)
       tAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACa                                                                 <  1:3233668/59‑1 (MQ=255)
       tAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACa                                                                 <  1:736657/59‑1 (MQ=255)
        aaaaTGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACa                                                                 <  1:2345315/58‑1 (MQ=255)
        aaaaTGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACATATGGCGGTAt                                                      <  1:547795/69‑1 (MQ=255)
         aaaTGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                      >  1:3378560/1‑68 (MQ=255)
         aaaTGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                      >  1:3568368/1‑68 (MQ=255)
                   tgatAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATg                                                             >  1:485661/1‑51 (MQ=255)
                   tgatAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAAt                                           >  1:829443/1‑69 (MQ=255)
                      tAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACATATGGCGGTATTTCACCAGAATGTc                                        >  1:1787805/1‑69 (MQ=255)
                       aaTTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCa                                       <  1:993420/69‑1 (MQ=255)
                           ccGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACa                                                                 <  1:1763164/39‑1 (MQ=255)
                           ccGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACATATGGCGGTATTTCACCAGAATGTCATTCt                                   >  1:2061434/1‑69 (MQ=255)
                           ccGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCAACATATGGCGGTATTTCACCAGAATGTCATTCt                                   >  1:1702395/1‑69 (MQ=255)
                               cTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCg                               <  1:1532451/69‑1 (MQ=255)
                                 tttttATCTGTCAGGATTCCGGTAGAACCAACATATGGCGGTATTTCACCAGAATGt                                         >  1:1733693/1‑57 (MQ=255)
                                        cTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCttt                      >  1:2711881/1‑69 (MQ=255)
                                           tCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtt                   >  1:1386454/1‑69 (MQ=255)
                                             aGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAttt                  >  1:295912/1‑68 (MQ=255)
                                             aGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAttt                  >  1:3030313/1‑68 (MQ=255)
                                             aGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAttt                  >  1:1435171/1‑68 (MQ=255)
                                             aGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAttt                  >  1:1540888/1‑68 (MQ=255)
                                             aGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtt                   >  1:822580/1‑67 (MQ=255)
                                                aTTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCa              <  1:3662631/69‑1 (MQ=255)
                                                 ttCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCaa             >  1:456011/1‑69 (MQ=255)
                                                 ttCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCaa             >  1:811405/1‑69 (MQ=255)
                                                      gTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATcc         <  1:422152/68‑1 (MQ=255)
                                                      gTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATcc         <  1:2256089/68‑1 (MQ=255)
                                                        ggAACCGACATATGGCGGTATTTCACAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca     <  1:1228056/69‑1 (MQ=255)
                                                         gAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca     <  1:2223364/69‑1 (MQ=255)
                                                         gAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca     <  1:2943961/69‑1 (MQ=255)
                                                           aCCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAg    <  1:61247/68‑1 (MQ=255)
                                                           aCCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAg    >  1:557963/1‑68 (MQ=255)
                                                            ccGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAc  <  1:2299554/69‑1 (MQ=255)
                                                            ccGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAc  <  1:2856685/69‑1 (MQ=255)
                                                             cGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCa              <  1:1884708/56‑1 (MQ=255)
                                                              |                                                                  
TAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAC  >  W3110S.gb/1101049‑1101177

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: