Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 443,246 G→A 22.6% H172H (CAC→CAT speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE443,2460GA22.6% 89.5 / 17.8 66H172H (CAC→CATspeFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base G (16/35);  new base A (7/8);  total (23/43)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

TTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCA  >  minE/443188‑443309
                                                          |                                                               
ttttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                                                       <  1:3270332/69‑1 (MQ=255)
ttttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                                                       <  1:993111/69‑1 (MQ=255)
ttttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                                                       <  1:2003341/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGTATcagcag                                                      <  1:2117492/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:1027132/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:3093638/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:2895602/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:2691573/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:629468/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:1373015/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                                                      <  1:112846/69‑1 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                                                       >  1:221263/1‑68 (MQ=255)
 tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATc                                                           >  1:2345477/1‑64 (MQ=255)
  ttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATc                                                           >  1:1429592/1‑63 (MQ=255)
  ttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATc                                                           >  1:2794263/1‑63 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:1188877/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:951395/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:732949/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:2589691/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:1468049/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:3334457/69‑1 (MQ=255)
       ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc                                                <  1:3036427/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:2517627/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:1458228/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:2981123/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:1461687/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:2089219/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:1841596/68‑1 (MQ=255)
         cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc                                               <  1:1765438/68‑1 (MQ=255)
             aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                          >  1:360799/1‑69 (MQ=255)
             aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                          >  1:2932650/1‑69 (MQ=255)
             aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                          >  1:787286/1‑69 (MQ=255)
             aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc                                          >  1:1457489/1‑69 (MQ=255)
              aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATcagaa                                                       >  1:1982993/1‑53 (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:2602378/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)
                  cTTTTGCCGCATGTTGCTGTGCAATGCATGTCGCGCCTTCGTGAATc                                                           >  1:1257114/1‑47 (MQ=255)
                  cTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGt                                     >  1:2940839/1‑69 (MQ=255)
                  cTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGt                                     >  1:2286845/1‑69 (MQ=255)
                                               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc        >  1:3417022/1‑69 (MQ=255)
                                               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc        >  1:1107246/1‑69 (MQ=255)
                                               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc        >  1:3161415/1‑69 (MQ=255)
                                               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc        >  1:2963793/1‑69 (MQ=255)
                                               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAgc        >  1:317494/1‑69 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:361372/69‑1 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:916189/69‑1 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:1830118/69‑1 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:1852498/69‑1 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:581962/69‑1 (MQ=255)
                                                  gcgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGcc     <  1:3301474/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:2558143/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:2669583/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:2874992/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:1768381/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:1764390/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:3056255/69‑1 (MQ=255)
                                                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:1122716/69‑1 (MQ=255)
                                                      cTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCa  <  1:1183350/68‑1 (MQ=255)
                                                          |                                                               
TTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCA  >  minE/443188‑443309

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: