Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 224,069 A→G 52.3% V342V (GTT→GTC pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE224,0690AG52.3% 3.1 / 56.7 67V342V (GTT→GTCpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (15/17);  new base G (17/18);  total (32/35)
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.73e-01

CGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATA  >  minE/224001‑224135
                                                                    |                                                                  
cGGATCAGGCAGTGAATTTCTATGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:3186797/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:1190912/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:2505399/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:1453786/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:538026/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:2042717/69‑1 (MQ=255)
cGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCa                                                                    <  1:2794605/69‑1 (MQ=255)
            tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                        <  1:344637/69‑1 (MQ=255)
            tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                        <  1:1739438/69‑1 (MQ=255)
            tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTg                                                        <  1:2887482/69‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:2309777/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:230761/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:2343708/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:2123258/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:2839255/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:2902742/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:1971727/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:294265/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:3061635/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:1585655/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:368695/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:853075/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:883151/66‑1 (MQ=255)
                        ttATTGTCGGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCg                                               <  1:132630/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:991492/1‑69 (MQ=255)
                          aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                          >  1:1491863/1‑69 (MQ=255)
                          aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGtt                                          >  1:1638182/1‑69 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:2147632/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:1433762/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:804585/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:1502096/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:1522504/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:3197014/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:1885387/69‑1 (MQ=255)
                                  tcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATca                                  <  1:2610103/69‑1 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:3422128/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:972795/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:2842390/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:1988690/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:596874/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:1291840/1‑69 (MQ=255)
                                   caTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATcac                                 >  1:3115915/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:1439515/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:207497/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:1493983/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:409882/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:701470/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:76017/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:2509527/1‑69 (MQ=255)
                                              cacCGACGTTCAGGGAGGTTTCGACCACACCTTTTGCCACGTCGGAGTTGCGGATCACACCGTTCgggg                      >  1:2981578/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:2036362/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:2556935/1‑68 (MQ=255)
                                                            gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag         >  1:840988/1‑68 (MQ=255)
                                                            gAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTcagcag         >  1:2969346/1‑68 (MQ=255)
                                                                 ttCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGaa    >  1:2592703/1‑68 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:1168803/69‑1 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:2691994/69‑1 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:3255515/69‑1 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:2962544/69‑1 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:2832539/69‑1 (MQ=255)
                                                                  tCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATa  <  1:2700986/69‑1 (MQ=255)
                                                                    |                                                                  
CGGATCAGGCAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATA  >  minE/224001‑224135

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: