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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,257,3760TC13.0% 200.8 / 19.8 77F7L (TTT→CTT) yecPpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base C (2/8);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTG  >  minE/1257307‑1257446
                                                                     |                                                                      
ccAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACtt                                                                       >  1:925010/1‑71 (MQ=255)
ccAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACtt                                                                       >  1:816527/1‑71 (MQ=255)
 cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                      <  1:2735571/71‑1 (MQ=255)
 cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                      <  1:3580764/71‑1 (MQ=255)
 cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                      <  1:1462675/71‑1 (MQ=255)
  aGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                      >  1:3723678/1‑70 (MQ=255)
   gTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTa                                                                    <  1:4110750/71‑1 (MQ=255)
    tGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACtt                                                                       >  1:1754632/1‑67 (MQ=255)
    tGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACtt                                                                       >  1:1048857/1‑67 (MQ=255)
    tGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                                   <  1:2819500/71‑1 (MQ=255)
    tGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                                   <  1:260184/71‑1 (MQ=255)
    tGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                                   <  1:1780242/71‑1 (MQ=255)
     gCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                      <  1:2348311/67‑1 (MQ=255)
      cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtc                                                                 >  1:4106718/1‑71 (MQ=255)
      cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtc                                                                 >  1:2257996/1‑71 (MQ=255)
          aaCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTctct                                                              <  1:288201/70‑1 (MQ=255)
           aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTctct                                                              <  1:3549258/69‑1 (MQ=255)
              ttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTctct                                                              <  1:120175/66‑1 (MQ=255)
                ggTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGc                                                        >  1:2470649/1‑70 (MQ=255)
                    cACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCaaa                                                    >  1:838876/1‑70 (MQ=255)
                     aCTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGCTAACTTTTa                                                                    <  1:3044617/53‑1 (MQ=255)
                        ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtc                                                                 <  1:1236698/53‑1 (MQ=255)
                        ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTctct                                                              <  1:1750163/56‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtct                                                                <  1:3112947/53‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtct                                                                <  1:1536368/53‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtct                                                                <  1:2228837/53‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAatc                                                <  1:3626723/69‑1 (MQ=255)
                           ggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTctct                                                              >  1:1145980/1‑53 (MQ=255)
                           ggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCt                                            >  1:361394/1‑71 (MQ=255)
                            gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACctt                                                                      <  1:3293534/44‑1 (MQ=255)
                            gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACctt                                                                      <  1:383499/44‑1 (MQ=255)
                            gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACctt                                                                      <  1:1756648/44‑1 (MQ=255)
                            gCATTAAAAGCAGAGGACGCTCCATGATCGACTTTTGTAACctt                                                                      <  1:214632/44‑1 (MQ=38)
                             cATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:1726009/71‑1 (MQ=255)
                             cATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:3945364/71‑1 (MQ=255)
                               ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:1863548/69‑1 (MQ=255)
                               ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:2485697/69‑1 (MQ=255)
                                tAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTca                                        >  1:3953559/1‑70 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACtt                                                                       >  1:319264/1‑38 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:3527549/1‑71 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:1414576/1‑71 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:1461827/1‑71 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:2407644/1‑71 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:313938/1‑71 (MQ=255)
                                        gagGACGCTGCTTGATCGACTTTGGTATGTTTTAt                                                                   <  1:906330/35‑1 (MQ=25)
                                          ggACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGa                                                            <  1:1585102/40‑1 (MQ=255)
                                          ggACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                             <  1:3643891/71‑1 (MQ=255)
                                           gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                             <  1:937862/70‑1 (MQ=255)
                                                tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          >  1:2219788/1‑52 (MQ=255)
                                                tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          >  1:862631/1‑52 (MQ=255)
                                                tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          >  1:1723460/1‑52 (MQ=255)
                                                 gCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTca                                        >  1:2503168/1‑53 (MQ=255)
                                                 gCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTca                                        >  1:3557578/1‑53 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:1798883/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:1226324/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:1305596/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:3652085/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:1792782/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:3873886/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaaa                                                    <  1:2278843/38‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatcat                                              <  1:410757/44‑1 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCt                                            >  1:3128820/1‑44 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCt                                            >  1:2354633/1‑44 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:738914/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:423041/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:2970144/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:1212705/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:1907653/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgc                 >  1:2702683/1‑71 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcg                  >  1:874116/1‑70 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcg                  >  1:838539/1‑70 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTATGATTGCCAAAAatcat                                              >  1:51322/1‑42 (MQ=255)
                                                       tCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatcat                                              <  1:447498/41‑1 (MQ=255)
                                                       tCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:870736/45‑1 (MQ=255)
                                                       tCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:1772089/45‑1 (MQ=255)
                                                         gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa              <  1:342719/71‑1 (MQ=255)
                                                         gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa              <  1:2040894/71‑1 (MQ=255)
                                                            tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:1720377/40‑1 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:3060878/1‑42 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:4092567/1‑42 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                      >  1:868536/1‑42 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcac                                       >  1:3170983/1‑41 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGACAAAAATCATCTTTcaca                                      >  1:875191/1‑42 (MQ=255)
                                                                gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt       <  1:1462119/71‑1 (MQ=255)
                                                                gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt       <  1:1113790/71‑1 (MQ=255)
                                                                 tAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTa      <  1:3857218/71‑1 (MQ=255)
                                                                 gaaCTTTTATTCTCTGATTGCCAAAAATCATCttt                                          <  1:2682795/34‑1 (MQ=255)
                                                                   aCTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAAc    <  1:2773737/71‑1 (MQ=255)
                                                                   aCTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAAc    <  1:3464387/71‑1 (MQ=255)
                                                                   aCTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAAc    <  1:3661786/71‑1 (MQ=255)
                                                                    cTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACt   <  1:3558670/71‑1 (MQ=255)
                                                                     ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa              <  1:683280/59‑1 (MQ=255)
                                                                     ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa              <  1:2138883/59‑1 (MQ=255)
                                                                     ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTg  >  1:2336597/1‑71 (MQ=255)
                                                                     |                                                                      
CCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTG  >  minE/1257307‑1257446

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: