Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,005,190 G→A 32.1% N135N (AAC→AAT pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,005,1900GA32.1% 76.9 / 43.1 87N135N (AAC→AATpntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (18/41);  new base A (22/6);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 3.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGC  >  minE/1005125‑1005238
                                                                 |                                                
aTTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTc                                               >  1:300208/1‑69 (MQ=255)
aTTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTc                                               >  1:2077038/1‑69 (MQ=255)
aTTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTc                                               >  1:2446057/1‑69 (MQ=255)
  ttGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGc                                             <  1:2462719/69‑1 (MQ=255)
  ttGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGc                                             <  1:3432698/69‑1 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGcc                                            >  1:906170/1‑69 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGcc                                            >  1:298412/1‑69 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGcc                                            >  1:2696105/1‑69 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGcc                                            >  1:2680242/1‑69 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGcc                                            >  1:2204175/1‑69 (MQ=255)
   tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGc                                             >  1:2848273/1‑68 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCa                                           >  1:1624577/1‑62 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAg                                      >  1:2660288/1‑67 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAg                                      >  1:1056274/1‑67 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAg                                      >  1:2404956/1‑67 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAg                                      >  1:794419/1‑67 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAg                                      >  1:1734478/1‑67 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:2274202/69‑1 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:191143/69‑1 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:1796381/69‑1 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:1393198/69‑1 (MQ=255)
           aaaGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:725214/69‑1 (MQ=255)
            aagaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCtt                                   <  1:309007/69‑1 (MQ=255)
            aagaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCt                                    <  1:2794472/68‑1 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:2834042/1‑69 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:2056257/1‑69 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:1585948/1‑69 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:3164750/1‑69 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:212901/1‑69 (MQ=255)
             agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTa                                  >  1:1188059/1‑69 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:1738027/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:2228729/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:531238/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:2472177/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:541037/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:353794/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:2938090/69‑1 (MQ=255)
                         aTTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTACGTCCAGCg                      <  1:1940453/69‑1 (MQ=255)
                          ttCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCg                      <  1:2578804/68‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:2205871/69‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:555458/69‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:3376590/69‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:1140180/69‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:2816637/69‑1 (MQ=255)
                           tCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:2820674/69‑1 (MQ=255)
                            cATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAtt                   >  1:2412029/1‑69 (MQ=255)
                            cATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                    <  1:2684363/68‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:138815/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:452579/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:3313042/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:3257720/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:702933/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:1612221/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:2431537/69‑1 (MQ=255)
                                  ccGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:208897/69‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:1817727/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:1714366/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:1639762/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:2829099/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:543056/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:2270897/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:613629/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:329080/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:1482921/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:334053/67‑1 (MQ=255)
                                          acaaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:2135637/67‑1 (MQ=255)
                                           caaTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:1653243/66‑1 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGAtt                   >  1:1744825/1‑53 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:577205/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:776197/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:622967/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:726596/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:1180428/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:1065958/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:1355164/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:444372/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:334847/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:1504576/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:3132635/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:303456/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:1981900/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:2094517/1‑69 (MQ=255)
                                            aaTGGCGCGATAACCGGCGATATTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg   >  1:2178057/1‑69 (MQ=255)
                                             aTGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgc  >  1:1404169/1‑69 (MQ=255)
                                             aTGGCGCGATAACCGGCGATGTTCACCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgc  >  1:2134559/1‑69 (MQ=255)
                                                  gcgATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc             <  1:733974/53‑1 (MQ=255)
                                                             cGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga       <  1:325161/48‑1 (MQ=255)
                                                                 |                                                
ATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGC  >  minE/1005125‑1005238

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: