Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 133,667 G→A 22.8% H225H (CAC→CAT panB ← 3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE133,6670GA22.8% 99.0 / 21.2 74H225H (CAC→CATpanB3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (35/22);  new base A (11/6);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ATACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCC  >  minE/133601‑133735
                                                                  |                                                                    
ataCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    <  1:172023/69‑1 (MQ=255)
ataCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    <  1:2387989/69‑1 (MQ=255)
ataCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    <  1:3376012/69‑1 (MQ=255)
ataCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    <  1:1072764/69‑1 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGGGAATTTAGGAATGTGa                                                                   >  1:2310827/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtgc                                                                   >  1:2038735/1‑68 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    >  1:3424369/1‑68 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgt                                                                     <  1:2799401/67‑1 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgt                                                                     <  1:2572229/67‑1 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:228127/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2279654/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2087705/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2402522/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2420758/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2454259/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2698521/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2478893/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2569973/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:265410/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1016487/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2710593/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2885652/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2899960/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:294779/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:3019070/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:421463/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:587515/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:670158/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:685255/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:77487/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1046338/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1486506/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1438063/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1415267/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1755591/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1405073/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1372713/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1138996/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1936770/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2010848/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:1122291/1‑69 (MQ=255)
 taCTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:2052771/1‑69 (MQ=255)
 taCTGCCGCACAGCCACGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGa                                                                   >  1:377213/1‑69 (MQ=255)
         cacaGCCGCGCGTATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:831421/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:1025143/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:2838098/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:1083309/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:1459753/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:159999/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:2625244/69‑1 (MQ=255)
         cacaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:1836560/69‑1 (MQ=255)
          acaGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTa                                                          <  1:2558379/69‑1 (MQ=255)
                  gcgGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCaaa                                                  <  1:626881/69‑1 (MQ=255)
                  gcgGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCaaa                                                  <  1:191674/69‑1 (MQ=255)
                   cgGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCaaa                                                  <  1:2073569/68‑1 (MQ=255)
                     gATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCaaa                                                  <  1:1396016/66‑1 (MQ=255)
                      aTGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAAtgtg                                                                    <  1:2880742/47‑1 (MQ=255)
                         tCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:327985/53‑1 (MQ=255)
                         tCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:600506/53‑1 (MQ=255)
                         tCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:2032742/53‑1 (MQ=255)
                         tCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:1680878/53‑1 (MQ=255)
                         tCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:381288/53‑1 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGATTTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:3106266/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAAt                                                        >  1:2822524/1‑52 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:37291/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:2251736/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:2313557/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:3235588/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:3185346/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:1005951/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:1179660/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:1905814/1‑68 (MQ=255)
                             ccGTTTCGGCGAGGAAATTTTTGGCGAATTTAGGAATATGACCGCCGGTAATGCCGAACGCGTCGTGc                                        >  1:2628821/1‑68 (MQ=255)
                                     gCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGt                                                           <  1:3249705/41‑1 (MQ=255)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:1061491/38‑2 (MQ=38)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:905682/38‑2 (MQ=38)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:368518/38‑2 (MQ=38)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:110853/38‑2 (MQ=38)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:2970030/38‑2 (MQ=38)
                                            aaTTTTTGGCGAATTTAGGAATATGACCGCCGGTAAtg                                                       <  1:1943256/38‑2 (MQ=38)
                                                          ttAGGAATGTGAGCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTg          <  1:78467/69‑1 (MQ=255)
                                                          ttAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTg          <  1:2320297/69‑1 (MQ=255)
                                                          ttAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTg          <  1:2053693/69‑1 (MQ=255)
                                                                  gtgACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgcc  <  1:404126/69‑1 (MQ=255)
                                                                  gtgACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgcc  <  1:3394497/69‑1 (MQ=255)
                                                                  gtgACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgcc  <  1:3055031/69‑1 (MQ=255)
                                                                  gtgACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgcc  <  1:2762769/69‑1 (MQ=255)
                                                                  gtgACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCgcc  <  1:840129/69‑1 (MQ=255)
                                                                  |                                                                    
ATACTGCCGCACAGCCGCGCGGATGTCGCCCGTTTCGGCGAGGAAATTTTTAGCGAATTTAGGAATGTGACCGCCGGTAATACCAAAGGCGTCGTGCATCACGAGGATCTGCCCGTCAGTGACGTTGCCTGCGCC  >  minE/133601‑133735

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: