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

AATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  minE/821506‑821643
                                                                       |                                                                   
aaTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     >  1:1624421/1‑71 (MQ=34)
 aTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:3822530/71‑1 (MQ=34)
 aTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     >  1:3195340/1‑70 (MQ=35)
 aTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     <  1:3179352/70‑1 (MQ=35)
    acaacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAcc                                                                 >  1:294815/1‑71 (MQ=34)
        caTTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCtctc                                                             <  1:1452934/71‑1 (MQ=11)
             tGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACg                                                         <  1:1620400/70‑1 (MQ=11)
              gCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:4028386/58‑1 (MQ=35)
               cGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                       <  1:3578907/70‑1 (MQ=35)
               cGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                       <  1:1608141/70‑1 (MQ=35)
                  cTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACct                                                                >  1:2128424/1‑58 (MQ=35)
                  cTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACct                                                                >  1:3029150/1‑58 (MQ=35)
                       cTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggggg                                                 <  1:1538781/68‑1 (MQ=35)
                        ttGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGtt                                               <  1:1778157/69‑1 (MQ=11)
                          gAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                             <  1:3080060/69‑1 (MQ=35)
                          gAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtc                                            <  1:3526143/70‑1 (MQ=35)
                             gtg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                    <  1:556394/43‑1 (MQ=38)
                              tg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                     <  1:3424958/41‑1 (MQ=38)
                              tg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTAAGGTTTTACCTCTCaa                                                           <  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)
                                                                       |                                                                   
AATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAA  >  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: ‑