Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,876,440 T→C 6.4% K328K (AAA→AAG eno ← enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,876,4400TC 6.4% 423.1 / 13.7 141K328K (AAA→AAGenoenolase
Reads supporting (aligned to +/- strand):  ref base T (74/58);  new base C (7/2);  total (81/60)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 156 total aligned reads displayed.

TCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGATTTTGTCGCCCAGAACTT  >  minE/1876372‑1876507
                                                                    |                                                                   
  ggTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTc                                                                  <  1:480855/70‑1 (MQ=255)
  ggTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTc                                                                  <  1:1583553/70‑1 (MQ=255)
      agagAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTc                                                                  <  1:4624269/66‑1 (MQ=255)
        agaACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAg                                                                >  1:2001127/1‑66 (MQ=255)
        agaACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAg                                                                >  1:4981690/1‑66 (MQ=255)
        agaACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCt                                                           >  1:2468298/1‑71 (MQ=255)
          aaCCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAg                                                                >  1:4603730/1‑64 (MQ=255)
            ccGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:1092509/1‑71 (MQ=255)
            ccGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:4944100/1‑71 (MQ=255)
            ccGATCTGGTTGAATTTGATAAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:2914062/1‑71 (MQ=255)
             cGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgtg                                                      >  1:558299/1‑71 (MQ=255)
               aTCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCttt                                                                   <  1:3816641/56‑1 (MQ=255)
               aTCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCttt                                                                   <  1:2629520/56‑1 (MQ=255)
               aTCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAgg                                                               <  1:5117033/60‑1 (MQ=255)
               aTCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAgg                                                               <  1:3779344/60‑1 (MQ=255)
                 cTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCa                                                                 >  1:3728913/1‑56 (MQ=255)
                    gTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTa                                               >  1:3280688/1‑71 (MQ=255)
                    gTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATAACTTCTTTCAGGATCTTg                                                         >  1:4029378/1‑61 (MQ=255)
                       gAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:3063040/70‑1 (MQ=255)
                       gAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:2848929/70‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTc                                                                  <  1:1867709/46‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:3598872/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:3414944/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:1409174/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:1727563/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:40257/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:1796853/67‑1 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:2097176/1‑69 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:5027547/1‑69 (MQ=255)
                          tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:4871698/1‑69 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTAc                                              >  1:3862310/1‑64 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:3706650/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:1150991/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:4491081/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:1686253/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:2695796/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:1075451/1‑67 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                                       >  1:4372961/1‑71 (MQ=255)
                            tGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAg                                        <  1:1302349/70‑1 (MQ=255)
                             gATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGa                                            <  1:4309621/65‑1 (MQ=255)
                             gATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGa                                            <  1:3742546/65‑1 (MQ=255)
                               tCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcg                                    <  1:840978/71‑1 (MQ=255)
                               tCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcg                                    <  1:353336/71‑1 (MQ=255)
                                cAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:4845378/1‑51 (MQ=255)
                                cAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:3222409/1‑51 (MQ=255)
                                cAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGgt                                                       >  1:2928097/1‑51 (MQ=255)
                                cAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   >  1:2014484/1‑71 (MQ=255)
                                cAGGATGGAGTTAGCGATACCTTTTTAGATACCTTCTTTCAGGATCTTGgt                                                       >  1:2605396/1‑51 (MQ=255)
                                  ggATCGAGTTGGCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAGgtcgtc                                  <  1:2755873/70‑1 (MQ=255)
                                   gATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAAAAGGTCGTCac                                <  1:4469081/71‑1 (MQ=255)
                                    aTGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:1798510/57‑1 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATc                                                            >  1:1166466/1‑41 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATc                                                            >  1:422539/1‑41 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATc                                                            >  1:648017/1‑41 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATc                                                            >  1:4928165/1‑41 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCacca                              <  1:4322975/71‑1 (MQ=255)
                                     tGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCacca                              <  1:4911037/71‑1 (MQ=255)
                                      ggAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtc                                     <  1:4488300/63‑1 (MQ=255)
                                      ggAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtc                                     <  1:1833948/63‑1 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTg                                                         >  1:1959907/1‑41 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   >  1:5025373/1‑63 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   >  1:3847053/1‑63 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   >  1:3989041/1‑63 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgtcg                                    >  1:4222619/1‑62 (MQ=255)
                                        aGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGgacgt                                   >  1:968472/1‑63 (MQ=255)
                                         gTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:1411472/1‑69 (MQ=255)
                                          ttAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg                         <  1:1934896/71‑1 (MQ=255)
                                          ttAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg                         <  1:4148953/71‑1 (MQ=255)
                                          ttAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg                         <  1:3896445/71‑1 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:3174283/2‑55 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:18579/2‑55 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:2686370/2‑55 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:5188816/2‑55 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:3215041/2‑55 (MQ=255)
                                            ggCGATGCCCTTCTGGATGCCTTCCTTCAGGATCTTGGTGTTGGTGACGAACAgg                                       >  1:900608/2‑55 (MQ=255)
                                            aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                                             <  1:1145920/49‑1 (MQ=255)
                                            aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGa                                            <  1:381944/50‑1 (MQ=255)
                                             gCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGGCACCAACCAg                         <  1:1321393/68‑1 (MQ=255)
                                                   ccTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGAtttt                <  1:2857281/71‑1 (MQ=255)
                                                    cTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCa                                 >  1:443854/1‑53 (MQ=255)
                                                    cTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCa                                 >  1:650838/1‑53 (MQ=255)
                                                    cTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCa                                 >  1:1587478/1‑53 (MQ=255)
                                                    cTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCa                                 >  1:473569/1‑53 (MQ=255)
                                                         tCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:2585999/1‑38 (MQ=255)
                                                         tCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:794740/1‑38 (MQ=255)
                                                         tCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:3521564/1‑38 (MQ=255)
                                                         tCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGaa                                           >  1:3701464/1‑38 (MQ=255)
                                                         tCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGa                                            >  1:488313/1‑37 (MQ=255)
                                                               ccTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGATTTTGTCGCCCAGaa     <  1:2216245/70‑1 (MQ=255)
                                                               ccTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGATTTTGTCGCCCAGaa     <  1:2437000/70‑1 (MQ=255)
                                                               ccTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGATTTTGTCGCCCAGaa     <  1:1303250/70‑1 (MQ=255)
                                                                cttcttTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   <  1:3104870/39‑1 (MQ=255)
                                                                cttcttTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   <  1:3460659/39‑1 (MQ=255)
                                                                cttcttTCAGGATCTTGGTGTTGGTTACGAACAGgtcgt                                   <  1:3593267/39‑1 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:2006134/1‑43 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:286959/1‑43 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:3109298/1‑43 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:348236/1‑43 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:931023/1‑43 (MQ=255)
                                                                   cttTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCaccaac                            >  1:1129601/1‑43 (MQ=255)
                                                                    |                                                                   
TCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAGCTGGATTTTGTCGCCCAGAACTT  >  minE/1876372‑1876507

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: