Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,044,263 T→C 57.0% V229V (GTT→GTC yqiK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,044,2630TC57.0% 24.8 / 58.4 79V229V (GTT→GTCyqiKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/23);  new base C (34/11);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAGCGCCGTCGCCGCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTT  >  minE/2044214‑2044328
                                                 |                                                                 
cAGAGCGCCGTCGCCGCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATc                                                >  1:540795/1‑69 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:1022997/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:3221523/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:1755652/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:225590/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:2374282/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:1478995/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:3023459/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:2469348/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:3002907/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:873387/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGCCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                 <  1:186851/69‑1 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2563649/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2934590/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:3056765/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2210772/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2873206/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2195939/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2165656/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:3113674/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:3092223/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2121326/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:3159343/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:3319535/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:379699/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:420559/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:534172/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:587205/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:614931/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:800275/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:1171446/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcg                                               >  1:1501528/1‑45 (MQ=255)
                         agAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                              >  1:2227568/2‑45 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCTCAAGTTGGAGATTg                  >  1:235295/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2753194/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2654327/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2872005/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:961334/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:959056/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2956381/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:303857/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:345994/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2656473/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2676581/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:2635116/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                  >  1:103551/1‑69 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:2585797/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:533079/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:928756/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:1247106/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:2105802/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                  >  1:1939834/1‑67 (MQ=255)
                                    aGGATGTAGAAGTTGCGGTGCGTGAAAAAAATcgc                                              <  1:1452943/35‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    <  1:140464/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    >  1:132695/1‑69 (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:1807895/1‑69 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    <  1:3283288/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    <  1:1844711/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    >  1:3049367/1‑69 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    <  1:1033423/69‑1 (MQ=255)
                                             aaGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg    <  1:3392373/68‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:611847/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:601409/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:983665/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:988800/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:590272/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:160048/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:417904/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:1773124/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:2159315/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:2242178/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:2413794/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:246174/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:2479879/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:1005451/69‑1 (MQ=255)
                                               gTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt  <  1:1505975/68‑1 (MQ=255)
                                                 |                                                                 
CAGAGCGCCGTCGCCGCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTT  >  minE/2044214‑2044328

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: