Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,527,117 G→A 39.9% G34G (GGG→GGA trmC → fused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,527,1170GA39.9% 29.1 / 42.3 65G34G (GGG→GGAtrmCfused 5‑methylaminomethyl‑2‑thiouridine forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base A (16/10);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCC  >  minE/1527057‑1527182
                                                            |                                                                 
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1750317/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1889511/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:3026482/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:2613969/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:2348161/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1796022/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1787324/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1776427/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1762118/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:569115/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:1528659/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:756651/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:797692/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:26210/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTCGAaga                                                           >  1:1011357/1‑69 (MQ=255)
tGCTGAGGGTACACCTGTTTCCAGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAaga                                                           >  1:2620146/1‑69 (MQ=255)
 gCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTTGAaga                                                           <  1:2321035/68‑1 (MQ=255)
 gCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAaga                                                           <  1:2315022/68‑1 (MQ=255)
          acacCTGTTTCCTGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGCAGAGACGCGatat                                                 >  1:1626070/1‑69 (MQ=255)
          acacCTGTTTCCTGAGATTTTACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATg                                                >  1:2230094/1‑69 (MQ=255)
          acacCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                 >  1:3092838/1‑69 (MQ=255)
          acacCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                 >  1:2730135/1‑69 (MQ=255)
          acacCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTAt                                                 >  1:775981/1‑69 (MQ=255)
                   tCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCt                                         >  1:3232994/1‑68 (MQ=255)
                   tCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCt                                         >  1:3015066/1‑68 (MQ=255)
                   tCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCt                                         >  1:472597/1‑68 (MQ=255)
                   tCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCt                                         >  1:2704776/1‑68 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:2606933/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:1143947/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:2807854/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:2130575/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:298858/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:2053813/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:32455/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:1637551/68‑1 (MQ=255)
                      cGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:1210891/68‑1 (MQ=255)
                       gagaTTTTGACGATGTCTATTTTTCCAACGATAACGGACTGGAAGAGACGCGTTATGTTTTTCTcgg                                      <  1:2651552/67‑1 (MQ=255)
                                     gTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTag                        >  1:2168777/1‑67 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     <  1:3058430/69‑1 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     <  1:3346109/69‑1 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     >  1:2558117/1‑69 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     >  1:2070570/1‑69 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     <  1:1809937/69‑1 (MQ=255)
                                      tCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAgg                     <  1:1820462/69‑1 (MQ=255)
                                        tATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTTTGTTTTTCTGGGAGGCAACCAATTAGAGGt                    >  1:645850/1‑68 (MQ=255)
                                        tATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCt                                         >  1:1859498/1‑47 (MQ=255)
                                        tATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTa                   >  1:1375563/1‑69 (MQ=255)
                                        tATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTa                   >  1:2099880/1‑69 (MQ=255)
                                                     aTAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGc      <  1:906408/69‑1 (MQ=255)
                                                     aTAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGc      <  1:1275732/69‑1 (MQ=255)
                                                     aTAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGc      <  1:2756322/69‑1 (MQ=255)
                                                     aTAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGc      <  1:1586100/69‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACTCtt              <  1:2630517/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:3184500/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:1894577/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:199770/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:593506/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:62409/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:634206/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:1422242/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:228163/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:2684762/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:949624/58‑1 (MQ=255)
                                                        aCGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:974038/58‑1 (MQ=255)
                                                         cGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCtt              <  1:705950/57‑1 (MQ=255)
                                                          gggCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATcc  <  1:186705/68‑1 (MQ=255)
                                                          gggCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATcc  <  1:3293364/68‑1 (MQ=255)
                                                          gggCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATcc  <  1:2370637/68‑1 (MQ=255)
                                                          gggCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATc   <  1:2089901/67‑1 (MQ=255)
                                                            |                                                                 
TGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCC  >  minE/1527057‑1527182

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: