Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,576 T→G 21.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,5760TG21.9% 57.5 / 17.5 32intergenic (‑113/+341)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (11/14);  new base G (3/4);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  minE/821506‑821643
                                                                       |                                                                   
aaTGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     >  1:1624421/1‑71 (MQ=34)
 aTGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:3822530/71‑1 (MQ=34)
 aTGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     >  1:3195340/1‑70 (MQ=35)
 aTGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     <  1:3179352/70‑1 (MQ=35)
    acaacaTTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAcc                                                                 >  1:294815/1‑71 (MQ=34)
        caTTTTGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCtctc                                                             <  1:1452934/71‑1 (MQ=11)
             tGCGGCTATTCTTGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACg                                                         <  1:1620400/70‑1 (MQ=11)
              gCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:4028386/58‑1 (MQ=35)
               cGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                       <  1:3578907/70‑1 (MQ=35)
               cGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                       <  1:1608141/70‑1 (MQ=35)
                  cTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACct                                                                >  1:2128424/1‑58 (MQ=35)
                  cTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACct                                                                >  1:3029150/1‑58 (MQ=35)
                       cTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggggg                                                 <  1:1538781/68‑1 (MQ=35)
                        ttGAATTGTTCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGtt                                               <  1:1778157/69‑1 (MQ=11)
                          gAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                             <  1:3080060/69‑1 (MQ=35)
                          gAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc                                            <  1:3526143/70‑1 (MQ=35)
                             gtgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:556394/43‑1 (MQ=38)
                              tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     <  1:3424958/41‑1 (MQ=38)
                              tgtCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTAAGGTTTTACCTCTCaa                                                           <  1:893060/51‑1 (MQ=255)
                                 tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACg                                                         >  1:103941/1‑51 (MQ=35)
                                  cTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCtctc                                                             <  1:4014743/46‑1 (MQ=38)
                                     gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCa                                 <  1:1063978/71‑1 (MQ=34)
                                        tCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                  >  1:4119763/1‑35 (MQ=40)
                                        tCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                  >  1:105460/1‑35 (MQ=40)
                                         cAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                              <  1:929802/70‑1 (MQ=25)
                                           aGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                     <  1:239473/45‑1 (MQ=38)
                                             aTTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                            >  1:591947/1‑68 (MQ=35)
                                                aGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc                                            >  1:3145238/1‑49 (MQ=38)
                                                 gCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccac                      <  1:3844609/70‑1 (MQ=35)
                                                 gCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccac                      <  1:1566700/70‑1 (MQ=35)
                                                  cccccGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcaac                      >  1:906797/1‑69 (MQ=25)
                                                     ccGTCCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc                                            <  1:3622396/44‑1 (MQ=255)
                                                      cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCaa                                >  1:3990135/1‑55 (MQ=35)
                                                       gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                             >  1:551319/1‑41 (MQ=38)
                                                        ttCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGAt               <  1:3810372/70‑1 (MQ=35)
                                                          cTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                            >  1:3268292/1‑55 (MQ=9)
                                                          cTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                            >  1:1626369/1‑55 (MQ=9)
                                                                  aGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCg     >  1:2422109/1‑70 (MQ=18)
                                                                    gTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa  <  1:188836/71‑1 (MQ=38)
                                                                       |                                                                   
AATGACAACATTTTGCGGCTATTCTTGAAGTGTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  minE/821506‑821643

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