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

GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCA  >  minE/821760‑821891
                                                                  |                                                                 
gCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGAt                                                                <  1:1084642/70‑1 (MQ=255)
gCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGAt                                                                <  1:3781259/70‑1 (MQ=255)
 cAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                               <  1:2591327/70‑1 (MQ=255)
 cAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATa                                                               <  1:5597983/70‑1 (MQ=255)
         tCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaa                                                       <  1:2494345/70‑1 (MQ=255)
           aTGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaaa                                                     <  1:4944308/70‑1 (MQ=255)
                  ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                             <  1:1054012/71‑1 (MQ=255)
                  ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                             <  1:1508998/71‑1 (MQ=255)
                    ctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaa                                                       <  1:1141625/59‑1 (MQ=255)
                     tGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAc                                                          >  1:34747/1‑55 (MQ=255)
                     tGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                             <  1:2228891/68‑1 (MQ=255)
                      ggATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGc                                          >  1:2077133/1‑70 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                             >  1:4257911/1‑66 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAgg                                           <  1:4582793/68‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAgg                                           <  1:2680119/68‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGc                                          <  1:3805761/69‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCt                                         <  1:3346073/70‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg                                        >  1:2476154/1‑71 (MQ=255)
                       gatgagGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCt                                         <  1:2602855/70‑1 (MQ=255)
                        atgaggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGt                                       <  1:354236/71‑1 (MQ=255)
                           aggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaa                                                       <  1:5235309/50‑1 (MQ=255)
                           aggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaa                                                       <  1:4533099/50‑1 (MQ=255)
                            tgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGc                                                 <  1:4883770/57‑1 (MQ=255)
                             gCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTAccc                                  <  1:2757126/71‑1 (MQ=255)
                             gCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                             >  1:666443/1‑44 (MQ=255)
                              cTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTAcccc                                 >  1:2328042/1‑71 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                             >  1:1169744/1‑39 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaa                                                      >  1:3300256/1‑46 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGc                              <  1:2159948/70‑1 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGc                              <  1:4517571/70‑1 (MQ=255)
                                  cAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCg                             >  1:5073868/1‑71 (MQ=255)
                                   aGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTAc                                    >  1:3388011/1‑63 (MQ=255)
                                    gCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                             <  1:4728784/37‑1 (MQ=255)
                                      tCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCg                         >  1:1938616/1‑71 (MQ=255)
                                         ggTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGtt                      >  1:4949673/1‑71 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCt                                         <  1:4833049/51‑1 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCt                                         <  1:1016865/51‑1 (MQ=255)
                                            tGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGt                                       >  1:3390871/1‑51 (MQ=255)
                                             gCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCg                             >  1:3036666/1‑60 (MQ=255)
                                             gCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCgg                  >  1:263875/1‑71 (MQ=255)
                                                  tatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa             <  1:218579/71‑1 (MQ=255)
                                                    tttACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGt                       >  1:4858454/1‑59 (MQ=255)
                                                           gtgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCAt    >  1:5444806/1‑71 (MQ=255)
                                                              gAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCa  <  1:5211129/70‑1 (MQ=255)
                                                               aaGTGATAATGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGTGCGCATATTGCCATc   <  1:901305/68‑1 (MQ=14)
                                                                  agaTAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCa  >  1:253090/1‑66 (MQ=255)
                                                                  |                                                                 
GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCA  >  minE/821760‑821891

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: