Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,790,638 C→T 13.2% intergenic (+33/‑742) pgi → / → ubiC glucosephosphate isomerase/chorismate pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,790,6380CT13.2% 148.0 / 11.1 53intergenic (+33/‑742)pgi/ubiCglucosephosphate isomerase/chorismate pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (5/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

TGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCC  >  minE/2790569‑2790706
                                                                     |                                                                    
tGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                     >  1:2505344/1‑71 (MQ=255)
tGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                     >  1:1585350/1‑71 (MQ=255)
 ggTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                     >  1:1536010/1‑70 (MQ=255)
  gTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                    <  1:2842301/70‑1 (MQ=255)
  gTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                    <  1:210972/70‑1 (MQ=255)
       aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                               <  1:2523028/70‑1 (MQ=255)
        ttAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                             <  1:3554798/71‑1 (MQ=255)
         tAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                            >  1:1344852/1‑71 (MQ=255)
          aaCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                           >  1:94173/1‑71 (MQ=255)
          aaCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                            <  1:680092/70‑1 (MQ=255)
          aaCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                            <  1:2380976/70‑1 (MQ=255)
              gCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                        >  1:1508970/1‑70 (MQ=255)
               cTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                      >  1:3767946/1‑71 (MQ=255)
               cTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                      >  1:3377727/1‑71 (MQ=255)
                  taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                    >  1:279988/1‑70 (MQ=255)
                  taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                     >  1:3302448/1‑69 (MQ=255)
                  taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCg                                                   <  1:1988246/71‑1 (MQ=255)
                    aaGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                 >  1:150484/1‑71 (MQ=255)
                     aGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                                >  1:3244876/1‑71 (MQ=255)
                      gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                  <  1:3910258/52‑1 (MQ=255)
                      gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                               >  1:1368499/1‑71 (MQ=255)
                      gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                               >  1:3815276/1‑71 (MQ=255)
                       cGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                  <  1:2047637/51‑1 (MQ=255)
                       cGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                          <  1:1822452/59‑1 (MQ=255)
                        gTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                             <  1:2609538/71‑1 (MQ=255)
                            cgcgGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                     >  1:1467790/1‑43 (MQ=255)
                                gTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                     <  1:2660733/71‑1 (MQ=255)
                                   aaTCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                 >  1:3708286/1‑40 (MQ=255)
                                         cgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGTTGCGACGCGGTc                               <  1:1892250/68‑1 (MQ=255)
                                          gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                            >  1:2580379/1‑70 (MQ=255)
                                             gATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                         >  1:447588/1‑70 (MQ=255)
                                               ttcgTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAACAGATGCCTGATGCGACGCTACCGCGTCTTAt                      <  1:382514/68‑1 (MQ=11)
                                               tatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                      <  1:3247193/71‑1 (MQ=255)
                                                 tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                     <  1:1509578/70‑1 (MQ=255)
                                                   tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                  <  1:380818/71‑1 (MQ=255)
                                                   tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                   >  1:1369066/1‑70 (MQ=255)
                                                    aGGACGGATAAGGCGTTTACGCCGCATCCGGCAACAGATGCCTGATGCGACGCTACCGCGTATTATCa                    >  1:2394589/1‑68 (MQ=14)
                                                     ggCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                  <  1:2896853/53‑1 (MQ=255)
                                                      gacGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCATGcc                >  1:1671754/3‑70 (MQ=11)
                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                >  1:1719234/1‑70 (MQ=255)
                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt               <  1:612785/71‑1 (MQ=255)
                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                             >  1:661493/1‑56 (MQ=255)
                                                       acGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCATGCCTa              >  1:3894761/2‑71 (MQ=11)
                                                       acGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCATGCCTa              <  1:1830808/70‑1 (MQ=11)
                                                        cGGATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCATGCCTa              >  1:2182873/1‑70 (MQ=11)
                                                         ggATAAGGCGTTTACGCTGCATCCGGCAACCTATGCCTGATGCGACGCTGAAGCGTCTTATCATGCCTACa            >  1:1131528/1‑71 (MQ=11)
                                                            tAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                      <  1:1557150/42‑1 (MQ=255)
                                                            tAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt         >  1:2354780/1‑71 (MQ=255)
                                                              aGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc        <  1:1910216/70‑1 (MQ=255)
                                                                 cGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa      >  1:219223/1‑69 (MQ=255)
                                                                 cGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt     >  1:1458042/1‑70 (MQ=255)
                                                                   ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc        <  1:2689770/65‑1 (MQ=255)
                                                                   ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc  >  1:3530579/1‑71 (MQ=255)
                                                                     |                                                                    
TGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCC  >  minE/2790569‑2790706

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: