Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,044,284 C→T 42.0% R236R (CGC→CGT yqiK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,044,2840CT42.0% 19.1 / 38.5 57R236R (CGC→CGTyqiKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/23);  new base T (13/11);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAT  >  minE/2044229‑2044331
                                                       |                                               
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:1022997/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:3221523/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:3023459/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:3002907/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:2469348/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:2374282/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:225590/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:873387/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:1478995/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:1755652/69‑1 (MQ=255)
gcgAACGTAACGAAGTTGAACAGGATGTGGAAGCCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:186851/69‑1 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCTCAAGTTGGAGATTg                     >  1:235295/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:961334/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:345994/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:303857/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:103551/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:2956381/1‑69 (MQ=255)
               ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:959056/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:2635116/1‑69 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:2585797/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:533079/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:1247106/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:1939834/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:928756/1‑67 (MQ=255)
                 gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:2105802/1‑67 (MQ=255)
                     aGGATGTAGAAGTTGCGGTGCGTGAAAAAAATcgc                                                 <  1:1452943/35‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:1296782/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:1807895/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:3328560/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:3283288/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:1033423/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:3049367/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:1844711/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:140464/69‑1 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:132695/1‑69 (MQ=255)
                             gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:1656290/1‑69 (MQ=255)
                              aaGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:3392373/68‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:1005451/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:988800/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:611847/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:417904/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:2479879/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:246174/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:2413794/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:2242178/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:2159315/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:1773124/69‑1 (MQ=255)
                               aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:160048/69‑1 (MQ=255)
                                gTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:1505975/68‑1 (MQ=255)
                                   gCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAt  <  1:2230436/68‑1 (MQ=255)
                                                       |                                               
GCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAT  >  minE/2044229‑2044331

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: