Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,302,292 A→G 43.5% S263S (AGT→AGC dacD ← D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,302,2920AG43.5% 6.4 / 33.7 46S263S (AGT→AGCdacDD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base A (6/20);  new base G (5/15);  total (11/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.84e-01

GTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATGAGATTA  >  minE/1302238‑1302358
                                                      |                                                                  
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:2298621/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:3358929/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:1135112/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:1210140/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:1992600/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:2625337/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:2518744/69‑1 (MQ=255)
gTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:743353/69‑1 (MQ=255)
 ttttGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAtt                                                     <  1:1171543/69‑1 (MQ=255)
 ttttGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAtt                                                     <  1:81453/69‑1 (MQ=255)
 ttttGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCAt                                                      <  1:1837964/68‑1 (MQ=255)
   ttgttgCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTAc                                                   <  1:2700608/69‑1 (MQ=255)
            ccAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCGCTGTCAGCCCCCATTACCACTGc                                             >  1:937130/1‑66 (MQ=255)
            ccAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCGCTGTCAGCCCCCATTACCACTGCa                                            >  1:2331845/1‑67 (MQ=255)
            ccAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCGCTGTCAGCCCCCATTACCACTGCa                                            >  1:704630/1‑67 (MQ=255)
            ccAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCGCTGTCAGCCCCCATTACCACTGCa                                            >  1:1369474/1‑67 (MQ=255)
            ccAACGCAGTAATTTTCTTGCATCTTCCTCACGACCTTTTGCGCTGTCAGCCCCCATTACCACTGCa                                            >  1:1462830/1‑67 (MQ=255)
                    gTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATGAgacgac                                  <  1:2859586/69‑1 (MQ=255)
                    gTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATGAgacgac                                  <  1:3208477/69‑1 (MQ=255)
                    gTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATGAgacgac                                  <  1:885145/69‑1 (MQ=255)
                       aTTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTGAGCCCCCATTACCACTGc                                             >  1:1122633/1‑55 (MQ=255)
                       aTTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGc                                             >  1:1731757/1‑55 (MQ=255)
                       aTTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGc                                             >  1:1621532/1‑55 (MQ=255)
                       aTTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGc                                             >  1:1536639/1‑55 (MQ=255)
                          ttCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcccc                                                         <  1:635208/40‑1 (MQ=255)
                          ttCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcccc                                                         <  1:2660202/40‑1 (MQ=255)
                          ttCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcccc                                                         <  1:2750049/40‑1 (MQ=255)
                          ttCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGcccc                                                         <  1:2750918/40‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCGCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:806935/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:96011/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:961975/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:261315/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:3009444/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:2859397/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:1078286/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:2519954/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:2419094/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:2148964/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:1846603/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:150838/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:1359639/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATga       <  1:1112068/65‑1 (MQ=255)
                                                   tGCGCTGTCAGCCCCCATTACCACTGCAATGAGACGAAGCTGCCCATCTACAGCCGAAGCAATga       <  1:409604/65‑1 (MQ=255)
                                                    gCACTGTCAGCCCCCATTACCACTGCAATGAGACGCCGCTGCCCATCTACAGCCGAAGCAATGAGAtt   >  1:2538360/1‑68 (MQ=255)
                                                    gCACTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATGAGAt    <  1:3272255/67‑1 (MQ=255)
                                                    gCACTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATGAGATTa  >  1:485721/1‑69 (MQ=255)
                                                      |                                                                  
GTTTTGTTGCCCCCAACGCAGTAATTTTCTTGCCTCTTCCTCACGACCTTTTGCACTGTCAGCCCCCATTACCACTGCAATGAGACGACGCTGCCCATCTACAGCCGAAGCAATGAGATTA  >  minE/1302238‑1302358

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: