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

GCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTC  >  minE/821493‑821603
                                          |                                                                     
gccgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTc                                          <  1:3827101/71‑1 (MQ=11)
  cgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtg                                        <  1:1149623/71‑1 (MQ=34)
   gccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgt                                       <  1:1230947/71‑1 (MQ=34)
   gccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgt                                       >  1:1341759/1‑71 (MQ=34)
        gttgtAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTc                                   <  1:766228/70‑1 (MQ=35)
        gttgtAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCACCGTTCtgttg                                     >  1:1255598/1‑68 (MQ=25)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              >  1:3045701/1‑71 (MQ=34)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              >  1:3050300/1‑71 (MQ=34)
            tAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGt                               <  1:371555/70‑1 (MQ=35)
               tGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTa                           >  1:1253158/1‑71 (MQ=11)
               tGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                           <  1:196800/71‑1 (MQ=34)
               tGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                           >  1:2540806/1‑71 (MQ=34)
                       ttttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                   >  1:2809023/1‑71 (MQ=34)
                        tttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAAc                  >  1:2311261/1‑71 (MQ=34)
                           gCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGt                <  1:869883/70‑1 (MQ=11)
                              gCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCgg            <  1:2183978/71‑1 (MQ=11)
                              gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCgg            <  1:770188/71‑1 (MQ=34)
                              gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCgg            <  1:3169187/71‑1 (MQ=34)
                               cTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCggg           >  1:3039768/1‑71 (MQ=11)
                               cTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCggg           >  1:3641791/1‑71 (MQ=11)
                               cTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                          >  1:626545/1‑56 (MQ=35)
                                tATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCa                                  <  1:3120075/47‑1 (MQ=11)
                                tATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACctct                      >  1:135042/1‑59 (MQ=11)
                                tATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                          >  1:3402600/1‑55 (MQ=35)
                                 aTTCTTGTAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTc                                   <  1:2972425/45‑1 (MQ=255)
                                  ttCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGc              <  1:318933/65‑1 (MQ=11)
                                  ttCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGAGGGt        >  1:3699716/1‑71 (MQ=21)
                                   tCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgttgt                                    <  1:3756497/42‑1 (MQ=38)
                                    cTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGttt      >  1:808936/1‑71 (MQ=34)
                                     ttGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAAc                  <  1:506088/58‑1 (MQ=11)
                                     ttGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGttt      >  1:110065/1‑70 (MQ=11)
                                     ttGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                            <  1:3690061/48‑1 (MQ=38)
                                      tGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc    >  1:209765/1‑71 (MQ=34)
                                        aaGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtt                              <  1:546272/43‑1 (MQ=38)
                                        aaGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtctc  <  1:2216441/71‑1 (MQ=34)
                                         aGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgttg                                     <  1:298772/35‑1 (MQ=40)
                                         aGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                             >  1:1725563/1‑43 (MQ=38)
                                          |                                                                     
GCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTC  >  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: ‑