Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,576 T→G 41.9% intergenic (‑113/+341) 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,5760TG41.9% 26.3 / 56.8 43intergenic (‑113/+341)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base G (7/11);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  minE/821508‑821643
                                                                     |                                                                   
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)
                 tATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACctct                                                              >  1:135042/1‑59 (MQ=11)
                 tATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                  >  1:3402600/1‑55 (MQ=35)
                   ttCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGc                                                      <  1:318933/65‑1 (MQ=11)
                   ttCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGAGGGt                                                >  1:3699716/1‑71 (MQ=21)
                     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)
                         aaGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtctc                                          <  1:2216441/71‑1 (MQ=34)
                          aGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     >  1:1725563/1‑43 (MQ=38)
                               tCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTc                                      <  1:962963/69‑1 (MQ=11)
                                   gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCtctc                                                             <  1:1157142/43‑1 (MQ=38)
                                   gTTTCAAGATTAGACCCCGTTCTGTTGTCAGGTTTTACCTCTCaa                                                           >  1:3900291/1‑45 (MQ=255)
                                       cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                              <  1:2743903/70‑1 (MQ=35)
                                        aaGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCa                             <  1:1599961/70‑1 (MQ=35)
                                        aaGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCa                             <  1:675432/70‑1 (MQ=35)
                                         aGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAg                                   >  1:3529007/1‑63 (MQ=25)
                                         aGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                           <  1:2480140/71‑1 (MQ=21)
                                             tAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                       >  1:380991/1‑71 (MQ=21)
                                              aGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                       >  1:1712516/1‑70 (MQ=35)
                                                cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                              <  1:3316450/61‑1 (MQ=35)
                                                cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                              <  1:877587/61‑1 (MQ=35)
                                                cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacca                    <  1:1957681/71‑1 (MQ=34)
                                                 ccccGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTtct                                           <  1:2517534/47‑1 (MQ=255)
                                                 ccccGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAg                   <  1:2118913/71‑1 (MQ=21)
                                                  cccGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                              <  1:2147854/59‑1 (MQ=255)
                                                     gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGa                <  1:3247547/70‑1 (MQ=35)
                                                      ttCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATa              <  1:1500648/71‑1 (MQ=34)
                                                      ttCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAGTGCCACCAGGGATa              <  1:2770583/71‑1 (MQ=255)
                                                       tCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAggg                 <  1:2471799/67‑1 (MQ=18)
                                                                aGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCg     >  1:246674/1‑70 (MQ=33)
                                                                  gTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa  >  1:3363216/1‑71 (MQ=38)
                                                                    tttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGccccc      >  1:594259/1‑65 (MQ=31)
                                                                     |                                                                   
TGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  minE/821508‑821643

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: