Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 775,925 C→T 11.3% N117N (AAC→AAT 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,9250CT11.3% 116.8 / 5.3 44N117N (AAC→AATtolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base C (26/13);  new base T (2/3);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ATAGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTA  >  W3110S.gb/775864‑775990
                                                             |                                                                 
aTAGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:914668/1‑71 (MQ=255)
aTAGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:3655038/1‑71 (MQ=255)
aTAGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:1130348/1‑71 (MQ=255)
  aGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTg                                                         >  1:599846/1‑70 (MQ=255)
  aGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTg                                                         >  1:1591106/1‑70 (MQ=255)
  aGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTg                                                         >  1:2733028/1‑70 (MQ=255)
  aGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTg                                                         >  1:587807/1‑70 (MQ=255)
         cgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGaaa                                                      >  1:2189071/1‑66 (MQ=255)
         cgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTg                                                 <  1:3723550/71‑1 (MQ=255)
         cgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTg                                                 >  1:757596/1‑71 (MQ=255)
         cgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTg                                                 >  1:2266147/1‑71 (MQ=255)
         cgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTg                                                 >  1:685589/1‑71 (MQ=255)
          gcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATTAAc                                                                   >  1:3252010/1‑52 (MQ=255)
           ccGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGa                                               <  1:1256650/71‑1 (MQ=255)
            cGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAAc                                                            <  1:1883843/57‑1 (MQ=255)
             ggAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAAc                                                            <  1:2224981/56‑1 (MQ=255)
              gAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:2250345/1‑57 (MQ=255)
              gAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:568252/1‑57 (MQ=255)
              gAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACtt                                                          >  1:416458/1‑57 (MQ=255)
               aaGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACg                                           <  1:2517166/71‑1 (MQ=255)
                       gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTTAACTTGAAAAt                                                    <  1:1170527/54‑1 (MQ=255)
                       gTGGAAGGGGCGTCCCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAAt                                                    <  1:3270425/54‑1 (MQ=255)
                          gAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGaaa                                             >  1:2953122/1‑58 (MQ=255)
                          gAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGaaa                                             >  1:664826/1‑58 (MQ=255)
                          gAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGtt                                <  1:525437/71‑1 (MQ=255)
                             ggggCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTcc                              >  1:1650808/1‑70 (MQ=255)
                               ggCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCg                                  <  1:1006763/64‑1 (MQ=255)
                                   tCGCGTGCTATGCGTATCTCCATGAACCGTTAACTTGAAAATCTGGa                                               <  1:1425626/47‑1 (MQ=255)
                                   tCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTAc                       >  1:3228372/1‑71 (MQ=255)
                                   tCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTAc                       >  1:764982/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:3567472/71‑1 (MQ=255)
                                                 tATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATCa         <  1:3282476/71‑1 (MQ=255)
                                                 tATTTCCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATCa         <  1:1038000/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)
                                                             |                                                                 
ATAGCCATGCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTA  >  W3110S.gb/775864‑775990

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