Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,414,454 G→C 16.4% V246L (GTC→CTC)  acrE → cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,414,4540GC16.4% 161.4 / 22.2 73V246L (GTC→CTC) acrEcytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (36/25);  new base C (6/6);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAA  >  W3110S.gb/3414388‑3414522
                                                                  |                                                                    
aGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCa                                                                   <  1:2729248/70‑1 (MQ=255)
 ggCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATg                                                                 >  1:3519754/1‑71 (MQ=255)
  gCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATg                                                                 <  1:2962497/70‑1 (MQ=255)
      aaGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGaa                                                              <  1:3333571/69‑1 (MQ=255)
       aGCAATCCGTAGAGCAATGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAc                                                           >  1:72467/1‑71 (MQ=255)
       aGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAc                                                           <  1:3199247/71‑1 (MQ=255)
          aaTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACgg                                                         <  1:2819700/70‑1 (MQ=255)
            tCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTc                                                                    <  1:1559460/57‑1 (MQ=255)
            tCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCa                                                                   <  1:725952/58‑1 (MQ=255)
              cGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTc                                                       >  1:818496/1‑68 (MQ=255)
              cGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTc                                                       >  1:1428814/1‑68 (MQ=255)
                  gagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCt                                                 <  1:2552198/70‑1 (MQ=255)
                   agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACg                                                          >  1:2254177/1‑60 (MQ=255)
                     cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATcc                                             <  1:123167/71‑1 (MQ=255)
                          aaaTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaa                                                     <  1:876456/58‑1 (MQ=255)
                           aaTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGaa                                       <  1:2317114/71‑1 (MQ=255)
                            aTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTgg                                                                      <  1:578012/39‑1 (MQ=255)
                             tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGt                                                                     >  1:455240/1‑39 (MQ=255)
                             tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGt                                                                     >  1:2358821/1‑39 (MQ=255)
                             tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGt                                                                     >  1:317476/1‑39 (MQ=255)
                             tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATc                                              >  1:3618233/1‑62 (MQ=255)
                              ttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAg                                     <  1:237479/70‑1 (MQ=255)
                               tGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTAt                                               >  1:227382/1‑59 (MQ=255)
                               tGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAgg                                    >  1:3509575/1‑70 (MQ=255)
                                  aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACg                                >  1:3386603/1‑71 (MQ=255)
                                    aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATg                                                                 >  1:2119395/1‑36 (MQ=255)
                                    aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAc                                                           >  1:2299702/1‑42 (MQ=255)
                                    aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAc                                                           >  1:1238510/1‑42 (MQ=255)
                                     aGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAc                                                           >  1:156383/1‑41 (MQ=255)
                                     aGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTg                             >  1:664464/1‑71 (MQ=255)
                                     aGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTg                             >  1:1194845/1‑71 (MQ=255)
                                      ggAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaa                                                     <  1:337088/46‑1 (MQ=255)
                                      ggAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTg                             >  1:2869091/1‑70 (MQ=255)
                                      ggAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGc                            <  1:2413119/71‑1 (MQ=255)
                                       gAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGaa                                       <  1:705086/59‑1 (MQ=255)
                                        cagagcgatgaaaacgtttcagtTTGCTCATGGAAAACgg                                                         >  1:2742474‑M2/24‑40 (MQ=255)
                                        cagagcgatgaaaacgtttcagtTTGCTCATGGAAAACGGt                                                        >  1:1162304‑M2/24‑41 (MQ=255)
                                         aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCaa                          <  1:1579364/70‑1 (MQ=255)
                                         aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCaa                          >  1:228905/1‑70 (MQ=255)
                                         aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCaa                          >  1:2822961/1‑70 (MQ=255)
                                          gagcgatgaaaacgtttcagtTTGCTCATGGAAAAc                                                           <  1:3196893‑M2/15‑1 (MQ=255)
                                          gagcgatgaaaacgtttcagtTTGCTCATGGAAAAc                                                           >  1:2121871‑M2/22‑36 (MQ=255)
                                          gagcgatgaaaacgtttcagtTTGCTCATGGAAAACg                                                          <  1:3017944‑M2/16‑1 (MQ=255)
                                          gagcgatgaaaacgtttcagtTTGCTCATGGAAAACGGt                                                        <  1:2610075‑M2/18‑1 (MQ=255)
                                          gagcgatgaaaacgtttcagtTTGCTCATGGAAAACGGt                                                        <  1:526245‑M2/18‑1 (MQ=255)
                                          aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaaa                                                    >  1:1988910/1‑43 (MQ=255)
                                          aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAtt                        >  1:3603305/1‑71 (MQ=255)
                                           agcgatgaaaacgtttcagtTTGCTCATGGAAAACgg                                                         >  1:3011527‑M2/21‑37 (MQ=255)
                                           agcgatgaaaacgtttcagtTTGCTCATGGAAAACgg                                                         <  1:1658513‑M2/17‑1 (MQ=255)
                                            gcgatgaaaacgtttcagtTTGCTCATGGAAAACgg                                                         >  1:258721‑M2/20‑36 (MQ=255)
                                            gcgatgaaaacgtttcagtTTGCTCATGGAAAACg                                                          <  1:2562813‑M2/16‑1 (MQ=255)
                                            cGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCaaa                                                    >  1:2125079/1‑41 (MQ=255)
                                            cGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAg                                     >  1:920766/1‑56 (MQ=255)
                                            cGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtct                      >  1:435959/1‑71 (MQ=255)
                                            cGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtct                      >  1:894329/1‑71 (MQ=255)
                                             cgatgaaaacgtttcagtTTGCTCATGGAAAACGGt                                                        >  1:1173006‑M2/19‑36 (MQ=255)
                                              ccaccaGCAACGTAGAGTTTGTCATGGAAAACGGTCAAACCTATCCCCTGAAAgg                                    <  1:127004/55‑1 (MQ=255)
                                              ccaccaGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                    <  1:3159877/71‑1 (MQ=255)
                                                 ccaGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGATGCAATTCTCCGAt                 >  1:1445869/1‑71 (MQ=255)
                                                  caGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATccc                                            >  1:2896589/1‑43 (MQ=255)
                                                  caGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGa                                        <  1:2235564/47‑1 (MQ=255)
                                                    gCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTa                                  >  1:1984759/1‑51 (MQ=255)
                                                    gCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgt               >  1:1419298/1‑70 (MQ=255)
                                                    gCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgt               >  1:2817602/1‑70 (MQ=255)
                                                     cAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgtg              <  1:2781750/70‑1 (MQ=255)
                                                      aaCGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGAc            <  1:3151206/71‑1 (MQ=255)
                                                         gTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGaaa                                      <  1:2143882/42‑1 (MQ=255)
                                                         gTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGt         <  1:745044/71‑1 (MQ=255)
                                                          tAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGtt        <  1:1576100/71‑1 (MQ=255)
                                                          tAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGtt        <  1:2442498/71‑1 (MQ=255)
                                                          tAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGtt        >  1:1995872/1‑71 (MQ=255)
                                                            gagTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                  >  1:14305/1‑59 (MQ=255)
                                                            gagTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtga      >  1:939919/1‑71 (MQ=255)
                                                            gagTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtga      >  1:2767639/1‑71 (MQ=255)
                                                                 ggTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgatga   >  1:3474707/1‑69 (MQ=255)
                                                                 ggTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGATGCAATTCTCCGATGTGACCGTTGATGaa  <  1:280554/70‑1 (MQ=255)
                                                                  |                                                                    
AGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAA  >  W3110S.gb/3414388‑3414522

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: