Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 160,302 T→C 31.1% G198G (GGT→GGC degP → serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE160,3020TC31.1% 65.6 / 36.8 77G198G (GGT→GGCdegPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base T (32/21);  new base C (12/12);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

AAGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCC  >  minE/160237‑160362
                                                                 |                                                            
aaGATGTCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTg                                                           >  1:699/1‑69 (MQ=255)
aaGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCt                                                            >  1:3059500/1‑68 (MQ=255)
aaGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCt                                                            >  1:2971326/1‑68 (MQ=255)
aaGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTg                                                           >  1:3362245/1‑69 (MQ=255)
aaGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTg                                                           >  1:3105240/1‑69 (MQ=255)
aaGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTg                                                           >  1:2174646/1‑69 (MQ=255)
     ggCGGATTCTGATGTACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                      >  1:3250915/1‑69 (MQ=255)
     ggCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGc                                                        >  1:2551810/1‑67 (MQ=255)
     ggCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                      >  1:1485344/1‑69 (MQ=255)
     ggCGGATTCTGATGCACTGCGCGTGGGTGATTACACAGTAGCGATTGGTAACCCGTTTGGTCTGGGCga                                                      >  1:97754/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:1119617/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:2963544/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:2763561/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:2187001/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:1843399/1‑69 (MQ=255)
           ttCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGt                                                >  1:1181201/1‑69 (MQ=255)
            tCTGCTGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTa                                               >  1:2782648/1‑69 (MQ=255)
            tCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTa                                               >  1:9450/1‑69 (MQ=255)
            tCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTa                                               >  1:3011616/1‑69 (MQ=255)
            tCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTa                                               >  1:1637256/1‑69 (MQ=255)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:516895/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:294070/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:2642991/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:3160327/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:3232421/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:2369009/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:1243867/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:699760/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:763701/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:833395/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGGCgaga                                                    <  1:196282/46‑1 (MQ=37)
                              ggTGATTACACCGTAGCAATTGGTAACCCGTTCGGCCTGGCCgaga                                                    <  1:86214/46‑1 (MQ=25)
                                    tACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTg                       >  1:2183200/1‑69 (MQ=255)
                                    tACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTg                       >  1:2363276/1‑69 (MQ=255)
                                    tACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTg                       >  1:1703304/1‑69 (MQ=255)
                                    tACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTg                       >  1:2915109/1‑69 (MQ=255)
                                     acacCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCTGGATTGTCTCTGCGCTgg                      >  1:94091/1‑69 (MQ=255)
                                     acacCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                      >  1:3034272/1‑69 (MQ=255)
                                     acacCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                      >  1:1049760/1‑69 (MQ=255)
                                     acacCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTgg                      >  1:1627781/1‑69 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:2876187/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:2245263/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:2186246/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:2179985/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:746971/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:2136866/68‑1 (MQ=255)
                                               gATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCg             <  1:1824118/68‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGTGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:364121/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:2098469/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:842190/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:2076570/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:682127/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:612967/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:577497/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:571055/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:458936/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:1560534/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:265481/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:2706649/69‑1 (MQ=255)
                                                   ggTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:1375261/69‑1 (MQ=255)
                                                    gTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTg        <  1:1361969/68‑1 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGt                 >  1:3241839/1‑56 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGaacg    >  1:3132218/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGaacg    >  1:3128204/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGaacg    >  1:972563/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGaacg    >  1:925056/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTCTGCGCTGGGACGTAGCGGCCTGaacg    >  1:2956930/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:261945/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:1635896/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:1422709/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:14399/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:2094666/1‑67 (MQ=255)
                                                       aCCCGTTCGGCCTGGGCGAGACAGTAACTTCCGGGATTGTCTATGCGCTGGGACGTAGCGGCCTGaacg    >  1:1642071/1‑67 (MQ=255)
                                                          cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGcc  >  1:483689/1‑68 (MQ=255)
                                                          cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGcc  >  1:967660/1‑68 (MQ=255)
                                                          cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGcc  >  1:1426054/1‑68 (MQ=255)
                                                          cGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGcc  >  1:2640666/1‑68 (MQ=255)
                                                                 |                                                            
AAGATGGCGGATTCTGATGCACTGCGCGTGGGTGATTACACCGTAGCGATTGGTAACCCGTTTGGTCTGGGCGAGACGGTAACTTCCGGGATTGTCTCTGCGCTGGGGCGTAGCGGCCTGAATGCC  >  minE/160237‑160362

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: