Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,413,415 T→C 30.7% D270D (GAT→GAC aslA → acrylsulfatase‑like enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,413,4150TC30.7% 93.7 / 50.5 104D270D (GAT→GACaslAacrylsulfatase‑like enzyme
Reads supporting (aligned to +/- strand):  ref base T (26/46);  new base C (13/19);  total (39/65)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CATGCGGTGCGCGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCA  >  minE/2413353‑2413481
                                                              |                                                                  
cATGCGGTGCGCGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGa                                                               >  1:1085322/1‑68 (MQ=255)
cATGCGGTGCGCGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGa                                                               >  1:2853497/1‑68 (MQ=255)
cATGCGGTGCGCGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGa                                                               >  1:2472065/1‑68 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:1258093/69‑1 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:2206433/69‑1 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:3086427/69‑1 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:3236561/69‑1 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:2416293/69‑1 (MQ=255)
        gcgcGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCtg                                                      <  1:3152097/69‑1 (MQ=255)
                 cgAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTa                                             <  1:376870/69‑1 (MQ=255)
                  gAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTa                                             <  1:989890/68‑1 (MQ=255)
                  gAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTa                                             <  1:1449873/68‑1 (MQ=255)
                  gAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTa                                             <  1:2481312/68‑1 (MQ=255)
                  gAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTa                                             <  1:3134112/68‑1 (MQ=255)
                                 gCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCt                              >  1:1203967/1‑68 (MQ=255)
                                          aCGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGAt                     >  1:2587003/1‑68 (MQ=255)
                                          aCGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATg                    >  1:1672406/1‑69 (MQ=255)
                                          aCGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATg                    >  1:2149187/1‑69 (MQ=255)
                                             ccGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATgg                   >  1:2703225/1‑67 (MQ=255)
                                             ccGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATgg                   >  1:2176046/1‑67 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:751269/69‑1 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:1529766/69‑1 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:1430434/69‑1 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:3396592/69‑1 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:1887882/69‑1 (MQ=255)
                                               gAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaa               <  1:373373/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCTTCGACAAGATGGCGAagag            <  1:1211681/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:419357/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:2971978/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:3007053/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:901376/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:3044467/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:3091368/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:1657946/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:1682876/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:329614/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:456472/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:483202/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:219580/69‑1 (MQ=255)
                                                  atatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:3396530/69‑1 (MQ=255)
                                                   tatatGTAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:830148/68‑1 (MQ=255)
                                                   tatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:678591/68‑1 (MQ=255)
                                                   tatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:2578360/68‑1 (MQ=255)
                                                   tatatGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAagag            <  1:846844/68‑1 (MQ=255)
                                                       tGGAAGACCTGGATCATCGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:3089889/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGtt                                 >  1:652804/1‑43 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:3411628/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:2868422/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:849273/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:310642/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:577184/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:3156021/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:3184975/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:1359169/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:2061287/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:2318459/1‑69 (MQ=255)
                                                       tGGAAGACCTGGATCAACGCTGGATGGAATATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGTGATa       >  1:2378867/1‑69 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:99729/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:571640/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:1839037/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:1647589/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:773225/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:3046736/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:1419916/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATaa      >  1:2618923/1‑68 (MQ=255)
                                                         gAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATa       >  1:78036/1‑67 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGaaa              >  1:2417583/1‑57 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:2454463/1‑58 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:2300705/1‑58 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:750623/1‑58 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:699833/1‑58 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:2176610/1‑58 (MQ=255)
                                                           aGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAag              >  1:592881/1‑58 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1638802/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2059199/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:57320/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1743325/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:622972/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1630644/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1440745/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1589381/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1888868/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2634461/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2508495/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1123932/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2436110/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3206324/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3270370/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3319531/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3361253/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3368072/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2177389/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2146955/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1064973/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1894818/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:55390/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCGCGACAAGATGGCGAAGAGCGATAAACCa  <  1:698390/69‑1 (MQ=255)
                                                            gATCTGGATCAACGCTGGATGGACTATGGCGTTAAGGTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:2333658/69‑1 (MQ=255)
                                                             aTCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:1386555/68‑1 (MQ=255)
                                                             aTCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCa  <  1:3042014/68‑1 (MQ=255)
                                                              |                                                                  
CATGCGGTGCGCGGCGGCGAACAACAGGCCATTGCCGACATTACGCCGAAATATATGGAAGATCTGGATCAACGCTGGATGGACTATGGCGTTAAGTTCCTCGACAAGATGGCGAAGAGCGATAAACCA  >  minE/2413353‑2413481

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: