Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,660,197 G→A Y69Y (TAC→TAT wzxE ← O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,660,1970GA83.8% 83.8 / 10.1 37Y69Y (TAC→TATwzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (0/31);  total (2/35)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAAGATGCCAGCCCCGGCAAGCA  >  W3110S.gb/3660137‑3660249
                                                            |                                                    
aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                                     <  1:1115478/62‑1 (MQ=255)
aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                                     <  1:2757703/62‑1 (MQ=255)
aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                                     <  1:2427949/62‑1 (MQ=255)
aCCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGt                                                     <  1:2218724/62‑1 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:997850/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:881019/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:82967/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:751623/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:68975/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:6118/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2837803/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2164052/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2160371/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1837672/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1502683/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:121680/1‑47 (MQ=255)
              cTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1637444/1‑47 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2961274/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:269243/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1399143/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2406384/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2282060/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1560667/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:2818909/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1143209/1‑37 (MQ=255)
                        gCCGACCACGCGGCGCAGCTGTTGCGGATTATCATgg                                                      >  1:1908173/1‑37 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:484140/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:425422/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:3038254/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:746019/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:485767/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:49962/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:53251/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:58074/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1462854/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1018239/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:13203/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:781762/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:887472/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:895275/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1679351/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2879394/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1754642/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1667166/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:194889/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2048984/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2122203/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1652899/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:1538401/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2606908/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2606913/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2621660/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2668881/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2671901/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:14929/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:2854432/50‑2 (MQ=255)
                                   ggcgCAGCTGTTGGGGATTATCATGATACTGGGCAACGTATTTGGTTacg                              <  1:3018547/50‑2 (MQ=255)
                                                   aTTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAAGATGCCAGCCCCGGCAAGc   >  1:204515/1‑61 (MQ=255)
                                                   aTTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAAGATGCCAGCCCCGGCAAGCa  >  1:2508652/1‑62 (MQ=255)
                                                            |                                                    
ACCAAGTACCATCGCTGATGAAGTGCCGACCACGCGGCGCAGCTGTTGCGGATTATCATGGTACTGGGCAACGTATTTGGTTACACCGTTAAAGATGCCAGCCCCGGCAAGCA  >  W3110S.gb/3660137‑3660249

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: