Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 647,950 A→G 45.6% A210A (GCA→GCG pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE647,9500AG45.6% 5.2 / 45.9 57A210A (GCA→GCGpepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (13/13);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAC  >  minE/647885‑648016
                                                                 |                                                                  
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:2864050/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:597660/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1861711/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1761963/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:3386729/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1537979/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:701460/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:1297621/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:2906822/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:742985/67‑1 (MQ=255)
ttggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTAACCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:275477/67‑1 (MQ=255)
gtggtggCTGGTGATTTTGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:342843/67‑1 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTgg                                                                >  1:2622762/1‑68 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:258749/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:2484557/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:343062/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:1102925/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:2030282/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:2010221/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:451210/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:1248359/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:1467933/1‑69 (MQ=255)
  ggtggCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGa                                                               >  1:1650299/1‑69 (MQ=255)
             attttGATGTACTGCGCGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCt                                                                  <  1:756278/53‑1 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:448900/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:707418/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:2915786/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1549525/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:243221/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:3330291/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:2323823/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1668907/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1712456/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1869964/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:1944714/1‑68 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTCTGGCCGTGAACTGGCGCTGGAGTTGTATGTCGATCGTGGCaa                                          >  1:2596225/1‑66 (MQ=255)
                          gcgcGATACCTTCACTACCCGCTATGGCCGTGAAGTGGCGCTGGAGTTGTATGTCGATCGTGGCAAcc                                        >  1:2069169/1‑68 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCTCGGCAACCtt                                      <  1:236491/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1898983/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:986495/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1132963/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1383609/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1470369/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:612270/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1867879/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2703997/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:391927/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2214016/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:2219895/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:3387716/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCtt                                      <  1:1083147/63‑1 (MQ=255)
                                 aCCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGGCGATCGCGGCAACCtt                                      <  1:2221256/63‑1 (MQ=255)
                                                 cTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGa                >  1:2510132/1‑69 (MQ=255)
                                                           aGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaa      <  1:197465/69‑1 (MQ=255)
                                                            gTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaaa     <  1:1191898/69‑1 (MQ=255)
                                                             tAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGaaa     <  1:721454/68‑1 (MQ=255)
                                                               gCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAc  >  1:858410/1‑69 (MQ=255)
                                                                 |                                                                  
CTGGTGGCAGGCGACTTTGATGTACTGCGCGATACCTTTACCACGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAAC  >  minE/647885‑648016

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: