Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,533,003 G→A 57.1% G230G (GGC→GGT prmB ← N5‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,533,0030GA57.1% 8.0 / 35.6 49G230G (GGC→GGTprmBN5‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/18);  new base A (22/6);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 9.55e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGA  >  minE/1532945‑1533061
                                                          |                                                          
cTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCg                                                  >  1:1377236/1‑69 (MQ=255)
 tGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCg                                                  >  1:81175/1‑68 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:1780471/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:2600999/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:3078580/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:359206/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:428468/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:471664/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:1863726/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:1698573/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:830593/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:1186644/69‑1 (MQ=255)
             gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:1347696/69‑1 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:1713197/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:833071/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:69819/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:1824874/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:865044/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:629469/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:1434617/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:2158315/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:1327621/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:2622736/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt                                    >  1:1216770/1‑69 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:918297/68‑1 (MQ=255)
              agaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:3392177/68‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGTGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:2454019/69‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:1363975/69‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:146431/69‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:168046/69‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:1926338/69‑1 (MQ=255)
                         cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg                         <  1:2079523/69‑1 (MQ=255)
                          gTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa                                     <  1:345114/56‑1 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGtg                      >  1:259405/1‑67 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTc                      >  1:2980474/1‑68 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:2585475/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:1673822/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:756519/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:1801012/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:2300387/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:3392044/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:3191484/1‑69 (MQ=255)
                             aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg                     >  1:2729403/1‑69 (MQ=255)
                             aGTCTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGTCGGTATTCGTTTGGCAGGTCg                     >  1:244740/1‑69 (MQ=255)
                                             tgccagatgccagaACGAGTTCCGGCTCGTGGCGGTATTCGTTTGTCAg                         <  1:694455/49‑1 (MQ=37)
                                               ccAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATc    >  1:385840/1‑68 (MQ=255)
                                                 aGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa  <  1:2318447/68‑1 (MQ=255)
                                                 aGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa  <  1:2775895/68‑1 (MQ=255)
                                                 aGATGCCAGGCCCAGTTCCGGCTCGTGACGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa  <  1:1731235/68‑1 (MQ=255)
                                                  gATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTCGGCAGGTCGGACATATCTTCCGCATCGa  <  1:2356513/67‑1 (MQ=255)
                                                          |                                                          
CTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGA  >  minE/1532945‑1533061

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: