Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,420,852 A→T 54.5% R79R (CGT→CGA rffT ← TDP‑Fuc4NAc:lipidII transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,420,8520AT54.5% 4.7 / 48.5 44R79R (CGT→CGArffTTDP‑Fuc4NAc:lipidII transferase
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base T (11/13);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGGGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCAGCGATTTTTTCCCA  >  minE/2420785‑2420912
                                                                   |                                                             
tggGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:1067098/69‑1 (MQ=255)
tggGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:2122997/69‑1 (MQ=255)
tggGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:1136496/69‑1 (MQ=255)
tggGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:1728389/69‑1 (MQ=255)
tggGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:478020/69‑1 (MQ=255)
 gggCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAAc                                                             <  1:154240/68‑1 (MQ=255)
      tAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACgctgct                                                       <  1:1612380/69‑1 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:712307/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:1511360/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:152866/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:1608112/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:3228948/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:166700/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:2630389/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:806730/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:829179/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:1877132/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACg                                                   >  1:1993513/1‑67 (MQ=255)
            caccacTCAGCAGAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAAATCGCTGCTGACg                                                   >  1:975028/1‑68 (MQ=255)
                    agcagAGCCAGCCACAGTTTTGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:1029834/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2073240/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:1285890/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:420154/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:352843/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:3253527/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:3064301/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2995508/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2810000/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2697143/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:1683996/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2018374/60‑1 (MQ=255)
                    agcagAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCtttt                                         >  1:2299151/1‑69 (MQ=255)
                    agcagAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCtttt                                         >  1:493543/1‑69 (MQ=255)
                    agcagAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCtttt                                         >  1:1344158/1‑69 (MQ=255)
                             agccaCAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTcgcg                                <  1:1395158/69‑1 (MQ=255)
                             agccaCAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTcgcg                                <  1:833211/69‑1 (MQ=255)
                               ccaCAGTTTGGGATTGAACTGACCGTGGAAAAAAAATCGCTGCTGACgg                                                  <  1:2506587/49‑1 (MQ=38)
                                    gtgtGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCg                         >  1:928875/1‑69 (MQ=255)
                                            ttGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGc                  >  1:1651355/1‑68 (MQ=255)
                                             tGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCa                >  1:2565093/1‑69 (MQ=255)
                                             tGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCa                >  1:1738152/1‑69 (MQ=255)
                                                   gACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCAGCGAtt          >  1:2498870/1‑69 (MQ=255)
                                                   gACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCAGCGAtt          >  1:1132094/1‑69 (MQ=255)
                                                            agaaAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCAGCGATTTTTTCCCa  >  1:353960/1‑68 (MQ=255)
                                                                   |                                                             
TGGGCTTAATGCCACCACTCAGCAGAGCCAGCCACAGTGTGGGATTGAACTGACCGTGGAAGAAAAAACGCTGCTGACGGTTAGCTTTTGCTTTCGCGATGACCGCTTCCGCCAGCGATTTTTTCCCA  >  minE/2420785‑2420912

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: