Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,337,984 C→T 22.8% S91S (TCG→TCA wcaE ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,337,9840CT22.8% 97.8 / 21.4 74S91S (TCG→TCAwcaEpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (20/37);  new base T (2/15);  total (22/52)
Fisher's exact test for biased strand distribution p-value = 7.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ACGTTATCTTTTTGCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTA  >  minE/1337919‑1338046
                                                                 |                                                              
aCGTTATCTTTTTGCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGaa                                                             >  1:2262951/1‑69 (MQ=255)
aCGTTATCTTTTTGCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGaa                                                             >  1:1055983/1‑69 (MQ=255)
             gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                 >  1:1577171/1‑68 (MQ=255)
             gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                 >  1:1386046/1‑68 (MQ=255)
             gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                 >  1:517398/1‑68 (MQ=255)
             gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                 >  1:386649/1‑68 (MQ=255)
             gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaa                                                  >  1:306066/1‑67 (MQ=255)
                 ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:2174220/69‑1 (MQ=255)
                 ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:2447819/69‑1 (MQ=255)
                 ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:2851986/69‑1 (MQ=255)
                 ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:668319/69‑1 (MQ=255)
                 ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:1091059/69‑1 (MQ=255)
                    tAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:802628/66‑1 (MQ=255)
                      aCTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTg                                        <  1:1313588/68‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2771230/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2505998/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2782034/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2756605/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2746826/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:3087067/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:2180335/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:1023828/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:3096682/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:1793977/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:1499250/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:1437758/69‑1 (MQ=255)
                                   aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:142918/69‑1 (MQ=255)
                                    tttGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:1538628/68‑1 (MQ=255)
                                    tttGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                          <  1:664830/68‑1 (MQ=255)
                                             tccTGATGAAAAATATCGCCTGAATTGAGAAACAGCGCGAACt                                          >  1:2422006/3‑43 (MQ=37)
                                             tccTGATGAAAAATATCGACTGAATTGAGAAACAGCGCGAACTTGCCTtg                                   >  1:647211/3‑50 (MQ=255)
                                              tCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                            <  1:2302390/40‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:3099760/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:3150088/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:3221528/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:1050338/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              >  1:1225801/1‑69 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:2743160/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:1768613/69‑1 (MQ=255)
                                               cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:2018488/69‑1 (MQ=255)
                                                tgatgaAAAATATCGCCCGAATTGAGAAACAACGCGAACTTg                                        <  1:3068316/42‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACGGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:660185/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:816265/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:495382/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:1730941/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:627007/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:242891/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:1214916/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:1119482/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:2477654/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:2843292/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:73127/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:2547468/69‑1 (MQ=255)
                                                 gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:2558788/69‑1 (MQ=255)
                                                 gatgaAAAATATAGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:2223816/69‑1 (MQ=255)
                                                  atgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa            <  1:3025518/68‑1 (MQ=255)
                                                           atCGCCCGAATTGAGAGACAACGCGAACTTGCCTTATGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  <  1:523982/69‑1 (MQ=255)
                                                           atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  <  1:743576/69‑1 (MQ=255)
                                                           atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  <  1:3158288/69‑1 (MQ=255)
                                                           atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  <  1:3239828/69‑1 (MQ=255)
                                                           atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTCCCATCGCAATGCCTTTGTTCATGGCGTCGTa  <  1:3379985/69‑1 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:1284568/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:3368545/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:2036926/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:3092651/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:298016/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:1127848/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:2898201/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:699847/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:2836149/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:2478823/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa  >  1:2639101/1‑68 (MQ=255)
                                                            tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTCGTTCATGGCGTCGTa  >  1:1093552/1‑68 (MQ=255)
                                                             cGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt              <  1:3021850/55‑1 (MQ=255)
                                                                 |                                                              
ACGTTATCTTTTTGCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTA  >  minE/1337919‑1338046

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: