Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,220,477 G→A 70.6% V102I (GTT→ATT)  yebS → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,220,4770GA70.6% 55.3 / 26.0 61V102I (GTT→ATT) yebSconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base A (16/27);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 9.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCT  >  minE/1220410‑1220537
                                                                   |                                                            
tATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGt                                                             <  1:801167/69‑1 (MQ=255)
         tGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAgg                                                    <  1:2022300/69‑1 (MQ=255)
         tGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAgg                                                    <  1:1491987/69‑1 (MQ=255)
         tGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAgg                                                    <  1:1941552/69‑1 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCAt                                                 >  1:792140/1‑68 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:2817242/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:313846/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1106701/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:3387810/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:2012313/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:257309/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1947029/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1884062/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:665112/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:2995030/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:2991287/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1345037/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1252712/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:1217557/1‑69 (MQ=255)
             tGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACATTATGCAAGGCATc                                                >  1:963208/1‑69 (MQ=255)
              ggggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATc                                                >  1:504505/1‑68 (MQ=255)
              ggggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCt                                               >  1:1927114/1‑69 (MQ=255)
              ggggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCt                                               >  1:1139868/1‑69 (MQ=255)
              ggggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCt                                               >  1:2277823/1‑69 (MQ=255)
                    aaCCGCTGTTGCATATCTGGCTTTTAGGCATTCGTATTGACGCCAACGTTATGCAAgg                                                    <  1:2319071/58‑1 (MQ=255)
                                 tatCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAgg                             >  1:2250560/1‑68 (MQ=255)
                                              ggCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:273418/69‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGTCCAATCAGGGCGATGCGATAACg               <  1:1133654/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCTATAACg               <  1:440185/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:3187654/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:3035842/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:3232877/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:705675/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:87958/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2951186/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2943803/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2871833/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:964013/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:251673/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2636687/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1604655/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1606499/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1616875/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1797509/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2734089/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1999826/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2034918/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2289293/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1191554/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2423105/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:1135916/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:2605203/68‑1 (MQ=255)
                                               gCATTCGTATTGACGCCAACATTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACg               <  1:262030/68‑1 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:898195/1‑69 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:1598953/1‑69 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:1964197/1‑69 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:736726/1‑69 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:2076277/1‑69 (MQ=255)
                                                    cGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCg         >  1:2420307/1‑69 (MQ=255)
                                                            cGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCt  <  1:2588462/68‑1 (MQ=255)
                                                            cGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGAGGCGATAACGGGGTAGATGTTCt  <  1:96828/68‑1 (MQ=255)
                                                                   |                                                            
TATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCT  >  minE/1220410‑1220537

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: