Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,509,148 C→T 66.3% L738L (CTG→TTG)  yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,509,1480CT66.3% 60.8 / 44.9 83L738L (CTG→TTG) yicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base T (29/26);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGCT  >  minE/2509085‑2509212
                                                               |                                                                
gcaCGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTgcgc                                                             <  1:1504965/69‑1 (MQ=255)
 caCGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCa                                                            <  1:1849430/69‑1 (MQ=255)
  aCGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCaa                                                           <  1:2505033/69‑1 (MQ=255)
  aCGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCaa                                                           <  1:3097848/69‑1 (MQ=255)
  aCGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCaa                                                           <  1:3053889/69‑1 (MQ=255)
   cGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCaa                                                           <  1:1481837/68‑1 (MQ=255)
   cGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAAt                                                          >  1:2037247/1‑69 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2734777/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:1080776/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2035368/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2738803/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2070728/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:289242/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2197650/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:1044409/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2278484/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:2474897/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:874757/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:1092534/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:3326869/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:853037/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:3422663/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:316369/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:600391/67‑1 (MQ=255)
            aaCACGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAATATTGTGa                                                   <  1:595979/67‑1 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:398036/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2674164/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:72634/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:702073/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:620845/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:409378/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2846957/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2995227/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2998929/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:3049378/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:3253052/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:3101167/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:3156095/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2303844/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1266846/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1265991/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1264489/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1399222/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1521006/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:157086/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1579935/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1612773/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1705050/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1765288/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:120970/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1161699/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:183170/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:2267496/1‑62 (MQ=255)
              cacGATTACTGTGACTGGTACGGGCGAGGCGAAGAACTGGACACTGTGCTTGCGCAAtattg                                                      >  1:1854833/1‑62 (MQ=255)
               acGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTa                                              <  1:3012197/69‑1 (MQ=255)
               acGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTa                                              <  1:1735739/69‑1 (MQ=255)
                  aTTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCaa                                                           <  1:3133712/53‑1 (MQ=255)
                       tgtgACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTc                                       <  1:3201124/68‑1 (MQ=255)
                                        aggcgaAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTAAATgg                                         <  1:2756211/49‑1 (MQ=255)
                                           cgaAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTc                                       <  1:866811/48‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:689126/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:644015/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:64027/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:1138494/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:769266/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:3206983/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:3226738/69‑1 (MQ=255)
                                                    ggACACTGTGCTTGCGCAATATTGTGAAAGTAAATGGTCTGCAAGGCGGTTCGCAGGCTGAAAGTGAGc         <  1:2852436/69‑1 (MQ=255)
                                                        acTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAggg     <  1:1340436/69‑1 (MQ=255)
                                                        acTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAggg     <  1:1387206/69‑1 (MQ=255)
                                                        acTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAggg     <  1:65572/69‑1 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGt                            >  1:2436586/1‑44 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:1353358/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:980953/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:1083494/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:2619424/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:2793972/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:756589/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:729533/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGc   >  1:1796415/1‑69 (MQ=255)
                                                          tgtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGAGAGCAgggg    >  1:181539/1‑68 (MQ=255)
                                                           gtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGCt  >  1:2068110/1‑69 (MQ=255)
                                                           gtgCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGCt  >  1:1212580/1‑69 (MQ=255)
                                                               |                                                                
GCACGTACTGGCAACACGATTACTGTGACTGGTGCGGGCGAGGCGAAGAACTGGACACTGTGCCTGCGCAATGTTGTGAAAGTAAATGGTCTGCAAGACGGTTCGCAGGCTGAAAGTGAGCAGGGGCT  >  minE/2509085‑2509212

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: