Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,839,488 A→G 41.4% L373L (CTT→CTC ampC ← beta‑lactamase/D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,839,4880AG41.4% 10.6 / 27.2 41L373L (CTT→CTCampCbeta‑lactamase/D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base A (11/13);  new base G (6/11);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCA  >  minE/2839421‑2839546
                                                                   |                                                          
gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                            >  1:2216741/1‑68 (MQ=255)
         gAAAAGCGGAGAAAAGGTCCGAATATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                   >  1:3427789/1‑68 (MQ=255)
         gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                  >  1:972076/1‑69 (MQ=255)
         gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                  >  1:2935723/1‑69 (MQ=255)
         gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                  >  1:2622779/1‑69 (MQ=255)
         gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                  >  1:3381693/1‑69 (MQ=255)
                            ctaaaaTTCGGACCCGATGGAATTTTTCTGTAGAGCGTTAAGAATCTGCCAGGCGGCGtg                                        <  1:2031654/58‑2 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:358097/1‑58 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:1616377/1‑58 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:2072472/1‑58 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:1412847/1‑58 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:626521/1‑58 (MQ=255)
                              aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        >  1:3190809/1‑58 (MQ=255)
                               aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                            <  1:196801/69‑1 (MQ=255)
                               aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                            <  1:3360471/69‑1 (MQ=255)
                                aaTTCGTACCCGATGGATTTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                            <  1:2424583/68‑1 (MQ=255)
                                aaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                        <  1:468141/56‑1 (MQ=255)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:406329/1‑48 (MQ=38)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:342263/1‑48 (MQ=38)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:2941116/1‑48 (MQ=38)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:2819132/1‑48 (MQ=38)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:1260785/1‑48 (MQ=38)
                                                aaTTTTACTGTAGAGCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCg                                >  1:2642217/1‑48 (MQ=38)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:450178/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:2705753/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:3175508/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:3206381/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:2473209/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:1681018/69‑1 (MQ=255)
                                                   tttACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:3424435/69‑1 (MQ=255)
                                                   tttACAGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg        <  1:1311422/69‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTTGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:2136937/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:1975702/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:1814121/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:2195627/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:2814618/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:1570881/66‑1 (MQ=255)
                                                            gagCGTTGAGAATCTGCCAGGCGGCGGCGACTCTCGCTGGATTGGGGTAGTTTTTGTTAGCCAGCa  <  1:1364654/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)
                                                                   |                                                          
GACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCA  >  minE/2839421‑2839546

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: