Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,044,293 A→T 36.4% L239L (CTA→CTT yqiK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,044,2930AT36.4% 50.8 / 54.8 66L239L (CTA→CTTyqiKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/24);  new base T (13/11);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGTT  >  minE/2044229‑2044356
                                                                |                                                               
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:1022997/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:225590/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:2374282/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:2469348/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:1755652/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:3002907/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:3023459/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:1478995/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:3221523/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:873387/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGCCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                                             <  1:186851/69‑1 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCTCAAGTTGGAGATTg                                              >  1:235295/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:303857/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:345994/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2956381/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2872005/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2753194/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2676581/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2656473/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2654327/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:959056/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:961334/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:2635116/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                                              >  1:103551/1‑69 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:533079/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:2105802/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:2585797/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:1247106/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:928756/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                                              >  1:1939834/1‑67 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                >  1:1656290/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:1296782/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:3328560/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:3283288/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                >  1:132695/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                >  1:3049367/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:140464/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:1844711/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                >  1:1807895/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:1033423/69‑1 (MQ=255)
                              aaGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg                                <  1:3392373/68‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:2242178/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:417904/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:988800/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:2159315/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:2413794/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:246174/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:590272/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:601409/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:983665/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:2479879/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:611847/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:1005451/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:160048/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:1773124/69‑1 (MQ=255)
                                gTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:1505975/68‑1 (MQ=255)
                                   gCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAt                           <  1:2230436/68‑1 (MQ=255)
                                                           gcgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt                              <  1:2692494/41‑1 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:3093610/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:1517398/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:2995482/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:2978116/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:2642862/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:2049186/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGtt  >  1:2325215/1‑68 (MQ=255)
                                                            cgcTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGt   >  1:3091072/1‑67 (MQ=255)
                                                                |                                                               
GCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTATGACGCTTGAGCAGGAGCAGCAGGTT  >  minE/2044229‑2044356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: