Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,108,643 G→C 36.5% L138L (CTG→CTC kbaY → tagatose 6‑phosphate aldolase 1, kbaY subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,108,6430GC36.5% 59.2 / 65.1 78L138L (CTG→CTCkbaYtagatose 6‑phosphate aldolase 1, kbaY subunit
Reads supporting (aligned to +/- strand):  ref base G (17/30);  new base C (8/20);  total (25/53)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTCCTGACCGATCCACAAGAAGCTA  >  minE/2108575‑2108710
                                                                    |                                                                   
cTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTg                                                                     <  1:1132669/69‑1 (MQ=255)
cTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTg                                                                     <  1:103442/69‑1 (MQ=255)
cTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTg                                                                     <  1:1337947/69‑1 (MQ=255)
cTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTg                                                                     <  1:2138809/69‑1 (MQ=255)
     gAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCCCTGGGCGgt                                                               <  1:1210538/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTTGGCGgt                                                               <  1:440964/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:303378/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:463472/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:330574/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:3303127/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:2162041/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:2831868/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:1981984/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:1922503/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:3095063/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:1568452/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               >  1:1558221/1‑69 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:1470811/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGgt                                                               <  1:1187335/69‑1 (MQ=255)
      aaaTCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTAGGCGgt                                                               <  1:1979680/69‑1 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:2470910/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:2404850/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:2846449/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:2991590/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:2255627/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:90784/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:1254521/1‑68 (MQ=255)
        aTCGGTTGTCGACTTCTGCCACTCCCAGGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTCGGCGgtg                                                              >  1:1352660/1‑68 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGaa                                                         <  1:1575329/68‑1 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGaa                                                         <  1:1084740/68‑1 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGaa                                                         <  1:559259/68‑1 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGaa                                                         <  1:609833/68‑1 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGaa                                                         <  1:1565135/68‑1 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAg                                                        >  1:347753/1‑69 (MQ=255)
             ttgttgACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAg                                                        >  1:2426429/1‑69 (MQ=255)
                  gACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    <  1:1269058/68‑1 (MQ=255)
                  gACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    <  1:2006763/68‑1 (MQ=255)
                  gACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    <  1:755041/68‑1 (MQ=255)
                  gACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    <  1:818254/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:2834508/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:1064556/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:247114/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:2393023/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:506496/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:2225265/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:2153834/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:178902/68‑1 (MQ=255)
                            aCTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGtt                                          <  1:1526320/68‑1 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:3285058/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:3261625/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:625085/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:861391/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:953462/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:1177718/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:2638575/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:2875672/1‑55 (MQ=255)
                               cacaAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGa                                                    >  1:2060965/1‑55 (MQ=255)
                                        gCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTg                             <  1:748033/69‑1 (MQ=255)
                                        gCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTg                             <  1:395746/69‑1 (MQ=255)
                                                gAAGCTGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:3374551/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:489854/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:1204972/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:872635/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:1663172/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:1962695/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:660451/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2033761/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2053746/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:584933/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2282237/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:488612/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:267498/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2738266/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2763336/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2833854/67‑1 (MQ=255)
                                                gAAGCAGAACTGGGCCGCCTAGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:434808/67‑1 (MQ=255)
                                                 aaGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:3378705/66‑1 (MQ=255)
                                                 aaGCAGAACTGGGCCGCCTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:1693008/66‑1 (MQ=255)
                                                               cgccTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:235322/52‑1 (MQ=255)
                                                                 ccTCGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTcc                       <  1:2644062/50‑1 (MQ=255)
                                                                   tGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTCCTGACCGATCCACAAGAAGCTa  <  1:3414964/69‑1 (MQ=255)
                                                                   tGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTCCTGACCGATCCACAAGAAGCTa  <  1:915564/69‑1 (MQ=255)
                                                                    |                                                                   
CTGGTGAAATCGGTTGTTGACTTCTGCCACTCACAAGATTGCAGCGTGGAAGCAGAACTGGGCCGCCTGGGCGGTGTTGAAGATGACATGAGCGTTGACGCCGAAAGTGCATTCCTGACCGATCCACAAGAAGCTA  >  minE/2108575‑2108710

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: