Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 488,534 T→C 23.9% F925F (TTT→TTC kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb488,5340TC23.9% 122.4 / 41.8 67F925F (TTT→TTCkefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base C (12/4);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 8.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCGTA  >  W3110S.gb/488469‑488601
                                                                 |                                                                   
tCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggttttgg                                                                   >  1:3279887/1‑68 (MQ=255)
        cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                          >  1:1987701/1‑69 (MQ=255)
        cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                          <  1:1930894/69‑1 (MQ=255)
        cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                          <  1:2972838/69‑1 (MQ=255)
                gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                                   >  1:1421931/1‑68 (MQ=255)
                gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                                   >  1:2571023/1‑68 (MQ=255)
                  gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                >  1:2362886/1‑69 (MQ=255)
                  gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                >  1:1220671/1‑69 (MQ=255)
                  gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                >  1:3327213/1‑69 (MQ=255)
                  gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                >  1:1428916/1‑69 (MQ=255)
                     aaaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                <  1:1962653/66‑1 (MQ=255)
                     aaaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                                <  1:169969/66‑1 (MQ=255)
                      aaCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACt                                            >  1:3358038/1‑69 (MQ=255)
                       aCTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGaaa                                                        <  1:2826014/56‑1 (MQ=255)
                        cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                                      <  1:1003560/57‑1 (MQ=255)
                        cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                                      <  1:2573322/57‑1 (MQ=255)
                        cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                                      <  1:79218/57‑1 (MQ=255)
                        cTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAAtt                                                      <  1:1242613/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)
                           cAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGtc                                       >  1:1829842/1‑69 (MQ=255)
                              tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                                  <  1:1090845/39‑1 (MQ=255)
                              tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                                  <  1:2232771/39‑1 (MQ=255)
                              tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                                  <  1:2862261/39‑1 (MQ=255)
                               ggctggcCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCg                                   >  1:1925575/1‑69 (MQ=255)
                                      cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTg                              <  1:1323852/67‑1 (MQ=255)
                                      cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTg                              <  1:1956525/67‑1 (MQ=255)
                                      cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                            <  1:927947/69‑1 (MQ=255)
                                      cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                            <  1:442615/69‑1 (MQ=255)
                                      cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                            <  1:848707/69‑1 (MQ=255)
                                       gcagcaTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                            <  1:2899112/68‑1 (MQ=255)
                                           caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc                       <  1:1249509/69‑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)
                                            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:3060973‑M1/3‑26 (MQ=255)
                                            gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                                 >  1:2450555‑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:2176695/1‑68 (MQ=255)
                                               aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt                    >  1:2328173/1‑68 (MQ=255)
                                               aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt                    >  1:2597303/1‑68 (MQ=255)
                                               aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt                    >  1:3046144/1‑68 (MQ=255)
                                                tCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                                 >  1:430822‑M1/1‑22 (MQ=255)
                                                tCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                                 >  1:443368‑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)
                                                        tCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTc           <  1:2536898/68‑1 (MQ=255)
                                                         tttggtttcggttgtgaaaagttttttggtatgat                                           >  1:367377‑M1/2‑13 (MQ=255)
                                                         tttggtttcggttgtgaaaagttttttggtatgat                                           >  1:642168‑M1/2‑13 (MQ=255)
                                                         cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTc                                          >  1:3220261/1‑36 (MQ=255)
                                                         cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCgg                                  <  1:3079877/44‑1 (MQ=255)
                                                         cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg         >  1:1507504/1‑69 (MQ=255)
                                                         cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg         >  1:3463949/1‑69 (MQ=255)
                                                         cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg         >  1:2654397/1‑69 (MQ=255)
                                                            ggttttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTg      >  1:643462/1‑69 (MQ=255)
                                                               tttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          >  1:678890/1‑62 (MQ=255)
                                                               tttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          >  1:674526/1‑62 (MQ=255)
                                                               tttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          >  1:3494262/1‑62 (MQ=255)
                                                               tttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCg    <  1:1677803/68‑1 (MQ=255)
                                                               tttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCGt   <  1:386027/69‑1 (MQ=255)
                                                                ttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          >  1:2499566/1‑61 (MQ=255)
                                                                ttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          <  1:1378511/61‑1 (MQ=255)
                                                                ttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTcc          <  1:2686265/61‑1 (MQ=255)
                                                                ttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCGt   <  1:2556633/68‑1 (MQ=255)
                                                                ttggtttACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCGTa  >  1:1289741/1‑69 (MQ=255)
                                                                 |                                                                   
TCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCGGTGCGTA  >  W3110S.gb/488469‑488601

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: