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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,257,3760TC26.9% 80.5 / 36.8 52F7L (TTT→CTT) yecPpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base C (5/9);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCT  >  minE/1257307‑1257441
                                                                     |                                                                 
ccAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACt                                                                   <  1:1413934/70‑1 (MQ=255)
 cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                 <  1:997033/71‑1 (MQ=255)
 cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                 <  1:681098/71‑1 (MQ=255)
  aGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACctt                                                                 >  1:223717/1‑70 (MQ=255)
     gCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                              >  1:336635/1‑70 (MQ=255)
      cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATtc                                                            <  1:1531343/71‑1 (MQ=255)
      cTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATAGACTTTGGTAACCTTTAtt                                                             <  1:1880113/70‑1 (MQ=255)
         tAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTctct                                                         >  1:1956749/1‑71 (MQ=255)
          aaCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTctct                                                         <  1:1247256/70‑1 (MQ=255)
           aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGa                                                       <  1:451966/71‑1 (MQ=255)
           aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGa                                                       <  1:1689455/71‑1 (MQ=255)
             tttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGa                                                       >  1:609898/1‑69 (MQ=255)
              ttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTg                                                    >  1:1191614/1‑71 (MQ=255)
                 gTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGcc                                                  <  1:1726217/70‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatca                                          <  1:770188/70‑1 (MQ=255)
                         gtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatca                                          <  1:1960892/70‑1 (MQ=255)
                             cATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCttt                                     <  1:1316987/71‑1 (MQ=255)
                               ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATtct                                                           <  1:822858/47‑1 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTa                                                               >  1:1453519/1‑41 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTgaca                                 >  1:1937168/1‑71 (MQ=255)
                                 aaaaGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTcaca                                 >  1:294509/1‑71 (MQ=255)
                                  aaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                              >  1:435537/1‑41 (MQ=255)
                                  aaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTAt                                                              >  1:1309312/1‑41 (MQ=255)
                                  aaaGCAGAGGACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTca                                   >  1:1825870/1‑68 (MQ=255)
                                     gCAGAGGACGCTGCATGATCGACTTTGGTAACttt                                                                 <  1:1335592/35‑1 (MQ=255)
                                          ggACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                        <  1:163534/71‑1 (MQ=255)
                                           gACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                        <  1:1144161/70‑1 (MQ=255)
                                           gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTg                                                        <  1:195771/38‑1 (MQ=255)
                                           gACGCTGCATGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGa                        >  1:1059900/1‑70 (MQ=255)
                                              gCTGCATGATCGACTTTTGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACg                    <  1:1985465/71‑1 (MQ=255)
                                              gCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACg                    <  1:695214/71‑1 (MQ=255)
                                              gCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACg                    <  1:473453/71‑1 (MQ=255)
                                                tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaa                                                <  1:1448476/41‑1 (MQ=255)
                                                tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCaa                                                <  1:316927/41‑1 (MQ=255)
                                                  cATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGc                >  1:1331755/1‑71 (MQ=255)
                                                   aTGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGcc               >  1:1258800/1‑71 (MQ=255)
                                                    tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAatc                                           <  1:783556/42‑1 (MQ=255)
                                                    tGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAatc                                           >  1:1440982/1‑42 (MQ=255)
                                                    tGATCGACTTTGGTAACCTTTATTCTCTGATTGCCAAAAatc                                           >  1:341449/1‑42 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACg                    >  1:1270333/1‑63 (MQ=255)
                                                      aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACg                    >  1:400264/1‑63 (MQ=255)
                                                       tCGACTATGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCcgcg           >  1:283351/1‑71 (MQ=255)
                                                        cGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTca                                   >  1:396451/1‑46 (MQ=255)
                                                         gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa         <  1:1304099/71‑1 (MQ=255)
                                                         gACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa         <  1:263943/71‑1 (MQ=255)
                                                         gACTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa         <  1:1397412/71‑1 (MQ=255)
                                                          aCTTTGGTAACCTTTATTCTCTGATTGCCAAAAATCATCttt                                     <  1:1767500/42‑1 (MQ=255)
                                                            tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcac                                  <  1:1564928/43‑1 (MQ=255)
                                                            tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcac                                  <  1:1087574/43‑1 (MQ=255)
                                                            tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTcac                                  <  1:1370577/43‑1 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCacac                                >  1:555171/1‑43 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCacac                                >  1:117504/1‑43 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa         >  1:2094162/1‑66 (MQ=255)
                                                              tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCa         >  1:2091597/1‑66 (MQ=255)
                                                              tGGTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGAt      <  1:31816/69‑1 (MQ=255)
                                                                gTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt  <  1:1051596/71‑1 (MQ=255)
                                                                gTAACCTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCt  <  1:522446/71‑1 (MQ=255)
                                                                     |                                                                 
CCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCT  >  minE/1257307‑1257441

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: