Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,766 G→A 41.3% intergenic (‑303/+151) 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,7660GA41.3% 58.0 / 69.5 63intergenic (‑303/+151)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (12/14);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

ATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAC  >  minE/821697‑821835
                                                                     |                                                                     
aTCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCg                                                                       <  1:743751/70‑1 (MQ=39)
           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:3878576/56‑1 (MQ=255)
                           cAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                                          <  1:1323280/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:2005732/1‑56 (MQ=21)
                                aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGga                                                     >  1:2556154/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:3387662/35‑1 (MQ=17)
                                         ttCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                                                 <  1:410323/35‑1 (MQ=17)
                                          tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa                            >  1:2617331/1‑71 (MQ=38)
                                          tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa                            >  1:1973720/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)
                                              gtttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAtt                        <  1:3696396/69‑1 (MQ=255)
                                              gtttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAtt                        <  1:1741522/69‑1 (MQ=255)
                                               tttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAttt                       <  1:3726787/70‑1 (MQ=255)
                                               cTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                                                    >  1:3242477/1‑42 (MQ=255)
                                                ttttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatga                                                  <  1:1992919/43‑1 (MQ=255)
                                                ttttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGtt                                 <  1:3699251/60‑1 (MQ=255)
                                                ttttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAat                           <  1:1490457/66‑1 (MQ=255)
                                                 tttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGtt                                 >  1:546941/1‑59 (MQ=255)
                                                 tttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGtt                                 >  1:3897959/1‑59 (MQ=255)
                                                  ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTc                                            >  1:345701/1‑47 (MQ=255)
                                                  ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTAc                     <  1:2325240/70‑1 (MQ=255)
                                                  ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTAc                     <  1:1920230/70‑1 (MQ=255)
                                                  ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACt                    >  1:3843171/1‑71 (MQ=255)
                                                    cAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACtt                   <  1:1277908/70‑1 (MQ=255)
                                                    cAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACtt                   <  1:2223815/70‑1 (MQ=255)
                                                        aGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAg                                         >  1:3075574/1‑44 (MQ=255)
                                                        aGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTgg               >  1:847166/1‑70 (MQ=255)
                                                          cTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAag            <  1:1648264/71‑1 (MQ=255)
                                                              tGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTgg               <  1:1008128/64‑1 (MQ=255)
                                                               acAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGAt        <  1:2122423/69‑1 (MQ=255)
                                                                 aGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa     <  1:105634/71‑1 (MQ=255)
                                                                 aGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAttt                       <  1:29036/53‑1 (MQ=255)
                                                                    cGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCgg                                    >  1:119851/1‑37 (MQ=255)
                                                                     gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAat                           >  1:3023663/1‑45 (MQ=255)
                                                                     gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAc  >  1:1774006/1‑70 (MQ=255)
                                                                     |                                                                     
ATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAC  >  minE/821697‑821835

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: