Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,695,295 G→A 37.9% Q239* (CAG→TAG)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,2950GA37.9% 22.4 / 28.0 29Q239* (CAG→TAG) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (4/14);  new base A (6/5);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTA  >  minE/1695234‑1695362
                                                             |                                                                   
tAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCaa                                                             <  1:1816466/70‑1 (MQ=255)
      gCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGc                                                      <  1:1203146/71‑1 (MQ=255)
         ttACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAg                                                   <  1:718734/71‑1 (MQ=255)
           aCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                                 >  1:602288/1‑71 (MQ=255)
             gcgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAg                                                   <  1:1424051/67‑1 (MQ=255)
                    aTCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATa                                                           <  1:583887/52‑1 (MQ=255)
                    aTCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCAAt                                                            <  1:181968/51‑1 (MQ=255)
                     tCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCaa                                                             <  1:593367/49‑1 (MQ=255)
                      cAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCaa                                                             <  1:1874575/48‑1 (MQ=255)
                         gCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAg                                                   <  1:1955090/55‑1 (MQ=255)
                            aCGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCaa                                                             <  1:78488/42‑1 (MQ=255)
                            aCGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc                                <  1:624598/71‑1 (MQ=255)
                             cGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                               <  1:1679561/71‑1 (MQ=255)
                               gCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                             >  1:1680913/1‑71 (MQ=255)
                               gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCaa                                                             >  1:1425497/1‑39 (MQ=255)
                               gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCaa                                                             >  1:209760/1‑39 (MQ=255)
                               gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCaa                                                             >  1:337713/1‑39 (MQ=255)
                               gCTGTTATTCGGCCTTGTTACCTTTATTCTACTGCTCaa                                                             >  1:446826/1‑39 (MQ=255)
                                   ttATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                           <  1:2167374/69‑1 (MQ=255)
                                     aTTCGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGata                                   >  1:1499657/1‑59 (MQ=255)
                                        cGGCCTTGTTACCTTTATTCTACTGCTCAATAGAAGCCACAATCCCGCGCAGGATATTCaa                              >  1:1542151/1‑61 (MQ=255)
                                             ttGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                                            <  1:1777049/42‑1 (MQ=255)
                                             ttGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                                            <  1:1132944/42‑1 (MQ=255)
                                             ttGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt                         <  1:1582591/61‑1 (MQ=255)
                                              tGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                              >  1:1206230/1‑55 (MQ=255)
                                              tGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                              >  1:921727/1‑55 (MQ=255)
                                                     tttATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt                    <  1:192834/58‑1 (MQ=255)
                                                      ttATTCTACTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg          <  1:938833/67‑1 (MQ=255)
                                                         ttCTACTGCTCAATAGAAGCCAGAATCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTa  <  1:787778/71‑1 (MQ=255)
                                                             |                                                                   
TAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTA  >  minE/1695234‑1695362

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: