Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,337,998 A→G 21.1% L87L (TTG→CTG)  wcaE ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,337,9980AG21.1% 108.6 / 21.8 80L87L (TTG→CTG) wcaEpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (24/39);  new base G (2/15);  total (26/54)
Fisher's exact test for biased strand distribution p-value = 4.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

GCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATC  >  minE/1337932‑1338058
                                                                  |                                                            
gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                             >  1:1577171/1‑68 (MQ=255)
gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                             >  1:1386046/1‑68 (MQ=255)
gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                             >  1:386649/1‑68 (MQ=255)
gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaac                                                             >  1:517398/1‑68 (MQ=255)
gCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAaacaa                                                              >  1:306066/1‑67 (MQ=255)
    ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                                        <  1:2447819/69‑1 (MQ=255)
    ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                                        <  1:2851986/69‑1 (MQ=255)
    ttttAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGaa                                                        <  1:2174220/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:2746826/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:2180335/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:2505998/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:2756605/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:2771230/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:2782034/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:1793977/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:1023828/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:3087067/69‑1 (MQ=255)
                      aTTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGc                                      <  1:3096682/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:2743160/69‑1 (MQ=255)
                                  cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt                          <  1:1050338/69‑1 (MQ=255)
                                  cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt                          <  1:1768613/69‑1 (MQ=255)
                                  cTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt                          >  1:1225801/1‑69 (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:627007/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:1730941/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:1119482/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:73127/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:2843292/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:2477654/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:816265/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:495382/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:2558788/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:2547468/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:242891/69‑1 (MQ=255)
                                    gatgaAAAATATCGCCTGAATTGAGAAACAGCGCGAACTTGCCTTGCGCCATTTCAATACCTTTGTTCa                        <  1:1214916/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:3239828/69‑1 (MQ=255)
                                              atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              <  1:743576/69‑1 (MQ=255)
                                              atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              <  1:3158288/69‑1 (MQ=255)
                                              atCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTCCCATCGCAATGCCTTTGTTCATGGCGTCGTa              <  1:3379985/69‑1 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:2478823/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:2639101/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:3368545/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:2036926/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:1284568/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:1127848/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:2836149/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:2898201/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:3092651/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:298016/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTa              >  1:699847/1‑68 (MQ=255)
                                               tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTCGTTCATGGCGTCGTa              >  1:1093552/1‑68 (MQ=255)
                                                cGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGtt                          <  1:3021850/55‑1 (MQ=255)
                                                       aTTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGtt     <  1:1471341/69‑1 (MQ=255)
                                                       aTTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGtt     <  1:2179397/69‑1 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGcgtcg                >  1:2198611/1‑55 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:348371/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:3323245/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:1178783/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:1448057/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:3134361/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:1595428/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:3005423/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:1931970/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:2814834/1‑69 (MQ=255)
                                                          gagaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATc  >  1:920521/1‑69 (MQ=255)
                                                                  |                                                            
GCATTTTTAACTTACGGACAAAATTTGCGGCATTCTGATGAAAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATC  >  minE/1337932‑1338058

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