Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 984,696 A→C 13.1% G151G (GGT→GGG)  ynfC ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE984,6960AC13.1% 183.1 / 18.3 91G151G (GGT→GGG) ynfChypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base C (6/6);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

Only 100 of 107 total aligned reads displayed.

CCGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCATTAT  >  minE/984629‑984763
                                                                   |                                                                   
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:1194690/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:1394449/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:1357764/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:1625734/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:2881195/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:1161890/69‑1 (MQ=255)
ccGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:472712/69‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:958790/64‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCTGCGGATACCCCTGa                                                               <  1:3195049/69‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCAGCGGATACCCCTGa                                                               <  1:589788/69‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCAGCGGATACCCCTGa                                                               <  1:2848222/69‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCAGCGGATACCCCTGa                                                               <  1:1615745/69‑1 (MQ=255)
     gATGGCGTGGCGCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCAGCGGATACCCCTGa                                                               <  1:1529648/69‑1 (MQ=255)
       tggcgtggcgCTGACAGATAATGTTTTATCGTTACTTTTCGTGGTTTTACCCAGCGGATACCCCTGa                                                               <  1:2473206/67‑1 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCTTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:737454/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:1015647/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:2440137/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:1871906/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:871173/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:1202216/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:1446910/1‑67 (MQ=255)
           gtggcgCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACACTGatcat                                                           >  1:2180456/1‑67 (MQ=255)
           gtggagCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcat                                                           >  1:1104295/1‑67 (MQ=255)
              gcgcTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATa                                                       >  1:1515680/1‑68 (MQ=255)
              gcgcTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATa                                                       >  1:757185/1‑68 (MQ=255)
              gcgcTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATa                                                       >  1:1827949/1‑68 (MQ=255)
                    aCAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAAc                                                                    <  1:2616058/49‑1 (MQ=255)
                        aTAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCTTAGCGATATTCCa                                            <  1:2336195/69‑1 (MQ=255)
                        aTAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCa                                            <  1:540805/69‑1 (MQ=255)
                        aTAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCa                                            <  1:310355/69‑1 (MQ=255)
                          aaTGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCAt                                           >  1:420981/1‑68 (MQ=255)
                          aaTGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATc                                          >  1:841793/1‑69 (MQ=255)
                          aaTGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATc                                          >  1:1368573/1‑69 (MQ=255)
                          aaTGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATc                                          >  1:472367/1‑69 (MQ=255)
                                       tACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTg                                        <  1:2600374/58‑1 (MQ=255)
                                       tACTTTTCGTGGTTTTACCCAGCGGATAACCCT‑ATCATCATAGCGATATTCCATCTg                                        <  1:1216111/57‑1 (MQ=255)
                                        tcTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgc                             >  1:2683822/2‑68 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGatcatc                                                          >  1:2013346/1‑39 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTg                                        <  1:1003463/57‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:3257547/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:1146647/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:3155472/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:3749/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2725031/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:437110/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:306073/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:1145452/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2976957/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2925271/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:55480/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:620494/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:708507/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:871157/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:1232132/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            >  1:1449318/1‑69 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:134231/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:1635093/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2095824/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2421554/69‑1 (MQ=255)
                                        aCTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2488242/69‑1 (MQ=255)
                                         cTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:3391161/68‑1 (MQ=255)
                                                 ggTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCttt                   <  1:1829203/69‑1 (MQ=255)
                                                   ttttACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:2036625/58‑1 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTTGctgctg                        >  1:742700/1‑60 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTTGCTGCTGGCTTTAATc               >  1:2576749/1‑69 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTTGCTGCTGGCTTTAATc               >  1:2614351/1‑69 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTTGCTGCTGGCTTTAATc               >  1:1233103/1‑69 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTTGCTGCTGGCTTTAATc               >  1:811220/1‑69 (MQ=255)
                                                     ttACCCAGCGGATACCCCTGATCATCATAGCGATATTCCATCTGCATCTGTTAGCTGCTGGCTTTAATc               >  1:925126/1‑69 (MQ=255)
                                                     ttACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATc               >  1:2551725/1‑69 (MQ=255)
                                                     ttACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATc               >  1:1460510/1‑69 (MQ=255)
                                                     ttACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATc               >  1:3408781/1‑69 (MQ=255)
                                                            gCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGc        <  1:1770526/69‑1 (MQ=255)
                                                            gCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGc        <  1:712090/69‑1 (MQ=255)
                                                            gCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGc        <  1:2468940/69‑1 (MQ=255)
                                                            gCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGc        <  1:2574039/69‑1 (MQ=255)
                                                            gCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGc        <  1:1023309/69‑1 (MQ=255)
                                                               gATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTtgct                            <  1:878280/46‑1 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:469758/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:2788087/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:1752592/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:1105723/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:601825/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:1948319/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:3135668/1‑68 (MQ=255)
                                                                 tAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCAtt    >  1:948854/1‑68 (MQ=255)
                                                                  aaCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCATTat  <  1:2588164/69‑1 (MQ=255)
                                                                  aaCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCATTat  <  1:2158244/69‑1 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:1472931/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:1531420/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:1754158/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:7346/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:654203/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:3259700/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:2079636/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:226418/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:2718812/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:1211097/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:1182183/1‑52 (MQ=255)
                                                                   aCCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTa                  >  1:2502510/1‑52 (MQ=255)
                                                                   |                                                                   
CCGTTGATGGCGTGGCGCTGACAGATAATGTTTTGTCGTTACTTTTCGTGGTTTTACCCAGCGGATAACCCTGATCATCATAGCGATATTCCATCTGCATTTGTTTGCTGCTGGCTTTAATCACGAAGCCATTAT  >  minE/984629‑984763

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: ‑