Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,826 A→T 38.6% intergenic (‑363/+91) 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,8260AT38.6% 32.5 / 48.2 44intergenic (‑363/+91)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base A (9/18);  new base T (7/10);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCC  >  minE/821760‑821887
                                                                  |                                                             
acAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGAt                                                            <  1:2122423/69‑1 (MQ=255)
  aGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                         <  1:105634/71‑1 (MQ=255)
      gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAc                                                      >  1:1774006/1‑70 (MQ=255)
       tttCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaa                                                    <  1:3172004/71‑1 (MQ=255)
        ttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                         >  1:773758/1‑65 (MQ=255)
         tCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                                                   >  1:278430/1‑70 (MQ=255)
           aTGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAAt                                                <  1:3394965/71‑1 (MQ=255)
            tGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATg                                               <  1:662431/71‑1 (MQ=255)
             gAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATg                                               <  1:2399610/70‑1 (MQ=255)
               gTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGcc                                            <  1:3825580/71‑1 (MQ=255)
                  ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGc                                          <  1:1552748/70‑1 (MQ=255)
                  ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                         <  1:2069879/71‑1 (MQ=255)
                    ctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                         <  1:1268934/69‑1 (MQ=255)
                    ctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATaaa                                                         >  1:3393867/1‑53 (MQ=255)
                    ctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCg                                           >  1:807849/1‑67 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg                                    >  1:668563/1‑71 (MQ=255)
                        atgatgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt                                   <  1:2425254/71‑1 (MQ=255)
                        atgatgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt                                   >  1:535585/1‑71 (MQ=255)
                         tgaggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGt                                   >  1:1029920/1‑70 (MQ=255)
                           atgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCt                                     <  1:2162208/66‑1 (MQ=255)
                             gCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                         <  1:2080475/60‑1 (MQ=255)
                                 ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaaaa                                                 <  1:3702351/48‑1 (MQ=255)
                                 ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGc                          >  1:3201175/1‑71 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCg                         <  1:209352/71‑1 (MQ=255)
                                    gCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCgg                        >  1:55312/1‑70 (MQ=255)
                                     cTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCg                         <  1:3689502/68‑1 (MQ=255)
                                       cGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGAc                                                      <  1:3634067/37‑1 (MQ=255)
                                        gggTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCg                                           <  1:3682159/47‑1 (MQ=255)
                                        gggTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCg                                           <  1:1124784/47‑1 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAAt                                                <  1:2931438/40‑1 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAAt                                                <  1:2992327/40‑1 (MQ=255)
                                           ttGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTc                <  1:1627388/71‑1 (MQ=255)
                                           ctGCCAATATTTACTTGTGGAAGTGATAAAGACaaaa                                                  <  1:1712164/36‑1 (MQ=39)
                                             gCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGc                                             >  1:2977712/1‑40 (MQ=255)
                                               cAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgc            >  1:3890915/1‑71 (MQ=255)
                                               cAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt                                   <  1:2076553/48‑1 (MQ=255)
                                                 atatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                         <  1:3503344/40‑1 (MQ=255)
                                                 atatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgcg           >  1:2478154/1‑70 (MQ=255)
                                                  tatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGc                                      >  1:165946/1‑42 (MQ=255)
                                                   atTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa         >  1:199331/1‑70 (MQ=255)
                                                       aCTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTg    >  1:2829459/1‑71 (MQ=255)
                                                           gtgGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGcc  <  1:864509/69‑1 (MQ=255)
                                                            tgGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTAc                                <  1:535924/38‑1 (MQ=255)
                                                               aaGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGcc                      <  1:29224/45‑1 (MQ=255)
                                                                  agaTAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGc                          >  1:848335/1‑38 (MQ=255)
                                                                  |                                                             
GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCC  >  minE/821760‑821887

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: