Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 224,081 G→T 44.2% A338A (GCC→GCA pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE224,0810GT44.2% 38.7 / 71.4 79A338A (GCC→GCApepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (17/27);  new base T (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

TGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGA  >  minE/224013‑224147
                                                                    |                                                                  
tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                                    <  1:2887482/69‑1 (MQ=255)
tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                                    <  1:344637/69‑1 (MQ=255)
tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                                    <  1:1739438/69‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:1585655/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2123258/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:230761/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2309777/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2343708/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:1971727/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2839255/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2902742/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:294265/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:3061635/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:368695/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:132630/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:853075/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:883151/66‑1 (MQ=255)
            ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:939391/66‑1 (MQ=255)
             ttttGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:3202044/63‑1 (MQ=255)
             tATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                                           <  1:2220054/65‑1 (MQ=255)
              aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:1638182/1‑69 (MQ=255)
              aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:1491863/1‑69 (MQ=255)
              aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                                      >  1:991492/1‑69 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:2610103/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:3197014/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:804585/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:2147632/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1502096/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1522504/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1433762/69‑1 (MQ=255)
                      tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                              <  1:1885387/69‑1 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:3422128/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:972795/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:1988690/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:3115915/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:2842390/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:596874/1‑69 (MQ=255)
                       caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                             >  1:1291840/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:409882/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2036362/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1493983/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2011496/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2336890/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:198830/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2981578/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:207497/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1439515/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:1902245/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:2509527/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:701470/1‑69 (MQ=255)
                                  cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                                  >  1:76017/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:840988/1‑68 (MQ=255)
                                                gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag                     >  1:2969346/1‑68 (MQ=255)
                                                gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag                     >  1:2556935/1‑68 (MQ=255)
                                                     ttCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGaa                >  1:2592703/1‑68 (MQ=255)
                                                      tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:1168803/69‑1 (MQ=255)
                                                      tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa              <  1:2962544/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:1165435/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2265793/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1168352/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2471120/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2603108/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:723280/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:695261/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2725793/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2815058/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1866948/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:388463/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:2921829/69‑1 (MQ=255)
                                                                  ttGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCga  <  1:1635060/69‑1 (MQ=255)
                                                                    |                                                                  
TGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGA  >  minE/224013‑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: