Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,220,863 T→C 22.8% L230L (CTT→CTC yebS → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,220,8630TC22.8% 58.6 / 12.4 48L230L (CTT→CTCyebSconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base C (6/5);  total (12/36)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGTCGCTGCC  >  minE/1220808‑1220918
                                                       |                                                       
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGAGGAGAAACTTCGTGTCTGTCTTg                                            <  1:83185/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:3333961/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:2589416/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:2639415/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:3059848/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:3415897/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:1457502/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:512064/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:1306460/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:1279082/69‑1 (MQ=255)
aaCTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:846284/69‑1 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:3003689/59‑1 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGc                                       >  1:2391073/1‑64 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAttt                                  <  1:2048924/69‑1 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAttt                                  <  1:2823988/69‑1 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAttt                                  <  1:2255395/69‑1 (MQ=255)
          gCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTCTCTTGGGTGc                                       >  1:403186/1‑64 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:2813784/69‑1 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:3087109/69‑1 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:2498015/69‑1 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:2505897/69‑1 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:107209/69‑1 (MQ=255)
           cGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:2171003/69‑1 (MQ=255)
                        aGCGCCCCGCTACGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGc                           <  1:79393/62‑1 (MQ=255)
                        aGCGCCCCGCTACGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGc                           <  1:633530/62‑1 (MQ=255)
                        aGCGCCCCGCTACGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGc                           <  1:331318/62‑1 (MQ=255)
                        aGCGCCCCGCTACGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGc                           <  1:1513441/62‑1 (MQ=255)
                              ccGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGc              <  1:70785/69‑1 (MQ=255)
                              ccGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGc              <  1:2237380/69‑1 (MQ=255)
                              ccGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGc              <  1:697773/69‑1 (MQ=255)
                              ccGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGc              <  1:324349/69‑1 (MQ=255)
                                gCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTg                                            <  1:3133200/37‑1 (MQ=255)
                                  tACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTgg          >  1:1974911/1‑69 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGt         <  1:1071902/69‑1 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGt         <  1:517663/69‑1 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGt         <  1:2474917/69‑1 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTAc                               >  1:1069114/1‑47 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCTGa                   >  1:3413768/1‑59 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCTGa                   >  1:3179657/1‑59 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCTGa                   >  1:2157792/1‑59 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCTGa                   >  1:769708/1‑59 (MQ=255)
                                   aCGCGTAGGGACGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCTGa                   >  1:2733046/1‑59 (MQ=255)
                                    cgcgTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGt         <  1:2710328/68‑1 (MQ=255)
                                       gTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCAtttt                                 <  1:146314/41‑1 (MQ=255)
                                          gggACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGTCGCTGcc  >  1:1329143/1‑69 (MQ=255)
                                          gggACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGTCGCTGcc  >  1:651424/1‑69 (MQ=255)
                                          gggACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGTCGCTGcc  >  1:773898/1‑69 (MQ=255)
                                              cGAGAAACTCCGTGTCTGTCTTGGGTGCCATTTTACCGGc                           <  1:1864004/40‑1 (MQ=255)
                                                       |                                                       
AACTGTGGGAGCGATTTTATCCGCAGCGCCCCGCTACGCGTAGGGACGAGAAACTTCGTGTCTGTCTTGGGTGCCATTTTACCGGCTATCCAGATCAGCGTGGTCGCTGCC  >  minE/1220808‑1220918

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: