Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,305,362 T→C 15.1% T233A (ACG→GCG)  yeeE ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,305,3620TC15.1% 125.6 / 13.7 53T233A (ACG→GCG) yeeEpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (3/5);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCG  >  minE/1305297‑1305427
                                                                 |                                                                 
cAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTgggg                                                              <  1:1873821/71‑1 (MQ=255)
cAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTgggg                                                              <  1:3487806/71‑1 (MQ=255)
cAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTgggg                                                              <  1:2178872/71‑1 (MQ=255)
cAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTgggg                                                              <  1:2828549/71‑1 (MQ=255)
cAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTggg                                                               >  1:2078642/1‑70 (MQ=255)
    aaaaCGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGt                                                          <  1:503956/71‑1 (MQ=255)
    aaaaCGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGt                                                          <  1:3355133/71‑1 (MQ=255)
         gCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTc                                                     <  1:4767852/71‑1 (MQ=255)
         gCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTc                                                     <  1:2862365/71‑1 (MQ=255)
          ccccAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGAt                                                       <  1:1238515/68‑1 (MQ=255)
            ccAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTa                                                  <  1:186714/71‑1 (MQ=255)
               gTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTgggg                                                              <  1:1481100/56‑1 (MQ=255)
                   aTGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:4427513/71‑1 (MQ=255)
                    tGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:3487624/70‑1 (MQ=255)
                      tATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCTAACAt                                        >  1:1193333/1‑71 (MQ=255)
                      tATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACAt                                        >  1:1975007/1‑71 (MQ=255)
                      tATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACAt                                        >  1:581691/1‑71 (MQ=255)
                      tATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACAt                                        >  1:2915586/1‑71 (MQ=255)
                        tttCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGc                                      >  1:781541/1‑71 (MQ=255)
                                     gcgACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           >  1:4992606/1‑53 (MQ=255)
                                     gcgACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCtt                          >  1:4326255/1‑70 (MQ=255)
                                        aCCAGAAATTGCAGAATATTGGCCGCTGGTGAAGTGATTCCTAACCCaaa                                           <  1:4908383/50‑1 (MQ=38)
                                        aCCAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCaaa                                           <  1:2888344/50‑1 (MQ=255)
                                        aCCAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCaaa                                           <  1:4648367/50‑1 (MQ=255)
                                        aCCAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCaaa                                           <  1:4564930/50‑1 (MQ=255)
                                        aCCAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCaaa                                           <  1:453073/50‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCTCTGGGGAAGTGATTCCTAACCCAAACa                                         >  1:434920/1‑50 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:2981440/48‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:25822/48‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:3542989/48‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaaa                                           <  1:2593280/48‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCa                    <  1:1100004/71‑1 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCAAACa                                         >  1:4360821/1‑50 (MQ=255)
                                          cAGAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCAAACa                                         >  1:398463/1‑50 (MQ=255)
                                            gAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAgg                  <  1:3970176/71‑1 (MQ=255)
                                            gAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAgg                  <  1:4509993/71‑1 (MQ=255)
                                             aaaTTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAggg                 <  1:447658/71‑1 (MQ=255)
                                               aTTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGc               <  1:3407456/71‑1 (MQ=255)
                                                  gCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCaa                                            >  1:4141191/1‑39 (MQ=255)
                                                  gCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAg            >  1:2686056/1‑71 (MQ=255)
                                                  gCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAg            >  1:558951/1‑71 (MQ=255)
                                                  gCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAg            >  1:1543503/1‑71 (MQ=255)
                                                     gAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGc          >  1:3862897/1‑70 (MQ=255)
                                                     gAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGc          >  1:1443054/1‑70 (MQ=255)
                                                      aaTATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCt         >  1:3310053/1‑70 (MQ=255)
                                                      aaTATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCt         >  1:4073393/1‑70 (MQ=255)
                                                          ttGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAg    <  1:3108894/70‑1 (MQ=255)
                                                           tGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGc                                      >  1:2913286/1‑36 (MQ=255)
                                                           tGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGtt                              >  1:2607566/1‑44 (MQ=255)
                                                           tGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTaaaa     >  1:4636626/1‑69 (MQ=255)
                                                            ggCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAggg                 <  1:3848132/56‑1 (MQ=255)
                                                            ggCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAgcg  >  1:2391326/1‑71 (MQ=255)
                                                              ccGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCa                    <  1:3704180/51‑1 (MQ=255)
                                                                 ttGGGGAAGTGATTCCTAACCCAAACATGCGTCCgg                                <  1:215481/36‑1 (MQ=255)
                                                                 ttGGGGAAGTGATTCCTAACCCAAACATGCGTCCgg                                <  1:3378493/36‑1 (MQ=255)
                                                                 ttGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTaa       <  1:3568354/61‑1 (MQ=255)
                                                                 |                                                                 
CAGGAAAACGCCCCAGTTAATGTATTTCATGTCACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCG  >  minE/1305297‑1305427

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