Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,290,892 T→G 52.5% intergenic (+7/+446) yeeN → / ← asnW conserved hypothetical protein/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,290,8920TG52.5% 16.9 / 88.5 82intergenic (+7/+446)yeeN/asnWconserved hypothetical protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base T (29/10);  new base G (23/20);  total (52/30)
Fisher's exact test for biased strand distribution p-value = 6.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATA  >  minE/1290830‑1290957
                                                              |                                                                 
cccTTGAAGATGACGACGATGTACAAAAAGTTTATTATAACGTCGCAAATCTCTAATTATCTTTTaaa                                                              >  1:2115101/1‑68 (MQ=255)
cccTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTaaa                                                              >  1:2080844/1‑68 (MQ=255)
cccTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTaaa                                                              >  1:2971974/1‑68 (MQ=255)
cccTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTaaa                                                              >  1:305271/1‑68 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTGTATGTTAAAGAAATc                                                        >  1:149635/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAAt                                                         >  1:462056/1‑68 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2598754/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2066341/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2860206/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:3035151/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:1780357/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:1759815/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:3116037/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:3291405/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:280171/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:1666632/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:43457/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:1614436/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:1570795/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:459613/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2235119/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2245582/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2301699/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:637531/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:231632/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:790272/1‑69 (MQ=255)
     gAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATc                                                        >  1:2424155/1‑69 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:2508404/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:2523603/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1063739/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:2887970/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:3324244/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:451760/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:560572/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:568012/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:603864/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:749528/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:823649/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:869825/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:921475/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1832628/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1368883/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1661736/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:250755/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1673529/68‑1 (MQ=255)
       aGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:1723232/68‑1 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:1187627/1‑69 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:748085/1‑69 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:3419650/1‑69 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:3267152/1‑69 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:1706753/1‑69 (MQ=255)
         aTGACGACGATGTACAAAAAGTTTATCAAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTc                                                    >  1:2831309/1‑68 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  >  1:817744/1‑69 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  >  1:1433271/1‑69 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  >  1:1935955/1‑69 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  <  1:1648275/69‑1 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  <  1:1026724/69‑1 (MQ=255)
           gacgacGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCtt                                                  >  1:1118492/1‑69 (MQ=255)
            acgacgATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCttt                                                 >  1:2865887/1‑69 (MQ=255)
            acgacgATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCttt                                                 >  1:1259278/1‑69 (MQ=255)
             cgacgaTGTACAATAAGTTTATCATAACGTCGCAAATCTCTAATTATATGTTAAAGAAATCt                                                       <  1:97565/62‑1 (MQ=255)
              gacgaTGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTAc                                               >  1:90024/1‑69 (MQ=255)
              gacgaTGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTAc                                               >  1:2016324/1‑69 (MQ=255)
                     tACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCttt                                                                  >  1:3192215/1‑43 (MQ=255)
                     tACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGAt                                        >  1:1207084/1‑69 (MQ=255)
                     tACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGAt                                        >  1:2695749/1‑69 (MQ=255)
                                tATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAAtct                             >  1:1901302/1‑69 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:3413348/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:1362381/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:1156497/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:2030714/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:2802953/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:2187033/68‑1 (MQ=255)
                                    aTAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCat                          <  1:2193275/68‑1 (MQ=255)
                                                   tctAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATaaaa           >  1:2605932/1‑68 (MQ=255)
                                                         tATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataa    >  1:2691738/1‑69 (MQ=255)
                                                         tATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataa    <  1:2499436/69‑1 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAtaatataata  >  1:906220/1‑68 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataata  >  1:1485846/1‑68 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataata  >  1:389625/1‑68 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataata  >  1:349767/1‑68 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataata  >  1:1075329/1‑68 (MQ=255)
                                                            cTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAaatataat   >  1:2758914/1‑67 (MQ=255)
                                                              |                                                                 
CCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATA  >  minE/1290830‑1290957

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: