Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,977 T→A 15.2% intergenic (‑33/+30) ampC ← / ← frdD penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9770TA15.2% 68.4 / 4.4 26intergenic (‑33/+30)ampC/frdDpenicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (15/7);  new base A (2/2);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG  >  NC_000913/4378836‑4379123
                                                                                                                                             |                                                                                                                                                   
ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                                                                                                                                              >  1:543724/1‑148 (MQ=255)
 tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGc                                                                                                                                            >  3:824092/1‑149 (MQ=255)
  ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGcc                                                                                                                                           <  4:824092/149‑1 (MQ=255)
                     tattcacaATATCGTTGATTTGTTGATGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAAGAAGGCGCAGAGCGTCGTTTTGAACATAGTGTCTGGTTTCCATACAAAACGGCCCGCCATAGTCGGGCCCGATTTACATTGGCGATg                                                                                                                         <  1:492381/144‑1 (MQ=255)
                     gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGc                                                                                                                        >  3:511278/1‑149 (MQ=255)
                                  gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACTGCCCGCCTATGGCGGGcc                                                                                                                                           <  3:904751/117‑1 (MQ=255)
                                  gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACTGCCCGCCTATGGCGGGcc                                                                                                                                           >  4:904751/1‑117 (MQ=255)
                                  gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                     >  4:886924/1‑139 (MQ=255)
                                  gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt                                                                                                                     <  3:886924/139‑1 (MQ=255)
                                   ttGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAAc                                                                                                           >  2:259030/1‑148 (MQ=255)
                                                gggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                             >  4:425725/1‑149 (MQ=255)
                                                gggCAGCAAATGTAGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGc                                                                                              >  1:685916/1‑148 (MQ=255)
                                                 ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                             <  2:543724/148‑1 (MQ=255)
                                                            tgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCTTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                 <  4:511278/149‑1 (MQ=255)
                                                                      ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagcag                                                                        <  1:259030/148‑1 (MQ=255)
                                                                                      gcgcAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAc                                                        <  2:685916/148‑1 (MQ=255)
                                                                                                                   tGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAg                          >  4:3204/1‑149 (MQ=255)
                                                                                                                     gTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGAt                        >  3:610387/1‑149 (MQ=255)
                                                                                                                                    ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt         >  3:729644/1‑149 (MQ=255)
                                                                                                                                     gCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGTCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt         >  1:589530/1‑148 (MQ=255)
                                                                                                                                       ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                 <  3:71003/74‑1 (MQ=255)
                                                                                                                                       ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa                                                                                 >  4:71003/1‑74 (MQ=255)
                                                                                                                                       ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa      >  4:49964/1‑149 (MQ=255)
                                                                                                                                       ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa      >  4:811996/1‑149 (MQ=255)
                                                                                                                                       ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa      >  3:350544/1‑149 (MQ=255)
                                                                                                                                           cATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTAACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACg  >  3:750089/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG  >  NC_000913/4378836‑4379123

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: