Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,884,080 G→A 11.9% L135L (CTG→TTG)  rumA ← 23S rRNA (uracil‑5)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,884,0800GA11.9% 124.9 / 4.2 66L135L (CTG→TTG) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base A (4/4);  total (36/30)
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.43e-01

GACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGT  >  minE/1884012‑1884144
                                                                    |                                                                
gACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAg                                                                  >  1:3284550/1‑69 (MQ=255)
gACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAg                                                                  >  1:3280897/1‑69 (MQ=255)
gACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAg                                                                  >  1:1585246/1‑69 (MQ=255)
gACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAg                                                                  >  1:3223843/1‑69 (MQ=255)
     tGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTa                                                             <  1:2066111/69‑1 (MQ=255)
     tGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTa                                                             <  1:2064839/69‑1 (MQ=255)
     tGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTa                                                             <  1:2907602/69‑1 (MQ=255)
     tGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTa                                                             <  1:400329/69‑1 (MQ=255)
      gTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTa                                                             <  1:1469654/68‑1 (MQ=255)
         aCTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAAcg                                                         >  1:465896/1‑69 (MQ=255)
         aCTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAAcg                                                         >  1:1345723/1‑69 (MQ=255)
         aCTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAAcg                                                         >  1:1194151/1‑69 (MQ=255)
         aCTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAAcg                                                         >  1:3029055/1‑69 (MQ=255)
             gAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAAcgcgc                                                      <  1:1630986/68‑1 (MQ=255)
               gCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAAcgcgcgcg                                                   >  1:2510535/1‑69 (MQ=255)
                            aaCCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAACGCGCGCGACGGCGATAGccc                                      <  1:37409/69‑1 (MQ=255)
                            aaCCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAACGCGCGCGACGGCGATAGccc                                      <  1:3237394/69‑1 (MQ=255)
                            aaCCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAACGCGCGCGACGGCGATAGccc                                      <  1:1319443/69‑1 (MQ=255)
                            aaCCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAAACTTAAACGCGCGCGACGGCGATAGccc                                      <  1:2935318/69‑1 (MQ=255)
                                ccATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAg                                   >  1:1762890/1‑68 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:2511953/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:964399/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:3163714/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:242017/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:2216860/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:87118/1‑45 (MQ=255)
                                         gTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGAc                                                 >  1:3301880/1‑45 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCgg                         >  1:3222684/1‑68 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:252022/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:2995935/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:2726562/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:1004034/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:506364/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:1812769/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:2270185/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:3368932/1‑69 (MQ=255)
                                          ttGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:3248570/1‑69 (MQ=255)
                                          ttGCTGAGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGc                        >  1:866/1‑69 (MQ=255)
                                              tgtgTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGAt                     <  1:710405/68‑1 (MQ=255)
                                                  ttttCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCAc                  <  1:1443729/67‑1 (MQ=255)
                                                  ttttCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCAc                  <  1:228525/67‑1 (MQ=255)
                                                     tCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGc                                              <  1:358465/36‑1 (MQ=255)
                                                     tCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGc                                              <  1:807956/36‑1 (MQ=255)
                                                     tCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGc                                              <  1:2685822/36‑1 (MQ=255)
                                                     tCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGc                                              <  1:2398776/36‑1 (MQ=255)
                                                     tCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGc                                              <  1:1826792/36‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:2940086/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:2763544/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:2736841/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:2214374/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:1068538/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:748581/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:1440012/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:3340140/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:1916190/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:1801348/69‑1 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:3103965/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:934119/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:1396066/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:62779/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:149737/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:1870987/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:3324491/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:1983150/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:2131981/1‑69 (MQ=255)
                                                              agtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgt    >  1:2763180/1‑69 (MQ=255)
                                                               gtTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACTTCgt    <  1:2002580/68‑1 (MQ=255)
                                                                ttCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCTGCGATCACTTCAGAGACTTCgt    <  1:1105190/67‑1 (MQ=255)
                                                                 tCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgtgt  >  1:148066/1‑68 (MQ=255)
                                                                 tCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCgtgt  >  1:104222/1‑68 (MQ=255)
                                                                    |                                                                
GACAATGTCACTGGAGCCCGCTTTGCGAAACCCCATCTGAAGTTGCTGTGTTTTCGGTAAGTAGTTCAGACTTAAACGCGCGCGACGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGT  >  minE/1884012‑1884144

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: