Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 443,210 T→A 19.6% A184A (GCA→GCT speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE443,2100TA19.6% 142.3 / 26.1 76A184A (GCA→GCTspeFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (23/38);  new base A (7/8);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

AAACCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGT  >  minE/443143‑443274
                                                                   |                                                                
aaaCCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACtttt                                                                  <  1:148212/68‑1 (MQ=255)
aaaCCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACtttt                                                                  <  1:2205925/68‑1 (MQ=255)
aaaCCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACtttt                                                                  <  1:36240/68‑1 (MQ=255)
aaaCCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACtttt                                                                  <  1:2592721/68‑1 (MQ=255)
     aCTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            >  1:1156555/1‑69 (MQ=255)
     aCTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            >  1:883306/1‑69 (MQ=255)
     aCTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            >  1:829464/1‑69 (MQ=255)
     aCTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            <  1:2495518/69‑1 (MQ=255)
     aCTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            >  1:39510/1‑69 (MQ=255)
     aCTTTGTTAGAAGATGAAGTCCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCa                                                            <  1:3397075/69‑1 (MQ=255)
           ttAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGtt                                                        <  1:3225010/67‑1 (MQ=255)
                agaTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGc                                                 <  1:138042/69‑1 (MQ=255)
                agaTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGc                                                 <  1:528290/69‑1 (MQ=255)
                agaTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGc                                                 <  1:3314136/69‑1 (MQ=255)
                agaTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGc                                                 <  1:295749/69‑1 (MQ=255)
                  aTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCaa                                               <  1:874399/69‑1 (MQ=255)
                  aTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCaa                                               <  1:2278566/69‑1 (MQ=255)
                  aTGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCaa                                               <  1:3188385/69‑1 (MQ=255)
                    gAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATg                                             <  1:2103727/69‑1 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:2411924/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:2013328/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:1490070/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:1741496/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:1617589/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:3256584/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:139451/1‑69 (MQ=255)
                            aTTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGgcgc                                     >  1:1109897/1‑69 (MQ=255)
                              ttAAAACGAAGTAGGTTTTATCGGCATTAAACACtttt                                                                  <  1:3162613/38‑1 (MQ=255)
                                  aaCGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCtt                                <  1:3200825/68‑1 (MQ=255)
                                  aaCGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCtt                                <  1:2019666/68‑1 (MQ=255)
                                  aaCGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCtt                                <  1:2252819/68‑1 (MQ=255)
                                  aaCGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCtt                                <  1:937563/68‑1 (MQ=255)
                                   aCGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCtt                                <  1:2870704/67‑1 (MQ=255)
                                            ggtttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCaa                                               <  1:664347/41‑1 (MQ=255)
                                             ttttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                    <  1:993111/69‑1 (MQ=255)
                                             ttttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagca                    <  1:3270332/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:1373015/69‑1 (MQ=255)
                                              tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                   <  1:629468/69‑1 (MQ=255)
                                              tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                   <  1:2895602/69‑1 (MQ=255)
                                              tttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATcagcag                   <  1:2691573/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:2794263/1‑63 (MQ=255)
                                               ttATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATc                        >  1:1429592/1‑63 (MQ=255)
                                                    ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc             <  1:732949/69‑1 (MQ=255)
                                                    ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc             <  1:2589691/69‑1 (MQ=255)
                                                    ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc             <  1:1188877/69‑1 (MQ=255)
                                                    ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc             <  1:3036427/69‑1 (MQ=255)
                                                    ggCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGcc             <  1:951395/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:1765438/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:1458228/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:1461687/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:1841596/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:2229102/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:2981123/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:2089219/68‑1 (MQ=255)
                                                      cATTAAACACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCAccc            <  1:2517627/68‑1 (MQ=255)
                                                          aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc       >  1:1457489/1‑69 (MQ=255)
                                                          aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc       >  1:360799/1‑69 (MQ=255)
                                                          aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc       >  1:2932650/1‑69 (MQ=255)
                                                          aaaCACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGc       >  1:787286/1‑69 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATcagaa                    >  1:1982993/1‑53 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:3301183/1‑68 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:1189236/1‑68 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:2524431/1‑68 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:2602378/1‑68 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:1318859/1‑68 (MQ=255)
                                                           aaCACTTTAGCCGCATGTTGCTGTGCAATGCATGGCGCACCTTCATGAATCAGAAGATCACCCATTGca      >  1:352689/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)
                                                                   |                                                                
AAACCACTTTGTTAGAAGATGAAGTGCCATTTAAAACGAAGTAGGTTTTATCGGCATTAAACACTTTTGCCGCATGTTGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGT  >  minE/443143‑443274

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: