Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,257,376 T→C 54.1% F7L (TTT→CTT)  yecP → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,257,3760TC54.1% 19.9 / 101.9 74F7L (TTT→CTT) yecPpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (14/20);  new base C (25/15);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

TCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACT  >  minE/1257306‑1257445
                                                                      |                                                                     
tCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACt                                                                       >  1:1202377/1‑71 (MQ=255)
tCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACt                                                                       >  1:2536861/1‑71 (MQ=255)
   aGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACcttt                                                                    <  1:3600196/71‑1 (MQ=255)
   aGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACcttt                                                                    <  1:2238648/71‑1 (MQ=255)
       cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtc                                                                >  1:3671387/1‑71 (MQ=255)
       cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtc                                                                >  1:1168638/1‑71 (MQ=255)
         ttAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtctc                                                              >  1:744925/1‑71 (MQ=255)
          tAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtct                                                               <  1:2603548/69‑1 (MQ=255)
          tAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtct                                                               <  1:2598433/69‑1 (MQ=255)
          tAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtct                                                               <  1:1633317/69‑1 (MQ=255)
           aaCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTctct                                                             <  1:2583366/70‑1 (MQ=255)
            aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGa                                                           <  1:1039540/71‑1 (MQ=255)
            aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGa                                                           >  1:4060543/1‑71 (MQ=255)
            aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGa                                                           >  1:716716/1‑71 (MQ=255)
               ttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTg                                                        >  1:2607080/1‑71 (MQ=255)
               ttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTg                                                        >  1:3585499/1‑71 (MQ=255)
                 ggTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGc                                                       >  1:3740382/1‑70 (MQ=255)
                 ggTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGc                                                       >  1:3353000/1‑70 (MQ=255)
                        tggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAat                                                >  1:1630079/1‑70 (MQ=255)
                         tgtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatca                                              >  1:1434989/2‑71 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTTGTAACCTTTAtt                                                                 <  1:654178/52‑1 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                                  >  1:3720466/1‑51 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                                  >  1:3351825/1‑51 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatta                                              >  1:880334/1‑69 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTAtt                                                                 <  1:3908972/52‑1 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTAtt                                                                 <  1:254244/52‑1 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTAtt                                                                 <  1:1261536/52‑1 (MQ=255)
                         ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTAtt                                                                 <  1:3342053/52‑1 (MQ=255)
                             gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCtt                                          >  1:3429216/1‑71 (MQ=255)
                             gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCtt                                          >  1:578620/1‑71 (MQ=255)
                              cATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCttt                                         <  1:1870959/71‑1 (MQ=255)
                                ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTctct                                                             <  1:1197695/49‑1 (MQ=255)
                                ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTctct                                                             <  1:1095295/49‑1 (MQ=255)
                                  aaaaGCAGAGGATGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                     >  1:3534894/1‑71 (MQ=255)
                                  aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGa                                                           >  1:1105967/1‑49 (MQ=255)
                                  aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatcatc                                            >  1:2728450/1‑64 (MQ=255)
                                  aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcact                                     >  1:2717238/1‑70 (MQ=255)
                                  aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                     >  1:1561596/1‑71 (MQ=255)
                                  aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                     >  1:750950/1‑71 (MQ=255)
                                    aaGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAatc                                               >  1:1547454/1‑59 (MQ=255)
                                      gCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTa                                                                   >  1:668629/1‑37 (MQ=255)
                                      gCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTa                                                                   >  1:3262207/1‑37 (MQ=255)
                                          agGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATAATCTTTCACCTGGCTCGa                            <  1:2287760/71‑1 (MQ=255)
                                           ggACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                            <  1:26888/71‑1 (MQ=255)
                                            gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                            >  1:577953/1‑70 (MQ=255)
                                            gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                            >  1:3547378/1‑70 (MQ=255)
                                            gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                            >  1:2004623/1‑70 (MQ=255)
                                             aCGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaaa                          >  1:3436205/1‑71 (MQ=255)
                                             aCGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaaa                          >  1:913113/1‑71 (MQ=255)
                                             aCGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaa                           >  1:2917145/1‑70 (MQ=255)
                                              cGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAAc                         <  1:518500/71‑1 (MQ=255)
                                                 tGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaaa                          >  1:3326435/1‑67 (MQ=255)
                                                 tGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCt                      >  1:1341767/1‑71 (MQ=255)
                                                 tGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCt                      >  1:2510383/1‑71 (MQ=255)
                                                   cATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCt                      <  1:1040815/69‑1 (MQ=255)
                                                   cATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGc                    >  1:2631216/1‑71 (MQ=255)
                                                    aTGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCtt                                          <  1:397891/48‑1 (MQ=255)
                                                     tGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaaa                          >  1:323844/1‑63 (MQ=255)
                                                     tGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGaaa                          >  1:2708821/1‑63 (MQ=255)
                                                          gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa             <  1:863111/71‑1 (MQ=255)
                                                          gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa             <  1:2149486/71‑1 (MQ=255)
                                                             tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                         <  1:3739777/40‑1 (MQ=255)
                                                             tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                         <  1:1971200/40‑1 (MQ=255)
                                                             tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa             <  1:1111704/68‑1 (MQ=255)
                                                             tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa             <  1:350423/68‑1 (MQ=255)
                                                             tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa             <  1:2540413/68‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                            <  1:2006508/49‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt      <  1:4103386/71‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt      <  1:3016333/71‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt      <  1:222187/71‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt      <  1:3089098/71‑1 (MQ=255)
                                                                 gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt      <  1:80495/71‑1 (MQ=255)
                                                                  tAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcg                 <  1:395609/59‑1 (MQ=255)
                                                                  tAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcg                 <  1:3456812/59‑1 (MQ=255)
                                                                  tAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTa     <  1:3880701/71‑1 (MQ=255)
                                                                  tAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGc       >  1:405844/1‑69 (MQ=255)
                                                                    aCTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAAc   >  1:190212/1‑71 (MQ=255)
                                                                     cTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACt  <  1:2615238/71‑1 (MQ=255)
                                                                      |                                                                     
TCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACT  >  minE/1257306‑1257445

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: