Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,044,281 T→C 56.3% N235N (AAT→AAC yqiK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,044,2810TC56.3% 21.7 / 60.8 80N235N (AAT→AACyqiKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/24);  new base C (34/11);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 1.05e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CAGAGCGCCGTCGCCGCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAT  >  minE/2044214‑2044331
                                                                   |                                                  
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:3002907/69‑1 (MQ=255)
               gcgAACGTAACGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTcgc                                    <  1:2469348/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:1171446/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:379699/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:800275/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:420559/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:614931/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:534172/1‑45 (MQ=255)
                         cGAAGTTGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACcgt                                                 >  1:587205/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:2654327/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:303857/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:2956381/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:345994/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:959056/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:961334/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:2872005/1‑69 (MQ=255)
                              ttGAACAGGATGTGGAAGTCGCAGTGCGTGAGAAAAACCGTGATGCACTTTCGCGCAAGTTGGAGATTg                     >  1:2753194/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:1939834/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:533079/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:2105802/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:1247106/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:2585797/1‑67 (MQ=255)
                                gAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTg                     >  1:928756/1‑67 (MQ=255)
                                    aGGATGTAGAAGTTGCGGTGCGTGAAAAAAATcgc                                                 <  1:1452943/35‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:3328560/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:140464/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:1656290/1‑69 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:132695/1‑69 (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:1296782/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:1033423/69‑1 (MQ=255)
                                            gAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       >  1:3049367/1‑69 (MQ=255)
                                             aaGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCg       <  1:3392373/68‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:988800/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:983665/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:611847/69‑1 (MQ=255)
                                              aGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGtt     <  1:601409/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)
                                                  gCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAt  <  1:2230436/68‑1 (MQ=255)
                                                                   |                                                  
CAGAGCGCCGTCGCCGCGAACGTAACGAAGTTGAACAGGATGTAGAAGTTGCGGTGCGTGAAAAAAATCGCGATGCGCTATCGCGCAAGCTGGAGATTGAACAGCAAGAAGCGTTTAT  >  minE/2044214‑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: