Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 778,653 T→G 5.1% intergenic (‑359/‑161) minC ← / → ycgJ cell division inhibitor/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE778,6530TG 5.1% 226.5 / 6.3 79intergenic (‑359/‑161)minC/ycgJcell division inhibitor/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/28);  new base G (2/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TTTTTTCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAGGT  >  minE/778586‑778723
                                                                   |                                                                      
ttttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTa                                                                      <  1:525095/70‑1 (MQ=255)
ttttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTAt                                                                     >  1:4567335/1‑71 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTAt                                                                     >  1:2573946/1‑69 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTAt                                                                     >  1:112673/1‑69 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTAt                                                                     >  1:4715025/1‑69 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTAt                                                                     >  1:3674252/1‑69 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATc                                                                    <  1:337859/70‑1 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCt                                                                   >  1:4216857/1‑71 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATAACGGGTATCt                                                                   >  1:2217585/1‑71 (MQ=255)
  ttttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATAACGGGTATCt                                                                   >  1:3142079/1‑71 (MQ=255)
    ttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCtt                                                                  <  1:1136818/70‑1 (MQ=255)
    ttCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATAACGGGTATCtg                                                                  <  1:1733582/70‑2 (MQ=255)
       accacGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATAACGGGTATCtg                                                                  <  1:3166201/67‑2 (MQ=255)
            gCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGc                                                         >  1:3969170/1‑71 (MQ=255)
             ccAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGc                                                         >  1:1804013/1‑70 (MQ=255)
               aGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCca                                                       >  1:556056/1‑70 (MQ=255)
                ggCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCca                                                       >  1:3613147/1‑69 (MQ=255)
                ggCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCca                                                       >  1:4743759/1‑69 (MQ=255)
                          aCTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGctc                                            >  1:4535597/1‑70 (MQ=255)
                          aCTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGctc                                            >  1:3291955/1‑70 (MQ=255)
                          aCTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGctc                                            >  1:4013157/1‑70 (MQ=255)
                           cTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGcc                                                        >  1:3678830/1‑57 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:4691493/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:3934900/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:811176/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:1216073/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:1399211/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:518381/1‑61 (MQ=255)
                            tttCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:1872579/1‑61 (MQ=255)
                              tCCATTTTTTGCGGCATCATTAAAGAACATATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       >  1:3825653/1‑71 (MQ=255)
                              tCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       >  1:2556654/1‑71 (MQ=255)
                              tCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       >  1:2973035/1‑71 (MQ=255)
                               ccATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       >  1:3716847/1‑70 (MQ=255)
                                cATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       >  1:630364/1‑69 (MQ=255)
                                  ttttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGa                                   <  1:3173538/71‑1 (MQ=255)
                                   tttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAat                                                               >  1:2572496/1‑42 (MQ=255)
                                   tttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAat                                                               >  1:4548897/1‑42 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTaa                                                                >  1:1819419/1‑40 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTaa                                                                >  1:570561/1‑40 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTaa                                                                >  1:2313587/1‑40 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCTCCACAGAATTGCTCTATCATAGaa                                  <  1:1780487/70‑1 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGaa                                  <  1:491144/70‑1 (MQ=255)
                                    ttttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGaa                                  <  1:1991353/70‑1 (MQ=255)
                                     tttGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGaaa                                 <  1:117948/70‑1 (MQ=255)
                                            catcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaa                         <  1:3147323/71‑1 (MQ=255)
                                            catcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaa                         <  1:3463048/71‑1 (MQ=255)
                                            catcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaa                         <  1:3091298/71‑1 (MQ=255)
                                             atcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaag                        <  1:97421/71‑1 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       <  1:1002966/55‑1 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       <  1:4226830/55‑1 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCa                                       <  1:138761/55‑1 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaaga                       >  1:168371/1‑71 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaaga                       >  1:4509603/1‑71 (MQ=255)
                                              tcatTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTACTCTATCa                                       <  1:4114202/55‑1 (MQ=255)
                                                atTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGAc                      <  1:1642819/70‑1 (MQ=255)
                                                 ttAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGt                    <  1:3252349/71‑1 (MQ=255)
                                                 ttAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGt                    <  1:4736901/71‑1 (MQ=255)
                                                  tAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGaa                                                  >  1:1085756/1‑40 (MQ=255)
                                                  tAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGaa                                                  >  1:736954/1‑40 (MQ=255)
                                                  tAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGaa                                                  >  1:2557916/1‑40 (MQ=255)
                                                  tAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGa                                                   >  1:4381275/1‑39 (MQ=255)
                                                         cAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACt            >  1:4378500/1‑71 (MQ=255)
                                                         cAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACt            <  1:4703686/71‑1 (MQ=255)
                                                          aGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACtt           >  1:4621857/1‑71 (MQ=255)
                                                          aGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACtt           <  1:4079734/71‑1 (MQ=255)
                                                          aGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACtt           >  1:4304873/1‑71 (MQ=255)
                                                          aGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACt            >  1:4246049/1‑70 (MQ=255)
                                                           gATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTc          >  1:3931842/1‑71 (MQ=255)
                                                            aTACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCagaa                         <  1:2313939/55‑1 (MQ=255)
                                                              aCCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCt         >  1:4482667/1‑69 (MQ=255)
                                                               cccGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGt       >  1:1915873/1‑70 (MQ=255)
                                                               cccGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGt       >  1:1817861/1‑70 (MQ=255)
                                                                ccGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGt                    <  1:2900329/56‑1 (MQ=255)
                                                                ccGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCa     <  1:2403012/71‑1 (MQ=255)
                                                                 cGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCaga                          <  1:2898528/49‑1 (MQ=255)
                                                                 cGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCag                           <  1:1602690/48‑1 (MQ=255)
                                                                 cGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAg    <  1:362990/71‑1 (MQ=255)
                                                                 cGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAg    <  1:4749688/71‑1 (MQ=255)
                                                                 cGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAg    >  1:281326/1‑71 (MQ=255)
                                                                   ttATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAGGt  >  1:4358980/1‑71 (MQ=255)
                                                                   ttATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAGGt  >  1:897745/1‑71 (MQ=255)
                                                                   |                                                                      
TTTTTTCACCACGCCAGGCTTTCTTAACTTTCCATTTTTTGCGGCATCATTAAAGAACAGATACCCGTTATCTTAATATCCGCCACAGAATTGCTCTATCATAGAAACCGCAGAAGACGTCTGATACTTCTGTCAGGT  >  minE/778586‑778723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: