Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,760 G→A 37.5% intergenic (‑297/+157) 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,7600GA37.5% 75.8 / 62.9 64intergenic (‑297/+157)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base A (12/12);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAG  >  minE/821692‑821825
                                                                    |                                                                 
aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                                                 <  1:914898/71‑1 (MQ=34)
     aTCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCg                                                             <  1:743751/70‑1 (MQ=39)
             ttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                                <  1:1271260/59‑1 (MQ=35)
              tcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                                <  1:2102745/58‑1 (MQ=35)
              tcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                                                <  1:1585193/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:3387662/35‑1 (MQ=17)
                                              ttCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                                       <  1:410323/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:2129108/1‑39 (MQ=255)
                                                 cGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                                >  1:2953427/1‑39 (MQ=255)
                                                  ggCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCt                             >  1:2577897/1‑57 (MQ=255)
                                                   gtttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAtt              <  1:1741522/69‑1 (MQ=255)
                                                   gtttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATAtt              <  1:3696396/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:3897959/1‑59 (MQ=255)
                                                      tttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGtt                       >  1:546941/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)
                                                                    |                                                                 
AATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAG  >  minE/821692‑821825

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: