Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,864,319 C→T 9.3% intergenic (+66/+82) ygcP → / ← ygcQ predicted anti‑terminator regulatory protein/predicted flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,864,3190CT 9.3% 308.4 / 14.5 107intergenic (+66/+82)ygcP/ygcQpredicted anti‑terminator regulatory protein/predicted flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (59/38);  new base T (4/6);  total (63/44)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCAT  >  minE/1864251‑1864389
                                                                    |                                                                      
tGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAgcgc                                                                      >  1:2331284/1‑71 (MQ=255)
tGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAgcgc                                                                      >  1:2134346/1‑71 (MQ=255)
  aCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCAGGCCACATCCGACATTCAGCGCCt                                                                    <  1:2484716/71‑1 (MQ=255)
  aCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGcc                                                                     <  1:1052084/70‑1 (MQ=255)
  aCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCt                                                                    >  1:2909533/1‑71 (MQ=255)
   cGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCt                                                                    >  1:773609/1‑70 (MQ=255)
         gcctgatttagcgtgattttgtaggtcggataaggcgttaatgccgcatccgacatcaaCGCCTGATGCGa                                                             >  1:1712698‑M2/60‑71 (MQ=255)
         aaTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGa                                                             >  1:1134952/1‑71 (MQ=255)
          cctgatttagcgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGAc                                                            >  1:279019‑M2/59‑71 (MQ=255)
          aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGAc                                                            <  1:3861122/71‑1 (MQ=255)
           ctgatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTg                                                                   <  1:3023827‑M2/6‑1 (MQ=255)
           ctgatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTg                                                                   <  1:573430‑M2/6‑1 (MQ=255)
           ctgatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTg                                                                   <  1:1515700‑M2/6‑1 (MQ=255)
            tgatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTGATg                                                                >  1:2714985‑M2/57‑65 (MQ=255)
             gatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTGATGCGACGCt                                                         >  1:800747‑M2/56‑71 (MQ=255)
             gatttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCt                                                         >  1:351329‑M2/56‑71 (MQ=255)
             gTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGc                                                          <  1:2173661/70‑1 (MQ=255)
             gTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCAACATCCGACATTCAGCGCCTGATGCGACGCt                                                         >  1:3788797/1‑71 (MQ=255)
              atttagcgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCtt                                                        <  1:3225617‑M2/17‑1 (MQ=255)
              atttagcgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCt                                                         >  1:541617‑M2/55‑70 (MQ=255)
              atttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCtt                                                        >  1:1298180‑M2/55‑71 (MQ=255)
              atttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCtt                                                        >  1:2856372‑M2/55‑71 (MQ=255)
              atttagcgtgattttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCtt                                                        <  1:3521988‑M2/17‑1 (MQ=255)
               tttagcgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGAt                                                                 >  1:189275‑M2/54‑61 (MQ=255)
               tttagcgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGctta                                                       <  1:873801‑M2/18‑2 (MQ=255)
               aaCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGa                                                                  >  1:1634287/1‑60 (MQ=255)
                 ccAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAc                                                     >  1:3527311/1‑71 (MQ=255)
                  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGa                                                             >  1:3033681/1‑62 (MQ=255)
                  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcg                                                    >  1:400592/1‑71 (MQ=255)
                   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgc                                                   <  1:840684/71‑1 (MQ=255)
                    cgtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCGGATGCGACGCTTAAcgc                                                   <  1:3529163‑M2/22‑1 (MQ=255)
                    aTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgcg                                                  >  1:603614/1‑71 (MQ=255)
                    aTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgcg                                                  >  1:351447/1‑71 (MQ=255)
                     ttGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGt                                                 >  1:2226349/1‑71 (MQ=255)
                     gtgattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGt                                                 >  1:1139882‑M2/48‑71 (MQ=255)
                      tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGa                                                      >  1:789276/1‑65 (MQ=255)
                       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                >  1:2433217‑M2/46‑70 (MQ=255)
                       gattttgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTc                                                >  1:1099010‑M2/46‑70 (MQ=255)
                       gAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGa                                                                  <  1:874481/52‑1 (MQ=255)
                       gAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGa                                                             >  1:3609912/1‑57 (MQ=255)
                        attttgtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCt                                               >  1:1342290‑M2/45‑70 (MQ=255)
                        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATg                                                                <  1:1193226/53‑1 (MQ=255)
                        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCt                                               >  1:2518750/1‑70 (MQ=255)
                         aTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTa                                             >  1:3074421/1‑71 (MQ=255)
                          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                            >  1:883636/1‑71 (MQ=255)
                          tttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                            >  1:360383/1‑71 (MQ=255)
                           ttgtaggtcggataaggcgtttacgccgcatccgacatcaaTGCCTGATGCGACGCTTGCCGCGTCTTAt                                            <  1:2019613‑M2/29‑1 (MQ=255)
                            tgtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATc                                           <  1:3765823‑M2/30‑1 (MQ=255)
                            tgtaggtcggataaggcgtttacgccgcatccgacatcaac                                                                        >  1:2163208‑M2/41‑41 (MQ=255)
                            tGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTAt                                            >  1:1662009/1‑69 (MQ=255)
                             gtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCa                                          <  1:1576870‑M2/31‑1 (MQ=255)
                             gtaggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCa                                          <  1:1235923‑M2/31‑1 (MQ=255)
                             gtaggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCa                                          <  1:2287383‑M2/31‑1 (MQ=255)
                              taggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAgg                                        <  1:2784971‑M2/33‑1 (MQ=255)
                              taggtcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                         <  1:2708781‑M2/32‑1 (MQ=255)
                              gTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                        >  1:2157464/1‑71 (MQ=255)
                              gTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAg                                         <  1:2989546/70‑1 (MQ=255)
                               tttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                       <  1:3588221/71‑1 (MQ=255)
                               aggtcggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                       >  1:1841386‑M2/38‑71 (MQ=255)
                                ttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                      >  1:1392949/1‑71 (MQ=255)
                                 tGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                      >  1:2537266/1‑70 (MQ=255)
                                  tcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                    >  1:1208360‑M2/35‑71 (MQ=255)
                                  tcggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTa                                    >  1:1896200‑M2/35‑71 (MQ=255)
                                  gATTTGTAGGCCGCACGTCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCt                                     <  1:111699/70‑1 (MQ=255)
                                   cggataaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAc                                   >  1:2838742‑M2/34‑71 (MQ=255)
                                    tttGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                   <  1:2371164/70‑1 (MQ=255)
                                    ggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTacg                                  >  1:1834397‑M2/33‑70 (MQ=255)
                                    ggataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTacg                                  <  1:2167983‑M2/39‑2 (MQ=255)
                                     gataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCtacgc                                 <  1:3770279‑M2/40‑3 (MQ=255)
                                     gataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCtacgc                                 <  1:3496256‑M2/40‑3 (MQ=255)
                                      ataaggcgtttatgccgcatccgacatcaaCGCCTGATGCGAc                                                            <  1:1628412‑M2/13‑1 (MQ=255)
                                       taaggcgtttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCtacgc                                 >  1:2336942‑M2/30‑67 (MQ=255)
                                       taaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCtt                                                        >  1:1947796‑M2/30‑46 (MQ=255)
                                       taaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtacgc                                 >  1:2597252‑M2/30‑67 (MQ=255)
                                       gTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAc                                                     >  1:1623181/1‑49 (MQ=255)
                                       gTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCa                                          >  1:1590822/1‑60 (MQ=255)
                                       gTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCa                                          >  1:3770073/1‑60 (MQ=255)
                                        aaggcgtttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtacgc                                 >  1:2146674‑M2/29‑66 (MQ=255)
                                          ggCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTccc                            <  1:1629385/71‑1 (MQ=255)
                                             cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAg                                         >  1:732232/1‑55 (MQ=255)
                                             cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                   >  1:3056204/1‑61 (MQ=255)
                                             cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCa                          >  1:2858827/1‑70 (MQ=255)
                                              tttatgccgcatccgacatcaaCGCCTGATGCGACGCTTAAcgcg                                                  <  1:243830‑M2/23‑1 (MQ=255)
                                              tttacgccgcatccgacatcaaTGCCTGATGCGACGCTTGCCGCGTCTTATCa                                          <  1:617188‑M2/31‑1 (MQ=255)
                                              tttacgccgcatccgacatcaaCGCCTGATGCGACGCTTGCcgcg                                                  <  1:1168714‑M2/23‑1 (MQ=255)
                                              gCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAgg                        >  1:2418114/1‑71 (MQ=255)
                                               cACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAGTCCCCAGGc                       <  1:2034244/70‑1 (MQ=255)
                                                 cGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGc                       >  1:2341777/1‑69 (MQ=255)
                                                    cacaTCCGACATTCAGCGCCTGATGCGACGCTTGAc                                                     <  1:2593669/36‑1 (MQ=255)
                                                      caTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAgg                        >  1:3739030/1‑63 (MQ=255)
                                                      caTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAAcc                >  1:3650948/1‑71 (MQ=255)
                                                      caTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAAcc                <  1:2492146/71‑1 (MQ=255)
                                                       aTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg               <  1:679890/71‑1 (MQ=255)
                                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGt              <  1:1094792/71‑1 (MQ=255)
                                                        atccgacatcaaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCctatc                                  <  1:2678970‑M2/39‑3 (MQ=255)
                                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCg               >  1:2859123/1‑69 (MQ=255)
                                                          cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAg            <  1:1254918/71‑1 (MQ=255)
                                                             cATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGCACCGTAgg           >  1:3225905/1‑69 (MQ=255)
                                                                    cgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatcat  >  1:1793332/1‑71 (MQ=255)
                                                                    cgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatcat  >  1:2132657/1‑71 (MQ=255)
                                                                    |                                                                      
TGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCAT  >  minE/1864251‑1864389

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: