Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,309,899 G→A 19.0% R31R (CGC→CGT yefM ← antitoxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,309,8990GA19.0% 144.5 / 26.3 63R31R (CGC→CGTyefMantitoxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (7/5);  total (36/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.35e-01

AAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCT  >  minE/1309835‑1309968
                                                                |                                                                     
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                                   >  1:841505/1‑69 (MQ=255)
 aGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGtaa                                                                <  1:51469/71‑1 (MQ=255)
 aGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGtaa                                                                >  1:2934828/1‑71 (MQ=255)
 aGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGta                                                                 <  1:392709/70‑1 (MQ=255)
   ccGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGtaata                                                              <  1:2364932/71‑1 (MQ=255)
      tctcTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                           <  1:3752674/71‑1 (MQ=255)
       ctctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGa                                                          >  1:1427803/1‑71 (MQ=255)
        tctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                         >  1:816267/1‑71 (MQ=255)
        tctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                         >  1:416646/1‑71 (MQ=255)
        actTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGAt                                                         <  1:3079396/70‑1 (MQ=255)
          ttCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                           <  1:1265221/67‑1 (MQ=255)
          ttCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGa                                                          <  1:3804850/68‑1 (MQ=255)
          ttCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGAt                                                         <  1:3848843/69‑1 (MQ=255)
          ttCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGAt                                                         <  1:2265764/69‑1 (MQ=255)
               gCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCACCATTCTGGCGAGTAATAAGGATCGGGGCa                                                  >  1:3920795/1‑71 (MQ=255)
                cGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCAt                                                 <  1:2464254/71‑1 (MQ=255)
                      gTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                         >  1:2962877/1‑57 (MQ=255)
                       tATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                          <  1:823806/71‑1 (MQ=255)
                       tATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                          <  1:1725498/71‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCaa                                        >  1:3182559/1‑71 (MQ=255)
                           cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGtaata                                                              <  1:4115021/47‑1 (MQ=255)
                            ttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCggg                                                     >  1:238495/1‑55 (MQ=255)
                            ttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                         >  1:846383/1‑51 (MQ=255)
                            ttcGAGGGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCg                                                       >  1:7725/1‑53 (MQ=255)
                               gAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                          <  1:3897882/63‑1 (MQ=255)
                                 gTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCttt                                 >  1:1964067/1‑70 (MQ=255)
                                    aCATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                              <  1:1933302/70‑1 (MQ=255)
                                     cATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatc                             <  1:206156/70‑1 (MQ=255)
                                      aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGAGCGGGGCATGATCTTCAACGGCTTtcatc                             <  1:44153/69‑1 (MQ=255)
                                       tCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                           <  1:1290463/38‑1 (MQ=255)
                                       tCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCgg                                                      <  1:3405641/43‑1 (MQ=255)
                                       tCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                <  1:2769034/49‑1 (MQ=255)
                                       tCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGAt                                              <  1:448531/51‑1 (MQ=255)
                                       tCAGAACACAAGCCTCTCCATTCTGACGAGtaata                                                              <  1:3161205/35‑1 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATc                                             >  1:286869/1‑50 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCt                                            >  1:3834369/1‑51 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCaa                                        >  1:2899844/1‑55 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCaa                                        >  1:811892/1‑55 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                           >  1:3301298/1‑68 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgt                        >  1:3952156/1‑71 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgt                        >  1:93541/1‑71 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttg                         <  1:940401/70‑1 (MQ=255)
                                           aaCACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgttg                      >  1:3457566/1‑71 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGa                  >  1:1897116/1‑71 (MQ=255)
                                                 aGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGAc                 >  1:1339638/1‑70 (MQ=255)
                                                 aGCCACTCCATTCTGGCGAGTAATAAGGATCGGGGCATGa                                               >  1:1962070/1‑40 (MQ=255)
                                                    ctctCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACaaa              >  1:850977/1‑70 (MQ=255)
                                                      ctcCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTc           >  1:3345463/1‑71 (MQ=255)
                                                       tcCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACaaa              >  1:2216262/1‑67 (MQ=255)
                                                        ccATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACaaa              >  1:410939/1‑66 (MQ=255)
                                                        ccATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCt          >  1:2660372/1‑70 (MQ=255)
                                                         cATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCt          <  1:1344552/69‑1 (MQ=255)
                                                         cATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGa        >  1:2231723/1‑71 (MQ=255)
                                                         cATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCt          <  1:3312259/69‑1 (MQ=255)
                                                         cATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTg         <  1:149327/70‑1 (MQ=255)
                                                         cATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGa        >  1:3253101/1‑71 (MQ=255)
                                                            tCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtc                                >  1:2308623/1‑44 (MQ=255)
                                                            tCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtc                                >  1:2308210/1‑44 (MQ=255)
                                                            tCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtc                                >  1:1026285/1‑44 (MQ=255)
                                                             cTGGCGAGTAATAAGGATCGGGGCATGATCTTCaa                                        >  1:1342691/1‑35 (MQ=255)
                                                             cTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGAcgcg    >  1:1597187/1‑71 (MQ=255)
                                                             cTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGAcgcg    >  1:2644187/1‑71 (MQ=255)
                                                             cTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGAcgc     <  1:702546/70‑1 (MQ=255)
                                                               ggCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCt  >  1:1633446/1‑71 (MQ=255)
                                                                gCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAAt             <  1:3074019/59‑1 (MQ=255)
                                                                |                                                                     
AAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCT  >  minE/1309835‑1309968

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