Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,535 G→T 29.7% intergenic (‑72/+382) 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,5350GT29.7% 71.8 / 19.6 37intergenic (‑72/+382)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base T (5/6);  total (18/19)
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

GCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTC  >  minE/821493‑821603
                                          |                                                                     
gccgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTc                                          <  1:3827101/71‑1 (MQ=11)
  cgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCtg                                        <  1:1149623/71‑1 (MQ=34)
   gccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCtgt                                       <  1:1230947/71‑1 (MQ=34)
   gccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCtgt                                       >  1:1341759/1‑71 (MQ=34)
        gttgtAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTc                                   <  1:766228/70‑1 (MQ=35)
        gttgtAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCACCGTTCtgttg                                     >  1:1255598/1‑68 (MQ=25)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              >  1:3045701/1‑71 (MQ=34)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              >  1:3050300/1‑71 (MQ=34)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGt                               <  1:371555/70‑1 (MQ=35)
               tGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTa                           >  1:1253158/1‑71 (MQ=11)
               tGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                           <  1:196800/71‑1 (MQ=34)
               tGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                           >  1:2540806/1‑71 (MQ=34)
                       ttttGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                   >  1:2809023/1‑71 (MQ=34)
                        tttGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAAc                  >  1:2311261/1‑71 (MQ=34)
                           gCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGt                <  1:869883/70‑1 (MQ=11)
                              gCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCgg            <  1:2183978/71‑1 (MQ=11)
                              gCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCgg            <  1:770188/71‑1 (MQ=34)
                              gCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCgg            <  1:3169187/71‑1 (MQ=34)
                               cTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCggg           >  1:3039768/1‑71 (MQ=11)
                               cTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCggg           >  1:3641791/1‑71 (MQ=11)
                               cTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                          >  1:626545/1‑56 (MQ=35)
                                tATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCa                                  <  1:3120075/47‑1 (MQ=11)
                                tATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACctct                      >  1:135042/1‑59 (MQ=11)
                                tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                          >  1:3402600/1‑55 (MQ=35)
                                 aTTCTTGTAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTc                                   <  1:2972425/45‑1 (MQ=255)
                                  ttCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGc              <  1:318933/65‑1 (MQ=11)
                                  ttCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGAGGGt        >  1:3699716/1‑71 (MQ=21)
                                   tCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCtgttgt                                    <  1:3756497/42‑1 (MQ=38)
                                    cTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGttt      >  1:808936/1‑71 (MQ=34)
                                     ttGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAAc                  <  1:506088/58‑1 (MQ=11)
                                     ttGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGttt      >  1:110065/1‑70 (MQ=11)
                                     ttGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                            <  1:3690061/48‑1 (MQ=38)
                                      tGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc    >  1:209765/1‑71 (MQ=34)
                                        aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              <  1:546272/43‑1 (MQ=38)
                                        aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtctc  <  1:2216441/71‑1 (MQ=34)
                                         aGTGTCTGGTTTCAAGATTAGCCCCCGTTCtgttg                                     <  1:298772/35‑1 (MQ=40)
                                         aGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                             >  1:1725563/1‑43 (MQ=38)
                                          |                                                                     
GCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTC  >  minE/821493‑821603

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: