Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,780,767 T→C 41.3% I65I (ATT→ATC ypjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,780,7670TC41.3% 12.1 / 31.6 46I65I (ATT→ATCypjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base C (15/4);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCCCGGTCTGCTGCAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAC  >  minE/1780704‑1780833
                                                               |                                                                  
tGTTCCCGGTCTGCTGCAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATCgcgc                                                                <  1:567339/68‑1 (MQ=255)
tGTTCCCGGTCTGCTGCAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATCgcgc                                                                <  1:2190790/68‑1 (MQ=255)
tGTTCCCGGGCTGCTGCAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATCgcgc                                                                <  1:496903/68‑1 (MQ=255)
               gcAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCa                                                  <  1:2145208/67‑1 (MQ=255)
               gcAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCa                                                  <  1:2052175/67‑1 (MQ=255)
               gcAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCa                                                  <  1:1438243/67‑1 (MQ=255)
                  aaaaaaCGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:3015904/69‑1 (MQ=255)
                  aaaaaaCGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:2145767/69‑1 (MQ=255)
                  aaaaaaCGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:1825790/69‑1 (MQ=255)
                  aaaaaaCGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:713844/69‑1 (MQ=255)
                   aaaaaCGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:2195848/68‑1 (MQ=255)
                        cggcggCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGc                                             <  1:743170/63‑1 (MQ=255)
                        cggcggCTGGCGGCGCATGGCTATTATTTCTGCGGTCATCgcgc                                                                <  1:2367841/44‑1 (MQ=255)
                                           gCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCCCGCAATCGCTCTGGAAGCCCGCATCCTGCCCGAc                    >  1:2731316/1‑69 (MQ=255)
                                           gCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGa                     >  1:3360703/1‑68 (MQ=255)
                                           gCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGAc                    >  1:2160081/1‑69 (MQ=255)
                                              attattTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGt                 <  1:2801466/69‑1 (MQ=255)
                                              attattTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGt                 <  1:202706/69‑1 (MQ=255)
                                              attattTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGt                 <  1:1928514/69‑1 (MQ=255)
                                                 attTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGa               >  1:1578553/1‑68 (MQ=255)
                                                 attTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGAt              >  1:851681/1‑69 (MQ=255)
                                                 attTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGAt              >  1:1355433/1‑69 (MQ=255)
                                                 attTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGAt              >  1:1191100/1‑69 (MQ=255)
                                                 attTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGAt              >  1:2638163/1‑69 (MQ=255)
                                                    tCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAg            >  1:630616/1‑68 (MQ=255)
                                                          gTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAg            >  1:2148472/1‑62 (MQ=255)
                                                          gTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaa     >  1:797961/1‑69 (MQ=255)
                                                          gTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaa     >  1:2608765/1‑69 (MQ=255)
                                                           tCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    <  1:1294418/69‑1 (MQ=255)
                                                           tCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaa     <  1:1547548/68‑1 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCGGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:85066/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:996315/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2983295/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:992136/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:864341/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2407096/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:1291655/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2363003/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2283536/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2185059/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2182705/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2140617/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:2055594/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:1879576/1‑69 (MQ=255)
                                                           tCATCGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaa    >  1:1419126/1‑69 (MQ=255)
                                                            cATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACaaaa   >  1:991370/1‑69 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCt                           >  1:386946/1‑42 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:1238808/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:250873/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:940014/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:875934/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:1534847/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:3374997/1‑67 (MQ=255)
                                                               tGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAc  >  1:2452378/1‑67 (MQ=255)
                                                               |                                                                  
TGTTCCCGGTCTGCTGCAAAAAAACGGCGGCTGGCGGCGCATGGCTATTATTTCTGCGGTCATTGCGCTGGTCTGCCACGCAATCGCTCTGGAAGCCCGCATCCTGCCCGACGGTGATAGCGGACAAAAC  >  minE/1780704‑1780833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: