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

ATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCA  >  minE/821507‑821642
                                                                      |                                                                  
aTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                   <  1:1909009/71‑1 (MQ=34)
aTGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                    >  1:4937720/1‑70 (MQ=35)
     aacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACctc                                                              >  1:4271205/1‑71 (MQ=34)
     aacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACctc                                                              >  1:1018866/1‑71 (MQ=34)
     aacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACctc                                                              >  1:762428/1‑71 (MQ=34)
      acaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACctc                                                              <  1:5348219/70‑1 (MQ=35)
       caTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACctc                                                              <  1:2749988/69‑1 (MQ=35)
        aTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCtctc                                                            >  1:4164091/1‑70 (MQ=35)
           ttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACg                                                        >  1:1730592/1‑71 (MQ=34)
              cGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGt                                                       <  1:753174/69‑1 (MQ=35)
                gCTATTCTTGAATTG‑TTCTGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACct                                                               <  1:1570890/59‑1 (MQ=11)
                gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                      <  1:4483300/68‑1 (MQ=35)
                gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                      <  1:2987968/68‑1 (MQ=35)
                gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                      <  1:1553419/68‑1 (MQ=35)
                gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                      <  1:1546760/68‑1 (MQ=35)
                gCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                      <  1:3472529/68‑1 (MQ=35)
                 cTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCggg                                                  >  1:1886804/1‑71 (MQ=34)
                 cTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCgg                                                   <  1:175163/70‑1 (MQ=35)
                    ttCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGt                                               >  1:1870206/1‑71 (MQ=11)
                    ttCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGt                                               >  1:1290788/1‑71 (MQ=11)
                    ttCTTGAATTG‑TTCAGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCaa                                                          <  1:1278158/60‑1 (MQ=11)
                     tCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtt                                              >  1:5693301/1‑71 (MQ=34)
                      cTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                 >  1:2826786/1‑51 (MQ=35)
                       ttGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCATGTTTTACCTCTCAACGTGCGGGGGtt                                              <  1:3721880/69‑1 (MQ=25)
                       ttGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGttt                                                                    <  1:1544248/47‑1 (MQ=38)
                       ttGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  <  1:3133947/49‑1 (MQ=38)
                       ttGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGttt                                             >  1:756040/1‑70 (MQ=35)
                           aGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTctct                                        <  1:459328/71‑1 (MQ=34)
                            gtg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                 >  1:759367/1‑45 (MQ=38)
                             tg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtct                                          <  1:1454530/67‑1 (MQ=35)
                             tg‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtct                                          <  1:402836/67‑1 (MQ=35)
                              g‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:2406855/1‑42 (MQ=38)
                              g‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:3040576/1‑42 (MQ=38)
                              g‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:3507434/1‑42 (MQ=38)
                              g‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:3918708/1‑42 (MQ=38)
                               ttCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTc                                     <  1:2163854/69‑1 (MQ=11)
                                tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                    <  1:3780151/55‑1 (MQ=35)
                                  tGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:4328844/1‑39 (MQ=40)
                                  tGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTa                                                                  >  1:4371516/1‑39 (MQ=40)
                                  tGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTtctc                                         >  1:3455138/1‑64 (MQ=35)
                                    gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAc                                                                 >  1:875181/1‑38 (MQ=40)
                                    gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCtctc                                                            <  1:3748817/43‑1 (MQ=38)
                                    gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                 >  1:1714470/1‑70 (MQ=35)
                                     gttCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                 <  1:4008734/68‑1 (MQ=255)
                                      ttCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                              <  1:4248633/71‑1 (MQ=21)
                                        cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                             >  1:2561353/1‑70 (MQ=35)
                                        cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                             <  1:848013/70‑1 (MQ=35)
                                        cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                             <  1:1348151/70‑1 (MQ=35)
                                        cAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCgg                                                   >  1:616431/1‑48 (MQ=255)
                                        cAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCgg                                                   >  1:1947692/1‑48 (MQ=255)
                                         aaGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCa                            <  1:3345401/70‑1 (MQ=25)
                                           gATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATg                         >  1:2076834/1‑71 (MQ=34)
                                           gATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATg                         <  1:2913665/71‑1 (MQ=21)
                                            aTTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                        >  1:2321658/1‑71 (MQ=21)
                                            aTTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                        <  1:2163303/71‑1 (MQ=21)
                                              tAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                      >  1:2765249/1‑71 (MQ=21)
                                              tAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                      >  1:1831334/1‑71 (MQ=21)
                                                gCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccac                     <  1:1062496/70‑1 (MQ=35)
                                                      gTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGTGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGa               >  1:4432181/1‑70 (MQ=255)
                                                       ttCTGTTGTCAGGTTGTACCTCTCAACGTTCggggg                                                <  1:4297806/36‑1 (MQ=14)
                                                       ttCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGAt              <  1:4036196/70‑1 (MQ=25)
                                                       ttCTGTTGTAAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCtt                                       >  1:863469/1‑45 (MQ=255)
                                                        tCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                              >  1:5281153/1‑53 (MQ=16)
                                                         cTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATa             >  1:559051/1‑69 (MQ=14)
                                                         cTGTTGTCAGGTTGTACCTATCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                        >  1:205662/1‑58 (MQ=11)
                                                                  ggTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCa  <  1:4729267/71‑1 (MQ=25)
                                                                  aggtttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGc   <  1:1827786/67‑1 (MQ=25)
                                                                      |                                                                  
ATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCA  >  minE/821507‑821642

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: ‑