Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 442,546 T→C 73.3% N406D (AAT→GAT)  speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE442,5460TC73.3% 66.0 / 20.6 60N406D (AAT→GAT) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (4/12);  new base C (14/29);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCATTTACCACACAATCCATCCACATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCGTTGTTC  >  minE/442479‑442607
                                                                   |                                                             
ccATTTACCACACAATCCATCCACATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATTc                                                              >  1:1538987/1‑69 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTTATATCCAGTGCGGCAAACAGCgg                                             <  1:156933/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1749499/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2625014/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2613241/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2764650/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2165857/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:216416/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:340393/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:3431484/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1969440/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2776968/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:485245/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1888121/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:2942192/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1827818/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:732016/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:3121905/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1098148/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:11172/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1309735/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1351515/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1743064/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:641158/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:619917/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1573029/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:16623/69‑1 (MQ=255)
                 catccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1723208/69‑1 (MQ=255)
                  atccaCATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:351809/68‑1 (MQ=255)
                            tACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACgcgc                                  >  1:2596837/1‑69 (MQ=255)
                                   gCTGACACCTTCATGCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCgg                                             <  1:1465142/51‑1 (MQ=255)
                                         acCTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTg                     >  1:3152856/1‑69 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:298715/69‑1 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:1839306/69‑1 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:3306991/69‑1 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:424695/69‑1 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:748965/69‑1 (MQ=255)
                                           cTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGca                   <  1:1113668/69‑1 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTg                     >  1:764278/1‑61 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:37383/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1116127/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:74988/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1387815/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1517652/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:541367/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1929268/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1929354/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:1938241/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:3407937/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:3389568/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:2322489/1‑69 (MQ=255)
                                                 gCATTTTGGCGTTGATATACAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:2295372/1‑69 (MQ=255)
                                                 gCATTTTGGAGTTGATATCCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATaa             >  1:145196/1‑69 (MQ=255)
                                                          cGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttg     <  1:2030221/68‑1 (MQ=255)
                                                          cGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttg     <  1:531703/68‑1 (MQ=255)
                                                          cGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttg     <  1:1521941/68‑1 (MQ=255)
                                                          cGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttg     <  1:3336962/68‑1 (MQ=255)
                                                           gTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttgtt   <  1:1989352/69‑1 (MQ=255)
                                                           gTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCgttgtt   <  1:3389806/69‑1 (MQ=255)
                                                            ttGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCGTTGTTc  >  1:3109899/1‑69 (MQ=255)
                                                                   |                                                             
CCATTTACCACACAATCCATCCACATATTACGACCGCTGACACCTTCATGCATTTTGGCGTTGATATTCAGTGCGGCAAACAGCGGATAGAACGGGCTGGTGGAGGCGTGCATCATAAAGGCGTTGTTC  >  minE/442479‑442607

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: