Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,953,382 G→A 46.9% L243L (CTG→TTG)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,953,3820GA46.9% 8.3 / 68.5 81L243L (CTG→TTG) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base G (35/8);  new base A (20/18);  total (55/26)
Fisher's exact test for biased strand distribution p-value = 8.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCG  >  minE/1953314‑1953445
                                                                    |                                                               
ttCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAg                                                                 <  1:1451152/69‑1 (MQ=255)
ttCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAg                                                                 <  1:1828250/69‑1 (MQ=255)
ttCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAg                                                                 <  1:2308587/69‑1 (MQ=255)
  cacaACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGaa                                                               >  1:2968709/1‑69 (MQ=255)
  cacaACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGaa                                                               >  1:1453883/1‑69 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:1539123/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:1987675/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:20576/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:2132997/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:714072/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:2538073/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:595123/69‑1 (MQ=255)
              caatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGtttt                                                   <  1:1114793/69‑1 (MQ=255)
               aatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCg                                                       >  1:1800486/1‑64 (MQ=255)
               aatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTc                                                  >  1:56279/1‑69 (MQ=255)
               aatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTc                                                  >  1:289378/1‑69 (MQ=255)
               aatcaaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTc                                                  >  1:1007134/1‑69 (MQ=255)
                   aaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCAcc                                               >  1:2147907/1‑68 (MQ=255)
                   aaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCAcc                                               >  1:2700247/1‑68 (MQ=255)
                   aaCACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCAc                                                >  1:2836099/1‑67 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:2407212/1‑69 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:1238949/1‑69 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:2918660/1‑69 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:1413795/1‑69 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:2695352/1‑69 (MQ=255)
                     cacCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCt                                            >  1:3349296/1‑69 (MQ=255)
                         aGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGcc                                         >  1:433681/1‑68 (MQ=255)
                         aGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGcc                                         >  1:1265779/1‑68 (MQ=255)
                         aGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGcc                                         >  1:1628064/1‑68 (MQ=255)
                         aGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGcc                                         >  1:2541141/1‑68 (MQ=255)
                          ggTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAg                                                                 <  1:1998373/43‑1 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAgg                        >  1:1975610/1‑67 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGa                       >  1:3176172/1‑68 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGa                       >  1:174312/1‑68 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGa                       >  1:955219/1‑68 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGa                       >  1:2717005/1‑68 (MQ=255)
                                           tCACACTCGTTTTCGGTGTAGTGCAGAATGAAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGa                       >  1:1930093/1‑68 (MQ=255)
                                                   gTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAgcg              >  1:819965/1‑69 (MQ=255)
                                                   gTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAgcg              >  1:2059373/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCcaacaa                                    >  1:2955493/1‑45 (MQ=38)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1131719/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:3405945/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:3018495/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1612396/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:933391/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:257655/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1176788/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:248489/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1723330/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1441913/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:742529/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1931780/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:820702/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:2649791/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:1981111/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgc             >  1:1393170/1‑68 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacg              >  1:27325/1‑67 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGCATCCGTTTTCACAGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:3113413/1‑69 (MQ=255)
                                                     tttCGGTGTAGTGCAAAATGAATCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAacgcg            >  1:3267320/1‑69 (MQ=255)
                                                      ttCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAgcgcgc           >  1:749506/1‑69 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:9048/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:998633/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:636534/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:609753/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:3274039/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2972208/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:1143734/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2869560/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:1214321/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2740246/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2677608/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2389771/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2323047/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:2203246/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:1328254/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:1848473/69‑1 (MQ=255)
                                                           tgtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:154025/69‑1 (MQ=255)
                                                            gtAGTGCAAAATGCAGCCGTTTTCACCGCTGCCAACAATGGTGTTATAGGAAGGGTAACGCGCACCGt      <  1:1056719/68‑1 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGc   >  1:8173/1‑68 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:1759660/1‑69 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:725433/1‑69 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:520699/1‑69 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:3798/1‑69 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:1115436/1‑69 (MQ=255)
                                                               gTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCg  >  1:3251082/1‑69 (MQ=255)
                                                                    |                                                               
TTCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCG  >  minE/1953314‑1953445

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: