Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,623,821 G→A 8.6% A292A (GCG→GCA ftsY → fused Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,623,8210GA 8.6% 217.1 / 8.3 70A292A (GCG→GCAftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (39/25);  new base A (3/3);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTGAAAGAAGAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGTAAGCTGGC  >  minE/2623756‑2623886
                                                                 |                                                                 
ctgAAAGAAGAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGtt                                                              >  1:3567320/1‑71 (MQ=255)
  gAAAGAAGAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtg                                                            >  1:2157909/1‑71 (MQ=255)
    aagaagAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTTAAGGCAAAGCGCCGttt                                                             <  1:1075418/68‑1 (MQ=255)
    aagaagAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCACCGTTtgtg                                                          <  1:1405413/71‑1 (MQ=255)
     agaagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtg                                                            >  1:1678123/1‑68 (MQ=255)
     agaagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGa                                                         >  1:3926541/1‑71 (MQ=255)
     agaagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGa                                                         >  1:117104/1‑71 (MQ=255)
     agaagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGa                                                         <  1:45744/71‑1 (MQ=255)
       aagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtg                                                            >  1:3626892/1‑66 (MQ=255)
       aagaGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATc                                                       >  1:202242/1‑71 (MQ=255)
        agagATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATcc                                                      >  1:1361989/1‑71 (MQ=255)
         gagaTGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCt                                                     >  1:2644139/1‑71 (MQ=255)
           gaTGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGa                                                   <  1:1686852/71‑1 (MQ=255)
            aTGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAt                                                  <  1:2480320/71‑1 (MQ=255)
            aTGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAt                                                  <  1:3400927/71‑1 (MQ=255)
              gggCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAtg                                                 >  1:3241878/1‑70 (MQ=255)
               ggCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTTAtggt                                               <  1:2782800/71‑1 (MQ=255)
               ggCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAtggt                                               <  1:3142293/71‑1 (MQ=255)
                gCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtgtg                                                          >  1:3370383/1‑59 (MQ=255)
                    gaTTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCg                                          >  1:2565754/1‑71 (MQ=255)
                    gaTTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCACCGTTTGTGATCCTGAtg                                                 >  1:1420529/1‑64 (MQ=255)
                     aTTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGt                                         >  1:1886931/1‑71 (MQ=255)
                      ttCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTc                                        <  1:3574073/71‑1 (MQ=255)
                        cTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATc                                                       >  1:2028806/1‑54 (MQ=255)
                        cTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCaa                                      >  1:230052/1‑71 (MQ=255)
                           gCGAAAGTCGATGAGCCTCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   <  1:806026/71‑1 (MQ=255)
                           gCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   >  1:2181137/1‑71 (MQ=255)
                           gCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   >  1:1279079/1‑71 (MQ=255)
                           gCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   >  1:4044785/1‑71 (MQ=255)
                           gCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   >  1:3219552/1‑71 (MQ=255)
                             gAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGgtg                                 <  1:2656290/71‑1 (MQ=255)
                              aaaGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtgt                                                           <  1:4119764/44‑1 (MQ=255)
                              aaaGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTtgt                                                           <  1:1754351/44‑1 (MQ=255)
                                aGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGa                                                   >  1:1615238/1‑50 (MQ=255)
                                  tCGATGAGCCGCTGAATGTTGAAGGCAAAGCACCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTggg                             >  1:1413920/1‑70 (MQ=255)
                                    gATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGgt                                  >  1:1512944/1‑63 (MQ=255)
                                     aTGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGgt                                  <  1:1988669/62‑1 (MQ=255)
                                     aTGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaa                          <  1:1673415/70‑1 (MQ=255)
                                     aTGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaa                          >  1:2658877/1‑70 (MQ=255)
                                       gAGCCTCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAAc                       >  1:3065261/1‑71 (MQ=255)
                                          ccGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATcc                                                      <  1:505333/37‑1 (MQ=255)
                                            gCTGAATGTTGAAGTCAAAGCACCGTTTGTGATCCTGATGGTGGGCGTCAACGgtg                                 <  1:4026502/56‑1 (MQ=255)
                                            gCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCacga                  >  1:1500215/1‑71 (MQ=255)
                                            gCTGAATGTTGAAGGCAAAGCACCGTTTGTGATCCTGATGGTGGGCGTCAACGgtg                                 <  1:2759386/56‑1 (MQ=255)
                                             cTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaaaa                        >  1:2916052/1‑64 (MQ=255)
                                             cTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaa                          <  1:2637882/62‑1 (MQ=255)
                                               gAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAtgg                                                <  1:304333/38‑1 (MQ=255)
                                               gAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCAcgacga               >  1:2055416/1‑71 (MQ=255)
                                               gAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCAcgacga               >  1:68042/1‑71 (MQ=255)
                                                aaTGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCAcgacga               <  1:2587058/70‑1 (MQ=255)
                                                aaTGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGGCAACGGTGTGGGTAAAACCAcgacga               <  1:701953/70‑1 (MQ=255)
                                                 aTGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACg                                    >  1:1483938/1‑48 (MQ=255)
                                                  tGTTGAAGGCAAAGCGCCGTTTGTGATCCTGAtgt                                                <  1:199255/35‑2 (MQ=255)
                                                  tGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGgt                                  >  1:2165388/1‑49 (MQ=255)
                                                  tGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaa                          <  1:3517909/57‑1 (MQ=255)
                                                  tGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTg            >  1:2475607/1‑71 (MQ=255)
                                                  tGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTg            >  1:2736743/1‑71 (MQ=255)
                                                  tGTTGAAGGCAAAGCACCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaaaa                        >  1:3211991/1‑59 (MQ=255)
                                                    ttGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCaa                                      >  1:2466006/1‑43 (MQ=255)
                                                    ttGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTgg           >  1:2465374/1‑70 (MQ=255)
                                                    ttGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTgg           >  1:242259/1‑70 (MQ=255)
                                                    ctGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTaaaa                        >  1:1623273/2‑57 (MQ=255)
                                                     tGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGt          <  1:1025800/70‑1 (MQ=255)
                                                      gAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGTa         <  1:1173799/70‑1 (MQ=255)
                                                         ggCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCacg                   <  1:3067373/57‑1 (MQ=255)
                                                           cAAAGCGCCGTTTGTTATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCacga                  >  1:3676393/1‑56 (MQ=255)
                                                           cAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACgg                                   >  1:2685781/1‑39 (MQ=255)
                                                           cAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCacgac                 >  1:3886908/1‑57 (MQ=255)
                                                           cAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGTAAGCTgg   >  1:4004999/1‑71 (MQ=255)
                                                            aaaGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGTAAGCTGgc  <  1:1991526/71‑1 (MQ=255)
                                                               gcgcCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCacga                  >  1:2815959/1‑52 (MQ=255)
                                                                 |                                                                 
CTGAAAGAAGAGATGGGCGAGATTCTGGCGAAAGTCGATGAGCCGCTGAATGTTGAAGGCAAAGCGCCGTTTGTGATCCTGATGGTGGGCGTCAACGGTGTGGGTAAAACCACGACGATTGGTAAGCTGGC  >  minE/2623756‑2623886

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