Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,629,040 G→T 49.3% R2R (CGC→CGA ypfG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,629,0400GT49.3% 9.6 / 83.0 75R2R (CGC→CGAypfGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/6);  new base T (24/13);  total (56/19)
Fisher's exact test for biased strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCAC  >  minE/1628975‑1629081
                                                                 |                                         
gcTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAgcg                                         <  1:13044/68‑1 (MQ=255)
  tGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:596254/69‑1 (MQ=255)
   ggTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:1288560/69‑1 (MQ=255)
   ggTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATg                                     <  1:3250872/69‑1 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACgc                           >  1:449204/1‑68 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:29625/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1907371/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3102749/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:873835/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1983485/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:892829/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3035912/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3102940/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2882371/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2832538/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2505191/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2725065/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:2717076/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3197360/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1109099/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1142622/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:470750/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:985774/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1382878/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:833812/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3379347/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:1605384/1‑69 (MQ=255)
             acacCAAAGACATCGGCAACAAAACGAAAAAAAAGAGAAGGAAAATGCGATATCGCATGGACGGTACGa                           >  1:3227228/1‑69 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2974952/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:3386579/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2954250/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2800407/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:509999/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2796437/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2722529/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:715590/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2667286/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:1532940/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2207634/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:2161193/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGAcccc                       >  1:1622843/1‑68 (MQ=255)
                  aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGAGGTACGAcccc                       >  1:1952535/1‑67 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1391063/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1323694/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:35977/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1297515/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1604147/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:3347799/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:3196388/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1787434/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:1935759/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:719441/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2248409/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2288851/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2398763/1‑68 (MQ=255)
                           ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTca              >  1:2465133/1‑68 (MQ=255)
                              aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:1673946/69‑1 (MQ=255)
                              aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:1285093/69‑1 (MQ=255)
                              aacaaAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTcaaaca          <  1:2445267/69‑1 (MQ=255)
                                 aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAgcg                                         >  1:3188361/1‑35 (MQ=255)
                                 aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:66064/1‑69 (MQ=255)
                                 aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:2608136/1‑69 (MQ=255)
                                 aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:3050649/1‑69 (MQ=255)
                                 aaaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCaaacaaac       >  1:2581472/1‑69 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:805772/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:898566/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:895447/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:994605/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:655544/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:419012/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:366493/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:3080040/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1671455/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1537532/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1358291/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGcac  <  1:1113967/63‑1 (MQ=255)
                                            aaaGAGAAGGAAAATGCGATATCGCATGGACGGTACAACCCCAGAAGTTCAAACAAACAGcac  <  1:2055067/63‑1 (MQ=255)
                                                                 |                                         
GCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCAC  >  minE/1628975‑1629081

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: