Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,396,068 G→T 20.6% M258I (ATG→ATT setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,396,0680GT20.6% 68.8 / 7.0 29M258I (ATG→ATTsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base G (19/4);  new base T (2/4);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGC  >  minE/1395999‑1396138
                                                                     |                                                                      
catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGatgatg                                                                        <  1:1046174/70‑1 (MQ=255)
catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                                       >  1:97344/1‑71 (MQ=255)
catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                                       >  1:2594023/1‑71 (MQ=255)
catcatCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgg                                                                       <  1:1777031/71‑1 (MQ=255)
 atcatcAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTTATGATggg                                                                      <  1:242835/71‑1 (MQ=255)
 atcatcAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                                      <  1:590192/71‑1 (MQ=255)
  tcatcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                                      >  1:1562378/1‑70 (MQ=255)
         cATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGTGAccgcc                                                               >  1:606880/1‑70 (MQ=255)
              cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt                                                         <  1:3109331/71‑2 (MQ=255)
              cGCTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCggt                                                         <  1:3437316/71‑2 (MQ=255)
                cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt                                                       <  1:1336594/71‑1 (MQ=255)
                cTATTTATTATTGATGAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCt                                                       <  1:2095777/71‑1 (MQ=255)
                   tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATgggg                                                                     <  1:3233829/54‑1 (MQ=255)
                   tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTgg                                                     >  1:1163031/1‑70 (MQ=255)
                     tattatCCACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccgccgc                                                             >  1:1807884/1‑60 (MQ=255)
                     tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGatgatg                                                                        <  1:1271501/49‑1 (MQ=255)
                     tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa                                                   <  1:2062447/70‑1 (MQ=255)
                     tattatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGaa                                                   <  1:1376619/70‑1 (MQ=255)
                          tCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAGATCCcg                                             <  1:1542191/71‑1 (MQ=255)
                             atgAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTAtcgctgccg                                                           >  1:2349334/3‑54 (MQ=255)
                             aCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAccg                                                                 >  1:116538/1‑48 (MQ=255)
                              cGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCTcgacga                                         <  1:506966/71‑1 (MQ=255)
                               gAACTGCATCTAACCGATAAACTGACTGGAGAAATGATTGGTATCGCTGCCGGTCTGGAAATTCCgatgat                                        >  1:1785542/1‑71 (MQ=255)
                                 aCTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGt                                      >  1:2811283/1‑71 (MQ=255)
                                    gCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGa                                   >  1:458160/1‑71 (MQ=255)
                                     cATCTTCCCGAGAAACTGGCCGGTGTGATGATggg                                                                      <  1:602895/35‑1 (MQ=255)
                                     cATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTg                                    >  1:2621546/1‑69 (MQ=255)
                                       tCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGA‑‑‑CCGCCGGGCTGGAAATcc                                               >  1:3070076/1‑53 (MQ=255)
                                        cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGAc                                                                   >  1:2547902/1‑35 (MQ=255)
                                        cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGc                               >  1:2723217/1‑71 (MQ=255)
                                               agaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATc                                                >  1:929507/1‑47 (MQ=255)
                                                gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGatat                        >  1:3039409/1‑70 (MQ=255)
                                                gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt                       >  1:2691213/1‑71 (MQ=255)
                                                gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATAtt                       >  1:2215856/1‑71 (MQ=255)
                                                     tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc                   >  1:2007992/1‑70 (MQ=255)
                                                     tGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGc                   >  1:954128/1‑70 (MQ=255)
                                                         cGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAAc              >  1:1135815/1‑71 (MQ=255)
                                                            tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa                           >  1:3357020/1‑55 (MQ=255)
                                                            tgtgATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGa                           >  1:557733/1‑55 (MQ=255)
                                                                atgatgGGGACCGCCGCCGGGCTGGAAATCCcgac                                           <  1:74870/35‑1 (MQ=255)
                                                                atgatgGGGACCGCCGCCGGGCTGGAAATACCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTggg       >  1:1859397/1‑71 (MQ=255)
                                                                 tgatgGGGGCCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTggg       >  1:2012507/1‑70 (MQ=255)
                                                                     ggggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGc  >  1:2506609/1‑71 (MQ=255)
                                                                     |                                                                      
CATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGC  >  minE/1395999‑1396138

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: