Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,101,395 G→A 10.2% R177R (CGC→CGT csgG ← outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,101,3950GA10.2% 132.9 / 5.0 49R177R (CGC→CGTcsgGouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

TCATAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTGCCCC  >  W3110S.gb/1101327‑1101463
                                                                    |                                                                    
tCATAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGAcgcg                                                                    >  1:3271864/1‑71 (MQ=255)
 cATAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGAcgcg                                                                    <  1:1600599/70‑1 (MQ=255)
 cATAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGAcgc                                                                     <  1:48599/69‑1 (MQ=255)
  aTAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCa                                                                  <  1:1902674/71‑1 (MQ=255)
  aTAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCa                                                                  <  1:1762179/71‑1 (MQ=255)
       aaaGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGCCGCGCAGGt                                                               >  1:1206657/1‑69 (MQ=255)
        aaGTATCGTCTTACTGGTGTTCACCGAGGAGAGGATCTCACCGGTACTCACATTGACGACACGCAGGTTCa                                                            >  1:2407433/1‑71 (MQ=255)
        aaGTATCGTCTTACTGGTGTTCACCGAGGAGAGGATCTCACCGGTACTCACATTGACGACACGCAGGTTCa                                                            >  1:2064372/1‑71 (MQ=255)
          gTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACg                                                          >  1:1569429/1‑71 (MQ=255)
           tATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAg                                                                 >  1:3424922/1‑63 (MQ=255)
            aTCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACgg                                                         >  1:1965651/1‑70 (MQ=255)
               gTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAg                                                                 >  1:2889059/1‑59 (MQ=255)
                      tGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCa                                                       <  1:1690788/62‑1 (MQ=255)
                        gtgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCaa                                                      >  1:3448478/1‑61 (MQ=255)
                        gtgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAAt                                                     >  1:170027/1‑62 (MQ=255)
                        gtgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGa                                                 >  1:3017996/1‑66 (MQ=255)
                        gtgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCg                                              <  1:1745874/69‑1 (MQ=255)
                         tgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGc                                           >  1:464098/1‑71 (MQ=255)
                         tgtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGc                                           >  1:327484/1‑71 (MQ=255)
                          gtTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCt                                          >  1:1452404/1‑71 (MQ=255)
                           ttCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCt                                          <  1:1702099/70‑1 (MQ=255)
                             cACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCa                                                            <  1:85671/50‑1 (MQ=255)
                             cACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGAt                                                >  1:3664889/1‑62 (MQ=255)
                             cACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTgg                                        >  1:910932/1‑70 (MQ=255)
                             cACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTgg                                        >  1:2013074/1‑70 (MQ=255)
                              aCCGAAGAAAGGATCTCGCCGGTACTCACATTGACGCCGCGCAGGTTCACGGCAATCTGATCGAGCTGGTa                                      >  1:2606300/1‑71 (MQ=255)
                              aCCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:2652531/70‑1 (MQ=255)
                                 gaagaaAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:1792847/67‑1 (MQ=255)
                                 gaagaaAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:1501965/67‑1 (MQ=255)
                                  aagaaAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGTCAATCTg                                                  <  1:3511222/55‑1 (MQ=255)
                                  aagaaAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTg                                                  <  1:2656628/55‑1 (MQ=255)
                                      aaGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGa                                             <  1:1988221/56‑1 (MQ=255)
                                       aGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTc                             <  1:1201106/71‑1 (MQ=255)
                                       aGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTc                             <  1:2199198/71‑1 (MQ=255)
                                        ggATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCgtgt                              >  1:1470246/1‑69 (MQ=255)
                                        ggATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCgtgt                              >  1:2033269/1‑69 (MQ=255)
                                           tctcGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCa                         >  1:1572147/1‑71 (MQ=255)
                                                   gTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:2412324/49‑1 (MQ=255)
                                                    tACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:158866/48‑1 (MQ=255)
                                                    tACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGcaa                <  1:673793/71‑3 (MQ=255)
                                                      cTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCaa               <  1:2602418/70‑1 (MQ=255)
                                                         acaTTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAata           >  1:2842451/1‑71 (MQ=255)
                                                         acaTTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAata           >  1:96629/1‑71 (MQ=255)
                                                         acaTTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAat            <  1:1605180/70‑1 (MQ=255)
                                                         acaTTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAat            >  1:92117/1‑70 (MQ=255)
                                                         acaTTGACGACACGCAGGTTCACGGCAATCTGATCAAGCTGGTAt                                     <  1:3135310/45‑1 (MQ=38)
                                                         acaTTGACGACACGCAGGTTCACGGCAATCTGATCAAGCTGGTAt                                     <  1:1249858/45‑1 (MQ=38)
                                                          caTTGACGACACGCAGGTTCACGGCAATCTGATCAAGCTGGTAt                                     <  1:1915196/44‑1 (MQ=38)
                                                            ttGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGt                                       <  1:399065/40‑1 (MQ=255)
                                                              gacgacGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCtt       <  1:3400981/70‑1 (MQ=255)
                                                              gacgacGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCtt       >  1:3187637/1‑70 (MQ=255)
                                                              gacgacGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTg      >  1:3439090/1‑71 (MQ=255)
                                                               acgacgCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCtt       <  1:2920769/69‑1 (MQ=255)
                                                                cgacgCGCAGGTTCACGGCAATCTGAGCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTg      <  1:1011729/69‑1 (MQ=255)
                                                                  acgCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTGcccc  <  1:591714/71‑1 (MQ=255)
                                                                    gcgcAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTg      >  1:521013/1‑65 (MQ=255)
                                                                    |                                                                    
TCATAGGAAAGTATCGTCTTACTGGTGTTCACCGAAGAAAGGATCTCGCCGGTACTCACATTGACGACGCGCAGGTTCACGGCAATCTGATCGAGCTGGTATTGCGTGTCGGCACCGATGCCAAAATATCTTGCCCC  >  W3110S.gb/1101327‑1101463

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