Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,788 T→G 37.8% intergenic (‑325/+129) 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,7880TG37.8% 16.4 / 62.2 53intergenic (‑325/+129)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base G (10/10);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGT  >  minE/821723‑821854
                                                                 |                                                                  
gCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCtt                                                                <  1:2229396/70‑1 (MQ=255)
 cAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTc                                                               <  1:57373/70‑1 (MQ=255)
    tAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGc                                                           >  1:2248697/1‑71 (MQ=255)
     gaaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTc                                                               >  1:1806389/2‑66 (MQ=21)
          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)
                tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa                                               >  1:2617331/1‑71 (MQ=38)
                tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCaa                                               >  1:1973720/1‑71 (MQ=38)
                tCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCa                                                >  1:2305213/1‑70 (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)
                      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:1920230/70‑1 (MQ=255)
                        ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTAc                                        <  1:2325240/70‑1 (MQ=255)
                        ttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACt                                       >  1:3843171/1‑71 (MQ=255)
                          cAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACtt                                      <  1:2223815/70‑1 (MQ=255)
                          cAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACtt                                      <  1:1277908/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)
                                            tttCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTtgt                                    <  1:1373198/54‑1 (MQ=255)
                                            tttCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaa                   <  1:3172004/71‑1 (MQ=255)
                                             ttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                        >  1:773758/1‑65 (MQ=255)
                                              tCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATAtt                                           >  1:403520/1‑45 (MQ=255)
                                              tCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                  >  1:278430/1‑70 (MQ=255)
                                              tCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTAc                                        >  1:3384180/1‑48 (MQ=255)
                                              tCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTAc                                        >  1:10341/1‑48 (MQ=255)
                                              tCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTgg                                  >  1:3869755/1‑54 (MQ=255)
                                                aTGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAAt               <  1:3394965/71‑1 (MQ=255)
                                                 tGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATg              <  1:662431/71‑1 (MQ=255)
                                                  gAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATg              <  1:2399610/70‑1 (MQ=255)
                                                    gTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGaa                                >  1:1656257/1‑50 (MQ=255)
                                                    gTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGcc           <  1:3825580/71‑1 (MQ=255)
                                                       ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGc         <  1:1552748/70‑1 (MQ=255)
                                                       ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa        <  1:2069879/71‑1 (MQ=255)
                                                         ctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa        <  1:1268934/69‑1 (MQ=255)
                                                         ctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATaaa                        >  1:3393867/1‑53 (MQ=255)
                                                         ctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCg          >  1:807849/1‑67 (MQ=255)
                                                            gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg   >  1:668563/1‑71 (MQ=255)
                                                             atgatgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt  >  1:535585/1‑71 (MQ=255)
                                                             atgatgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt  <  1:2425254/71‑1 (MQ=255)
                                                              tgaggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGt  >  1:1029920/1‑70 (MQ=255)
                                                                atgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGaa                                <  1:3756016/38‑1 (MQ=255)
                                                                atgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCt    <  1:2162208/66‑1 (MQ=255)
                                                                 |                                                                  
GCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGT  >  minE/821723‑821854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: