Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 157,243 A→G 55.3% L15L (CTT→CTC btuF ← vitamin B12 transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE157,2430AG55.3% 13.4 / 51.9 67L15L (CTT→CTCbtuFvitamin B12 transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (18/12);  new base G (13/24);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 5.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGT  >  minE/157177‑157307
                                                                  |                                                                
ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa                                                                 >  1:867947/1‑68 (MQ=255)
ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa                                                                 >  1:1115255/1‑68 (MQ=255)
ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa                                                                 >  1:2025753/1‑68 (MQ=255)
ggCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCaa                                                                 >  1:1681464/1‑68 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:452317/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:2895811/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:2902772/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:1701905/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:1618794/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:280808/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:589132/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:195580/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:1307685/69‑1 (MQ=255)
      ttCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:2242587/69‑1 (MQ=255)
      ttCAGTGTTCGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:1889126/69‑1 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGaaaa                                                            >  1:1486485/1‑65 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:949111/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:1952064/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:2888518/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:2408796/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:3255873/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:865300/1‑69 (MQ=255)
        cAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACa                                                        >  1:1103349/1‑69 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:1789629/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:3056123/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:329390/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:1349291/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:120550/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:2331715/69‑1 (MQ=255)
                            gTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:2272093/69‑1 (MQ=255)
                             tgatgaCGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAAc                                    <  1:2322932/68‑1 (MQ=255)
                              gatgaCGCGCGGCGCGGCATTGAGCCACAGCGGCGCGAGAAAAGa                                                          <  1:811066/45‑1 (MQ=38)
                                       cggcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGc                         <  1:3079822/69‑1 (MQ=255)
                                        ggcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGc                         <  1:2878466/68‑1 (MQ=255)
                                         gcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCa                       <  1:1497878/69‑1 (MQ=255)
                                          cgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCa                       <  1:1081336/68‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2883458/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:455007/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:3233735/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:1742025/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:1425628/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:1113678/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2358906/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2042231/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2190159/63‑1 (MQ=255)
                                               cattGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:21503/63‑1 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2476514/63‑1 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:469923/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:2443549/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:460738/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:755123/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:3337974/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:80751/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:3250186/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:3244723/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:1038039/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:2206324/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:1808901/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:2387336/2‑67 (MQ=255)
                                                attGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGcca                >  1:2041870/2‑67 (MQ=255)
                                                  tGAGCCACAGCGGCGCGAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCAtg                     <  1:2852459/62‑1 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:166052/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:84383/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:528192/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:1463372/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:1492078/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:1702071/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:1931471/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:2723695/1‑69 (MQ=255)
                                                   gAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGtt             >  1:2211823/1‑69 (MQ=255)
                                                              gcgcAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGt  >  1:1458254/1‑69 (MQ=255)
                                                                  aaGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGc                         <  1:2236836/42‑1 (MQ=255)
                                                                  |                                                                
GGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGT  >  minE/157177‑157307

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