Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 488,526 C→T 23.2% L923F (CTT→TTT)  kefA → fused mechanosensitive channel proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb488,5260CT23.2% 137.7 / 38.9 69L923F (CTT→TTT) kefAfused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (12/4);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GTGTTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCG  >  W3110S.gb/488460‑488594
                                                                  |                                                                    
gtgtTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCtt                                                                    <  1:1410315/69‑1 (MQ=255)
   ttCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggt                                                                 >  1:865522/1‑69 (MQ=255)
         tCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggttttgg                                                            >  1:3279887/1‑68 (MQ=255)
                 cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                   <  1:1930894/69‑1 (MQ=255)
                 cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                   <  1:2972838/69‑1 (MQ=255)
                 cGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGa                                                   >  1:1987701/1‑69 (MQ=255)
                         gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                            >  1:2571023/1‑68 (MQ=255)
                         gggATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTc                                            >  1:1421931/1‑68 (MQ=255)
                           gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                         >  1:1428916/1‑69 (MQ=255)
                           gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                         >  1:3327213/1‑69 (MQ=255)
                           gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                         >  1:1220671/1‑69 (MQ=255)
                           gATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGt                                         >  1:2362886/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: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)
                                    cAGTGGCTGGCCGCAGCATTATCCGTAGGTCttggt                                                                 >  1:3292166/1‑36 (MQ=255)
                                    cAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGtc                                >  1:1829842/1‑69 (MQ=255)
                                       tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                           <  1:2862261/39‑1 (MQ=255)
                                       tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                           <  1:1090845/39‑1 (MQ=255)
                                       tggctggcCGCAGCATTATCCGTAGGTCttggttttggt                                                           <  1:2232771/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:442615/69‑1 (MQ=255)
                                               cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                     <  1:848707/69‑1 (MQ=255)
                                               cGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGAt                     <  1:927947/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:1249509/69‑1 (MQ=255)
                                                    caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc                <  1:2806199/69‑1 (MQ=255)
                                                    caTTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTc                <  1:2024464/69‑1 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:2988575‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          <  1:3242849‑M1/40‑17 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:2450555‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:797325‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:3060973‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:379962‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:1171766‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:870414‑M1/3‑26 (MQ=255)
                                                     gatatCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttgg                                          >  1:897869‑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:927043‑M1/1‑24 (MQ=255)
                                                       tATCCGTTGGTTTTGGTTTCGGTTgtgaaaagttttttggt                                         >  1:1678224‑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:2328173/1‑68 (MQ=255)
                                                        aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt             >  1:3046144/1‑68 (MQ=255)
                                                        aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt             >  1:2853629/1‑68 (MQ=255)
                                                        aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt             >  1:2597303/1‑68 (MQ=255)
                                                        aTCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTAt             >  1:2176695/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)
                                                                  cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTc                                   >  1:3220261/1‑36 (MQ=255)
                                                                  cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCgg                           <  1:3079877/44‑1 (MQ=255)
                                                                  cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg  >  1:3463949/1‑69 (MQ=255)
                                                                  cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg  >  1:1507504/1‑69 (MQ=255)
                                                                  cTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCg  >  1:2654397/1‑69 (MQ=255)
                                                                  |                                                                    
GTGTTCGGATCGCTGGGCGTCTCTTGGGATAAACTCCAGTGGCTGGCCGCAGCATTATCCGTAGGTCTTGGTTTTGGTTTACAAGAAATTTTCGGTAACTTCGTCTCCGGTTTGATCATTCTATTCGAACGTCCG  >  W3110S.gb/488460‑488594

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: