Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,023,466 G→T 43.2% V99V (GTC→GTA uidA ← beta‑D‑glucuronidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,023,4660GT43.2% 71.9 / 117.2 127V99V (GTC→GTAuidAbeta‑D‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base G (10/62);  new base T (34/21);  total (44/83)
Fisher's exact test for biased strand distribution p-value = 2.09e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

ACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGG  >  minE/1023398‑1023520
                                                                    |                                                      
aCATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTg                                                        >  1:604431/1‑69 (MQ=255)
aCATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTg                                                        >  1:3211234/1‑69 (MQ=255)
aCATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTg                                                        >  1:3282362/1‑69 (MQ=255)
aCATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTg                                                        >  1:3286084/1‑69 (MQ=255)
aCATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTg                                                        >  1:436126/1‑69 (MQ=255)
                cGGCTTCAATTGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGt                                        <  1:2417848/69‑1 (MQ=255)
                cGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGt                                        <  1:2770486/69‑1 (MQ=255)
                cGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGt                                        <  1:912615/69‑1 (MQ=255)
                    ttCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATg                                    <  1:2928534/69‑1 (MQ=255)
                     tCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCTCACTTTGCCGTAATGa                                   <  1:2710678/69‑1 (MQ=255)
                     tCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGa                                   <  1:905733/69‑1 (MQ=255)
                     tCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGa                                   <  1:2776031/69‑1 (MQ=255)
                     tCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGa                                   <  1:2935554/69‑1 (MQ=255)
                         aTGGCGTATAGCGGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTg                                <  1:2321867/68‑1 (MQ=255)
                         aTGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTg                                <  1:2022487/68‑1 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTaa                                      >  1:941760/1‑56 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:57009/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2581890/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:281525/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:891369/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2667327/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2496246/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1200556/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:865815/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1334696/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2958836/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2368869/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:831197/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1078627/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:2208866/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:539178/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:15962/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:3398234/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:34175/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:374412/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1981513/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:386548/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1785119/1‑57 (MQ=255)
                               tataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCAACACTTTGCCGTAAt                                     >  1:2074710/1‑57 (MQ=255)
                               tataGCAGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:441482/1‑57 (MQ=255)
                               cataGCCGCCCTGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAAt                                     >  1:1678164/2‑57 (MQ=255)
                                 taGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATc                        >  1:2815349/1‑68 (MQ=255)
                                 taGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATc                        >  1:2872531/1‑68 (MQ=255)
                                 taGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATc                        >  1:2963613/1‑68 (MQ=255)
                                 taGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATc                        >  1:1026802/1‑68 (MQ=255)
                                 taGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATc                        >  1:2065825/1‑68 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:989715/68‑1 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:1147890/68‑1 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:421521/68‑1 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:1377398/68‑1 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:2409931/68‑1 (MQ=255)
                                     cgccCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGaaa                    <  1:648187/68‑1 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAATCgcagt              <  1:435674/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:3190573/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:3222019/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:363517/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:1591817/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:2886632/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:2500590/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:2841348/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:257597/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:1221850/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:2690151/69‑2 (MQ=255)
                                          tGATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:115353/69‑2 (MQ=255)
                                           gATGCTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:2028055/68‑2 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTTACCGCATCGAAACGCAGCACGa          <  1:2923500/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:1521407/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2725881/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:984644/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:948339/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:1011445/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2758211/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2593326/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:85522/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:253578/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:735862/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:698369/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:3210389/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:1920451/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:226443/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:3167717/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2297863/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2203127/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2134470/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:414288/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:1944242/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2026746/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:2695694/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:3159257/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc            <  1:3395625/67‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:1665726/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:960785/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:1129309/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:631651/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:1400173/69‑1 (MQ=255)
                                              gCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:2461060/69‑1 (MQ=255)
                                               cTCCATCACTTCCTGATTATTTACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagt              <  1:57724/64‑2 (MQ=255)
                                               cTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:2742613/68‑1 (MQ=255)
                                               cTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGa          <  1:3004882/68‑1 (MQ=255)
                                                      aCTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTgg  <  1:2444329/69‑1 (MQ=255)
                                                                    |                                                      
ACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGG  >  minE/1023398‑1023520

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: