Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 775,958 G→T 10.3% P128P (CCG→CCT tolQ → membrane spanning protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb775,9580GT10.3% 146.0 / 2.8 48P128P (CCG→CCTtolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base T (2/3);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCAC  >  W3110S.gb/775890‑776027
                                                                    |                                                                     
gAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGtt                                                                     <  1:525437/71‑1 (MQ=255)
   ggggCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTcc                                                                   >  1:1650808/1‑70 (MQ=255)
     ggCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCg                                                                       <  1:1006763/64‑1 (MQ=255)
         tCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTAc                                                            >  1:764982/1‑71 (MQ=255)
         tCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTAc                                                            >  1:3228372/1‑71 (MQ=255)
           gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTc                                                                    >  1:2201320/1‑61 (MQ=255)
           gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCgg                                                               >  1:1887679/1‑66 (MQ=255)
            cgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                                                         <  1:2754076/71‑1 (MQ=255)
                   tGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   <  1:3015056/70‑1 (MQ=255)
                   tGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   >  1:3606344/1‑70 (MQ=255)
                   tGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   >  1:3118004/1‑70 (MQ=255)
                       tATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATCa                                              <  1:1038000/71‑1 (MQ=255)
                       tATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATCa                                              <  1:3282476/71‑1 (MQ=255)
                       tATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATCa                                              <  1:3567472/71‑1 (MQ=255)
                           tcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAg                                             >  1:1977900/1‑68 (MQ=255)
                           tcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGccc                                          >  1:2362009/1‑71 (MQ=255)
                             cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGt                                                                      >  1:3845927/1‑41 (MQ=255)
                             cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg                                         >  1:3255098/1‑70 (MQ=255)
                              aTGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGta                                       <  1:1617961/71‑1 (MQ=255)
                                 aaCCGTGAACTTGAAAATCTGGAAACGCACATTCCGtt                                                                     >  1:3527176/1‑38 (MQ=255)
                                   ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   <  1:3829040/54‑1 (MQ=255)
                                   ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   <  1:1668835/54‑1 (MQ=255)
                                   ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                                                   <  1:1385745/54‑1 (MQ=255)
                                     gTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGt                                >  1:363043/1‑71 (MQ=255)
                                        aaCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                                                            <  1:3058461/40‑1 (MQ=255)
                                        aaCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                                                         >  1:3599059/1‑43 (MQ=255)
                                        aaCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                                                         >  1:1281431/1‑43 (MQ=255)
                                            tGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTatat                                    >  1:2595615/1‑60 (MQ=255)
                                             gAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata                                     >  1:976105/1‑58 (MQ=255)
                                             gAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTGTGTTTgg                        >  1:2974054/1‑71 (MQ=255)
                                             gAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGttt                          <  1:2697179/69‑1 (MQ=255)
                                             gAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTGTgg                        >  1:1463936/1‑71 (MQ=255)
                                              aaaaTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTgg                        <  1:284905/70‑1 (MQ=255)
                                                aaTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata                                     <  1:3355382/55‑1 (MQ=255)
                                                aaTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTTGTCTGttt                          <  1:3759364/66‑1 (MQ=255)
                                                 aTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACg                    <  1:1329221/71‑1 (MQ=255)
                                                    tGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACgg                   <  1:3188432/69‑1 (MQ=255)
                                                    tGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTc                 >  1:1818741/1‑71 (MQ=255)
                                                      gAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCt                <  1:878161/70‑1 (MQ=255)
                                                      gAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCt                <  1:2066994/70‑1 (MQ=255)
                                                          cGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTggggg           <  1:1091028/71‑1 (MQ=255)
                                                          cGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTggggg           <  1:2364160/71‑1 (MQ=255)
                                                            cacaTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGAt         <  1:2075966/71‑1 (MQ=255)
                                                               aTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGttt                          <  1:2212787/51‑1 (MQ=255)
                                                               aTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGttt                          <  1:814828/51‑1 (MQ=255)
                                                                ttCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCt                <  1:2996272/60‑1 (MQ=255)
                                                                ttCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATg     <  1:3678976/71‑1 (MQ=255)
                                                                ttCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATg     <  1:26010/71‑1 (MQ=255)
                                                                   cGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCAc  >  1:2451975/1‑71 (MQ=255)
                                                                    |                                                                     
GAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCAC  >  W3110S.gb/775890‑776027

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: