Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,257,879 G→A 58.2% A174A (GCG→GCA yecP → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,257,8790GA58.2% 10.2 / 29.7 43A174A (GCG→GCAyecPpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base A (12/13);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTC  >  minE/1257831‑1257941
                                                |                                                              
gtgCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTa                                             <  1:3223495/68‑1 (MQ=255)
gtgCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTa                                             <  1:2620751/68‑1 (MQ=255)
gtgCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTa                                             <  1:1798089/68‑1 (MQ=255)
gtgCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTa                                             <  1:2002597/68‑1 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTTTTACCGTTAGGt                                              >  1:82431/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:1260617/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:3156403/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:2908251/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:520903/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:524883/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:227375/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:785011/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:849036/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:91491/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:1769885/1‑66 (MQ=255)
 tgCCAGTTTGAAGCAGTGCGTAAACTGCTGGGTAACGATCAGCGCGCACATTTGTTACCGTTAGGt                                              >  1:1505485/1‑66 (MQ=255)
         tGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAgg                                               <  1:13696/57‑1 (MQ=255)
         tGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACt                                   >  1:3287905/1‑69 (MQ=255)
         tGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACt                                   >  1:475724/1‑69 (MQ=255)
                      aaaCTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTa                                             <  1:1134256/46‑1 (MQ=255)
                              gggTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCttt                <  1:1362677/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGa                                         >  1:2023828/1‑39 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:1647354/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:1598106/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:1962775/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:640055/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:721607/67‑1 (MQ=255)
                                 tAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATAc           >  1:3082598/1‑69 (MQ=255)
                                      aTCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCt      >  1:3394148/1‑69 (MQ=255)
                                      aTCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCt      >  1:3271603/1‑69 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:543117/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:557763/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:96569/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:3231875/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2884370/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2739047/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2543278/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:243650/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2164488/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2127516/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:191607/56‑1 (MQ=255)
                                        cAGCGCGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:1875335/56‑1 (MQ=255)
                                        cAGCGAGCACATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCtt                 <  1:2930588/56‑1 (MQ=255)
                                           cgcgcgcATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTc  <  1:3144619/68‑1 (MQ=255)
                                           cgcgcgcATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTc  <  1:2520542/68‑1 (MQ=255)
                                           cgcgcgcATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTc  <  1:2454179/68‑1 (MQ=255)
                                           cgcgcgcATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTc  <  1:1742155/68‑1 (MQ=255)
                                               cgcATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGAt             <  1:2977356/53‑1 (MQ=255)
                                                |                                                              
GTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTC  >  minE/1257831‑1257941

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: