Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 323,233 G→T 19.0% P120P (CCG→CCT recR → gap repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE323,2330GT19.0% 254.4 / 40.0 126P120P (CCG→CCTrecRgap repair protein
Reads supporting (aligned to +/- strand):  ref base G (66/36);  new base T (9/15);  total (75/51)
Fisher's exact test for biased strand distribution p-value = 2.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 126 total aligned reads displayed.

TTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTC  >  minE/323175‑323299
                                                          |                                                                  
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGgatg                                                                >  1:2470084/1‑63 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGatat                                                            >  1:3180172/1‑67 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATc                                                           >  1:2083509/1‑68 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATc                                                           >  1:119098/1‑68 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATc                                                           >  1:1539203/1‑68 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATc                                                           >  1:1318006/1‑68 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:584216/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:1713459/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:68884/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:393277/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:2016845/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:807063/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:2027159/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:2267187/1‑69 (MQ=255)
tttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:2880139/1‑69 (MQ=255)
tttCAGGTCGTTAATTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATc                                                           >  1:3337741/1‑68 (MQ=255)
tttAAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:1542589/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGTGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:1945273/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:1996476/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:526904/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:2940792/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:2704059/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:919366/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:1639592/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:711206/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:1423353/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:732161/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCgg                                                         >  1:843765/1‑69 (MQ=255)
 ttCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCg                                                          >  1:917374/1‑68 (MQ=255)
      gTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:2943272/69‑1 (MQ=255)
      gTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:2368907/69‑1 (MQ=255)
      gTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:2135878/69‑1 (MQ=255)
         gTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATc                                                 <  1:12019/69‑1 (MQ=255)
           tATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATc                                                 <  1:2212769/67‑1 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:2273442/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:1142922/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:3039922/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:2989909/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:2448508/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:654803/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATCGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:2864375/1‑67 (MQ=255)
                 gtgtTGATGGGACACCTGTCACCGCTGGACGGCATAGGTCCTGATGATATCGGGCTTGATCGTCTgg                                           >  1:1861318/1‑67 (MQ=255)
                                cTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCaga                           >  1:2844391/1‑68 (MQ=255)
                                cTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTAGATCGTCTGGAACAGCGTCTGGCaga                           >  1:1504451/1‑68 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCagag                          >  1:557279/1‑68 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCagag                          >  1:1874912/1‑68 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCagag                          >  1:608847/1‑68 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCagag                          >  1:874609/1‑68 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCaga                           >  1:1000322/1‑67 (MQ=255)
                                 tGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATAGTCTGGAACAGCGTCTGGCagag                          >  1:2851671/1‑68 (MQ=255)
                                 agTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAAc                                        >  1:1024680/2‑54 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:3373335/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:3337441/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:2819052/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:2523431/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:3140026/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:287087/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:41631/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:1786439/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                     >  1:2140161/1‑69 (MQ=255)
                                     aCCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaa                      >  1:2803073/1‑68 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGTCTTg                                                    <  1:897075/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:1628824/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:3415993/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:1821372/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:2997139/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:2637821/35‑1 (MQ=255)
                                        gCTGGACGGCATCGGTCCGGATGATATCGGGCTTg                                                    <  1:1058262/35‑1 (MQ=255)
                                              cGGCATCGGTCCGGATGATATCGGGGTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg             <  1:3430062/68‑1 (MQ=255)
                                              cGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTg             <  1:3088197/68‑1 (MQ=255)
                                                  aTCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg        >  1:691655/1‑69 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:2039607/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:467513/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:1304181/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:1281352/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:1146266/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:2790810/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:2510457/69‑1 (MQ=255)
                                                   tCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGa       <  1:3055292/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:321467/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:2682980/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:740435/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:1646971/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:624634/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:3065489/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:20084/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:2667187/69‑1 (MQ=255)
                                                     ggTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGGCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:2208683/69‑1 (MQ=255)
                                                     ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTg        <  1:2042097/66‑1 (MQ=255)
                                                     ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     >  1:1974518/1‑69 (MQ=255)
                                                     ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     >  1:1602607/1‑69 (MQ=255)
                                                     ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     >  1:1286847/1‑69 (MQ=255)
                                                     ggTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     >  1:1118715/1‑69 (MQ=255)
                                                      gTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:3055105/68‑1 (MQ=255)
                                                      gTCCTGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATc     <  1:2114541/68‑1 (MQ=255)
                                                        ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc  <  1:1885545/69‑1 (MQ=255)
                                                        ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc  <  1:2158109/69‑1 (MQ=255)
                                                        ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc  <  1:2337412/69‑1 (MQ=255)
                                                        ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc  <  1:835960/69‑1 (MQ=255)
                                                        ccGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc  <  1:2568519/69‑1 (MQ=255)
                                                          |                                                                  
TTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTC  >  minE/323175‑323299

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: