Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 224,087 A→G 44.3% D336D (GAT→GAC pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE224,0870AG44.3% 13.7 / 63.7 79D336D (GAT→GACpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (17/27);  new base G (17/18);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGA  >  minE/224025‑224147
                                                              |                                                            
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:853075/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:3061635/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:294265/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2902742/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:1971727/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2123258/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:883151/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:939391/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:1585655/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:230761/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2309777/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:132630/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2343708/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:368695/66‑1 (MQ=255)
ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2839255/66‑1 (MQ=255)
 ttttGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:3202044/63‑1 (MQ=255)
 tATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2220054/65‑1 (MQ=255)
  aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:1491863/1‑69 (MQ=255)
  aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:1638182/1‑69 (MQ=255)
  aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:991492/1‑69 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:2610103/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:2147632/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:804585/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:3197014/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1885387/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1522504/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1502096/69‑1 (MQ=255)
          tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1433762/69‑1 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:1988690/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:3115915/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:3422128/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:2842390/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:1291840/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:596874/1‑69 (MQ=255)
           caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:972795/1‑69 (MQ=255)
                    cacaCCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCGACGTCg                                                           <  1:613716/46‑1 (MQ=37)
                      cacCGACGTTCAGGGAGGTTTCGCCCACACCTTTTGCCACGTCGGAGTTGAGGATCACACCGTTCgggg                                  >  1:2797708/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGTTCACACCGTTCgggg                                  >  1:2388271/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2336890/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1493983/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1902245/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:198830/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2011496/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2981578/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2509527/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2036362/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1439515/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:207497/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:701470/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:76017/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:409882/1‑69 (MQ=255)
                      cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCggg                                   >  1:2968440/1‑68 (MQ=255)
                      cacCGAAGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:439860/1‑69 (MQ=255)
                                    gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCgg                                    >  1:2384751/1‑53 (MQ=255)
                                    gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag                     >  1:2969346/1‑68 (MQ=255)
                                    gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag                     >  1:840988/1‑68 (MQ=255)
                                    gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag                     >  1:2556935/1‑68 (MQ=255)
                                         ttCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGaa                >  1:2592703/1‑68 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:2962544/69‑1 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:1168803/69‑1 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:3255515/69‑1 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:2832539/69‑1 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:2700986/69‑1 (MQ=255)
                                          tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:2691994/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATc      >  1:1269815/1‑65 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATc      >  1:111847/1‑65 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:723280/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2603108/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:695261/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:388463/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2921829/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2815058/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2725793/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2471120/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2265793/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1866948/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1635060/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1168352/69‑1 (MQ=255)
                                                      ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1165435/69‑1 (MQ=255)
                                                              |                                                            
TTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGA  >  minE/224025‑224147

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: