Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 488,522 A→T 23.2% V921V (GTA→GTT kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb488,5220AT23.2% 131.9 / 42.6 69V921V (GTA→GTTkefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base T (12/4);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 8.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

ACGGTGTTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAT  >  W3110S.gb/488457‑488583
                                                                 |                                                             
aCGGTGTTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGt                                                            >  1:1651090/1‑69 (MQ=255)
   gtgtTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCtt                                                         <  1:1410315/69‑1 (MQ=255)
      ttCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggt                                                      >  1:865522/1‑69 (MQ=255)
            tCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggttttgg                                                 >  1:3279887/1‑68 (MQ=255)
                    cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                        >  1:1987701/1‑69 (MQ=255)
                    cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                        <  1:2972838/69‑1 (MQ=255)
                    cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                        <  1:1930894/69‑1 (MQ=255)
                            gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                 >  1:2571023/1‑68 (MQ=255)
                            gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                 >  1:1421931/1‑68 (MQ=255)
                              gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              >  1:3327213/1‑69 (MQ=255)
                              gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              >  1:2362886/1‑69 (MQ=255)
                              gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              >  1:1220671/1‑69 (MQ=255)
                              gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              >  1:1428916/1‑69 (MQ=255)
                                 aaaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              <  1:169969/66‑1 (MQ=255)
                                 aaaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              <  1:1962653/66‑1 (MQ=255)
                                  aaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACt                          >  1:3358038/1‑69 (MQ=255)
                                   aCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGaaa                                      <  1:2826014/56‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                    <  1:79218/57‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                    <  1:1242613/57‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                    <  1:2573322/57‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                    <  1:1003560/57‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCTGTAACTTc                        <  1:2820154/69‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              <  1:2535051/63‑1 (MQ=255)
                                    cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTc                        <  1:981728/69‑1 (MQ=255)
                                       cAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggt                                                      >  1:3292166/1‑36 (MQ=255)
                                       cAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGtc                     >  1:1829842/1‑69 (MQ=255)
                                          tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                <  1:1090845/39‑1 (MQ=255)
                                          tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                <  1:2862261/39‑1 (MQ=255)
                                          tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                <  1:2232771/39‑1 (MQ=255)
                                           ggctggcCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCg                 >  1:1925575/1‑69 (MQ=255)
                                                  cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTg            <  1:1956525/67‑1 (MQ=255)
                                                  cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTg            <  1:1323852/67‑1 (MQ=255)
                                                  cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt          <  1:927947/69‑1 (MQ=255)
                                                  cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt          <  1:848707/69‑1 (MQ=255)
                                                  cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt          <  1:442615/69‑1 (MQ=255)
                                                   gcagcaTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt          <  1:2899112/68‑1 (MQ=255)
                                                       caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:228123/69‑1 (MQ=255)
                                                       caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:1917903/69‑1 (MQ=255)
                                                       caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:2024464/69‑1 (MQ=255)
                                                       caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:2806199/69‑1 (MQ=255)
                                                       caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:1249509/69‑1 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:797325‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:870414‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:379962‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:897869‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:1171766‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               <  1:3242849‑M1/40‑17 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:2988575‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:2450555‑M1/3‑26 (MQ=255)
                                                        gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:3060973‑M1/3‑26 (MQ=255)
                                                        aTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     <  1:95966/68‑1 (MQ=255)
                                                         ttATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCt    >  1:1848944/1‑68 (MQ=255)
                                                         ttATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTa   >  1:2641476/1‑69 (MQ=255)
                                                          tATCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttggt                              >  1:1678224‑M1/1‑24 (MQ=255)
                                                          tATCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttggt                              >  1:927043‑M1/1‑24 (MQ=255)
                                                          tATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTa   >  1:2015671/1‑68 (MQ=255)
                                                           aTCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttggt                              <  1:3281144‑M1/40‑18 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                              >  1:840139/1‑40 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTa   >  1:2258341/1‑67 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt  >  1:2853629/1‑68 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt  >  1:2328173/1‑68 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt  >  1:2176695/1‑68 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt  >  1:2597303/1‑68 (MQ=255)
                                                           aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt  >  1:3046144/1‑68 (MQ=255)
                                                            tCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:443368‑M1/1‑22 (MQ=255)
                                                            tCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                               >  1:430822‑M1/1‑22 (MQ=255)
                                                            tCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc     >  1:1964438/1‑64 (MQ=255)
                                                             ccGTTGGTTTTGGTTTCGGTTgtgaaaagtttttt                                 <  1:3031426‑M1/35‑15 (MQ=255)
                                                              cGTTGGTTTTGGTTTCGGTTgtgaaaagttttttggt                              <  1:2979144‑M1/37‑18 (MQ=255)
                                                                 |                                                             
ACGGTGTTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAT  >  W3110S.gb/488457‑488583

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: