Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,839,507 T→G 38.7% D367A (GAC→GCC)  ampC ← beta‑lactamase/D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,839,5070TG38.7% 15.0 / 37.3 44D367A (GAC→GCC) ampCbeta‑lactamase/D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base G (6/11);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATA  >  minE/2839449‑2839575
                                                          |                                                                    
ctaaaaTTCGGACCCGATGGAATTTTTCTGTAGAGCGTTAAGAATCTGCCAGGCGGCGtg                                                                     <  1:2031654/58‑2 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:3190809/1‑58 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:358097/1‑58 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:1412847/1‑58 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:1616377/1‑58 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:626521/1‑58 (MQ=255)
  aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     >  1:2072472/1‑58 (MQ=255)
   aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                                                         <  1:3360471/69‑1 (MQ=255)
   aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                                                         <  1:196801/69‑1 (MQ=255)
    aaTTCGTACCCGATGGATTTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                                                         <  1:2424583/68‑1 (MQ=255)
    aaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                                                     <  1:468141/56‑1 (MQ=255)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:2642217/1‑48 (MQ=38)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:406329/1‑48 (MQ=38)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:342263/1‑48 (MQ=38)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:2819132/1‑48 (MQ=38)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:2941116/1‑48 (MQ=38)
                    aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                                             >  1:1260785/1‑48 (MQ=38)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:3175508/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:2705753/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:3206381/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:3424435/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:2473209/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:450178/69‑1 (MQ=255)
                       tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:1681018/69‑1 (MQ=255)
                       tttACAGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                     <  1:1311422/69‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTTGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:2136937/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:765093/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:840717/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:871082/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:2905296/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:2814618/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:1364654/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:1570881/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:1814121/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:2195627/66‑1 (MQ=255)
                                gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa                               <  1:1975702/66‑1 (MQ=255)
                                           aTCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttt                                            >  1:1225736/1‑42 (MQ=255)
                                           aTCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGc                 >  1:605736/1‑69 (MQ=255)
                                           aTCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGc                 >  1:393036/1‑69 (MQ=255)
                                           aTCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGc                 >  1:77637/1‑69 (MQ=255)
                                           aTCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGc                 >  1:827891/1‑69 (MQ=255)
                                            tCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGct                >  1:2123671/1‑69 (MQ=255)
                                            tCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGct                >  1:656656/1‑69 (MQ=255)
                                            tCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGct                >  1:750675/1‑69 (MQ=255)
                                            tCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGct                >  1:330767/1‑69 (MQ=255)
                                                          tCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATa  <  1:741345/69‑1 (MQ=255)
                                                          tCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATa  <  1:2307558/69‑1 (MQ=255)
                                                          tCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATa  <  1:2616897/69‑1 (MQ=255)
                                                          |                                                                    
CGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATA  >  minE/2839449‑2839575

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: