Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 443,261 G→A 24.8% G167G (GGC→GGT speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE443,2610GA24.8% 68.2 / 16.7 56G167G (GGC→GGTspeFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base G (14/28);  new base A (6/8);  total (20/36)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAA  >  minE/443195‑443326
                                                                  |                                                                 
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:951395/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:3036427/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:1188877/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:732949/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:2589691/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:3334457/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:1468049/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGAAGATCGcc                                                                 <  1:1335858/69‑1 (MQ=255)
ggCATTAAACACTTTTGCCGCATGTAGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                                 <  1:844562/69‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:2229102/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:2089219/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:2517627/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:1841596/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:1765438/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:1461687/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:1458228/68‑1 (MQ=255)
  cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                                                <  1:2981123/68‑1 (MQ=255)
      aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                                           >  1:360799/1‑69 (MQ=255)
      aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                                           >  1:1457489/1‑69 (MQ=255)
      aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                                           >  1:787286/1‑69 (MQ=255)
      aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                                           >  1:2932650/1‑69 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:2602378/1‑68 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:2524431/1‑68 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:3301183/1‑68 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:352689/1‑68 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:1189236/1‑68 (MQ=255)
       aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca                                                          >  1:1318859/1‑68 (MQ=255)
         cacTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTAc                                                        <  1:2519555/69‑1 (MQ=255)
           cTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGt                                                      >  1:2940839/1‑69 (MQ=255)
           cTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGt                                                      >  1:2286845/1‑69 (MQ=255)
                                        ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc                         >  1:1107246/1‑69 (MQ=255)
                                        ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc                         >  1:3417022/1‑69 (MQ=255)
                                        ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc                         >  1:2963793/1‑69 (MQ=255)
                                        ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc                         >  1:3161415/1‑69 (MQ=255)
                                        ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc                         >  1:317494/1‑69 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:1852498/69‑1 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:1830118/69‑1 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:361372/69‑1 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:581962/69‑1 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:916189/69‑1 (MQ=255)
                                           gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc                      <  1:3301474/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:2558143/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:3056255/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:2874992/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:2669583/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:1768381/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:1764390/69‑1 (MQ=255)
                                              ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:1122716/69‑1 (MQ=255)
                                               cTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:1183350/68‑1 (MQ=255)
                                                         cagcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa                   <  1:3371298/58‑1 (MQ=255)
                                                                tCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACaaa  >  1:3334999/1‑68 (MQ=255)
                                                                tCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACaaa  >  1:1794275/1‑68 (MQ=255)
                                                                tCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACaaa  >  1:145008/1‑68 (MQ=255)
                                                                 cGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAg             <  1:1911529/56‑1 (MQ=255)
                                                                 cGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAg             <  1:1400702/56‑1 (MQ=255)
                                                                 cGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAg             <  1:1161933/56‑1 (MQ=255)
                                                                  |                                                                 
GGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAA  >  minE/443195‑443326

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: