Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,744 C→T 33.3% intergenic (‑281/+173) ldrA ← / ← ldrC toxic polypeptide, small/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE821,7440CT33.3% 91.0 / 50.5 63intergenic (‑281/+173)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (12/9);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAA  >  minE/821679‑821809
                                                                 |                                                                 
gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                               >  1:804117/1‑70 (MQ=40)
gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                               >  1:509117/1‑70 (MQ=40)
gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCttt                                                                >  1:87792/1‑69 (MQ=40)
gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCttt                                                                >  1:1489830/1‑69 (MQ=40)
  tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtt                                                                 <  1:198576/66‑1 (MQ=40)
  tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                            <  1:2913701/71‑1 (MQ=39)
    ctacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                               >  1:66514/1‑66 (MQ=40)
       cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAc                                                           >  1:2593556/1‑67 (MQ=40)
       cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAg                                                         <  1:2751842/69‑1 (MQ=40)
       cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                                        <  1:586455/70‑1 (MQ=40)
       cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCt                                                       >  1:3906341/1‑71 (MQ=39)
        gTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTg                                                      <  1:3762849/71‑1 (MQ=39)
        gTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTg                                                      <  1:2767886/71‑1 (MQ=39)
         tctcAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGt                                                     >  1:1724755/1‑71 (MQ=39)
             aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                            >  1:3542756/1‑60 (MQ=35)
             aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                                        <  1:3135994/64‑1 (MQ=35)
             aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                                 <  1:914898/71‑1 (MQ=34)
                  aTCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCg                                             <  1:743751/70‑1 (MQ=39)
                   tCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                                              <  1:1725448/52‑1 (MQ=35)
                    cttcacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                               <  1:2922922/50‑1 (MQ=38)
                          ttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                <  1:1271260/59‑1 (MQ=35)
                           tcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                <  1:1585193/58‑1 (MQ=35)
                           tcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                <  1:2102745/58‑1 (MQ=35)
                             aTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGc                                 <  1:2768597/71‑1 (MQ=14)
                              tCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGct                                <  1:1387584/71‑1 (MQ=14)
                                ccTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                       <  1:1129744/62‑1 (MQ=11)
                                 cTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTg                             <  1:2463129/71‑1 (MQ=17)
                                 cTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTg                             <  1:3720935/71‑1 (MQ=255)
                                  tGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAg                                    >  1:583189/1‑63 (MQ=14)
                                  tGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTg                             >  1:3138346/1‑70 (MQ=255)
                                   gAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGga                           >  1:3123666/1‑71 (MQ=17)
                                     aaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgat                          >  1:3682894/1‑70 (MQ=17)
                                     aaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgat                          >  1:240717/1‑70 (MQ=17)
                                     aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGctc                               >  1:23727/1‑65 (MQ=255)
                                     aaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatg                         <  1:3007666/71‑1 (MQ=255)
                                      aTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatga                        <  1:1700420/71‑1 (MQ=255)
                                      aTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatga                        <  1:1275452/71‑1 (MQ=255)
                                       tGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatgag                       <  1:3071787/71‑2 (MQ=18)
                                            gCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCtt                   <  1:2229396/70‑1 (MQ=255)
                                             cAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                <  1:1323280/56‑1 (MQ=255)
                                             cAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                <  1:3878576/56‑1 (MQ=255)
                                             cAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTc                  <  1:57373/70‑1 (MQ=255)
                                                tgaaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatg                         <  1:3052251/58‑1 (MQ=21)
                                                tAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGc              >  1:2248697/1‑71 (MQ=255)
                                                 gaaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTc                  >  1:1806389/2‑66 (MQ=21)
                                                  aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatga                        >  1:1164136/1‑59 (MQ=21)
                                                  aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGga                           >  1:2556154/1‑56 (MQ=21)
                                                  aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGga                           >  1:2005732/1‑56 (MQ=21)
                                                  aaGAGTGCATTCCGGTATTTCAACAGCTGTTACAGTCATTTCATGAGAGCTCTGgatga                        >  1:2361553/1‑59 (MQ=21)
                                                  aaGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGtt                                           >  1:3570292/1‑40 (MQ=40)
                                                     agTGCATTCCGGTTTTTCAACAGCTGTTACAGTCAttt                                          <  1:540275/38‑1 (MQ=11)
                                                      gTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCtt                   <  1:3690362/60‑1 (MQ=255)
                                                         cATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCGAGCTCggg         >  1:161017/1‑67 (MQ=255)
                                                          aTTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCt                    >  1:661338/1‑55 (MQ=25)
                                                          aTTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCt                    <  1:2248109/55‑1 (MQ=255)
                                                           ttCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                       <  1:410323/35‑1 (MQ=17)
                                                           ttCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                       <  1:3387662/35‑1 (MQ=17)
                                                            tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa  >  1:1973720/1‑71 (MQ=38)
                                                            tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa  >  1:2617331/1‑71 (MQ=38)
                                                            tCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCa   >  1:2305213/1‑70 (MQ=255)
                                                              cGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                >  1:2953427/1‑39 (MQ=255)
                                                              cGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                >  1:2129108/1‑39 (MQ=255)
                                                               ggCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCt             >  1:2577897/1‑57 (MQ=255)
                                                                 cTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                          >  1:3242477/1‑42 (MQ=255)
                                                                 |                                                                 
GATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAA  >  minE/821679‑821809

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: