Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,953,119 G→A 36.2% G330G (GGC→GGT pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,953,1190GA36.2% 99.7 / 87.1 144G330G (GGC→GGTpepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base G (52/39);  new base A (30/22);  total (82/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 144 total aligned reads displayed.

CCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTC  >  minE/1953054‑1953186
                                                                 |                                                                   
ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:2741059/69‑1 (MQ=255)
ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:743981/69‑1 (MQ=255)
ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:1425469/69‑1 (MQ=255)
ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:1224196/69‑1 (MQ=255)
ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:3401201/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:919845/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:2490287/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:2902680/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:3038430/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:794045/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGt                                                               <  1:1189214/69‑1 (MQ=255)
   tGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACACCACCTTTCAGGATGCCGAGt                                                               <  1:53224/69‑1 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:15224/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:2041129/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:2417483/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:2645216/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:2848845/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:1610809/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:2938112/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:1204715/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:3133328/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGttt                                                             >  1:1276740/1‑69 (MQ=255)
     cATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGtt                                                              >  1:497004/1‑68 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:689720/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:362279/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:882284/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:2691527/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:427499/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1084717/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:2179240/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:2108121/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:2048262/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1993231/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1897817/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:187322/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1297032/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1237882/69‑1 (MQ=255)
        aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1165156/69‑1 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:1400470/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:353272/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:1414939/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:3120681/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:423608/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:1264122/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:1686050/1‑69 (MQ=255)
          agaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTACCa                                                        >  1:841493/1‑69 (MQ=255)
                gACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATACCGAGTTTTAc                                                          <  1:1861778/61‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCATGATTCCGAGTTTTACCAGGCCGCTa                                                <  1:2378254/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:227313/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:2341720/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:491875/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:503709/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:2731206/69‑1 (MQ=255)
                  cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTa                                                <  1:1030170/69‑1 (MQ=255)
                       ggCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCat                                           <  1:1173874/69‑1 (MQ=255)
                       ggCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCat                                           <  1:2631436/69‑1 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCg                                                                  <  1:2872097/40‑1 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCatgatg                                       >  1:2259304/1‑67 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCatgatg                                       >  1:1618185/1‑67 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCatgatg                                       >  1:2562771/1‑67 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCatgatg                                       >  1:2624830/1‑67 (MQ=255)
                             cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCatgatg                                       >  1:2610515/1‑67 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGGTGCGCACCACTTc                            >  1:2583594/1‑69 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACtt                             >  1:66004/1‑68 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTc                            >  1:2977285/1‑69 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTc                            >  1:360692/1‑69 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTc                            >  1:1668675/1‑69 (MQ=255)
                                      cAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTc                            >  1:1824645/1‑69 (MQ=255)
                                         ttCATCAACATCACCTTTCATGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCAc                          <  1:2441530/68‑1 (MQ=255)
                                         ttCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCAc                          <  1:473645/68‑1 (MQ=255)
                                         ttCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCAc                          <  1:2081098/68‑1 (MQ=255)
                                         ttCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCAc                          <  1:3114729/68‑1 (MQ=255)
                                         ttCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCAc                          <  1:1514677/68‑1 (MQ=255)
                                          tcatcaACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCa                        >  1:2766994/1‑69 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:3302567/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:2157377/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:1102428/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:3193146/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:1292119/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:743248/1‑63 (MQ=255)
                                                  aTCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGt                      >  1:3035992/1‑63 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg             >  1:3317443/1‑69 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg             >  1:2329864/1‑69 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg             >  1:2660414/1‑69 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg             >  1:602968/1‑69 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg             >  1:1485661/1‑69 (MQ=255)
                                                     aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACAACTTCACCAGAGACTTCCAg             >  1:487831/1‑69 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCTCACCACTTCACCAGTGACTTcctg             >  1:3093860/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:1241949/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:686200/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:147231/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:773602/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:1850551/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:823809/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:3097238/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:326217/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGCACCACTTCACCAGTGACTTcctg             >  1:3397371/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGCGAACCACTTCACCAGTGACTTcctg             >  1:476724/1‑66 (MQ=255)
                                                      ccTTTCAGGATACCGAGTTTTACCAGACCGCTAACCATGATGAGCACCACTTCACCAGTGACTTcctg             >  1:2958230/1‑66 (MQ=255)
                                                                tGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTc  >  1:92933/1‑69 (MQ=255)
                                                                 |                                                                   
CCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTC  >  minE/1953054‑1953186

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: