Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,818,451 T→C E162G (GAA→GGA)  mltB ← lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,818,4510TC80.0% 14.1 / 0.4 10E162G (GAA→GGA) mltBlytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (1/7);  total (2/8)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGA  >  NZ_CP009273/2818380‑2818531
                                                                       |                                                                                
ggTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCcga                                                                <  2:160832/90‑1 (MQ=255)
           gACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTcc                                                     <  2:124718/90‑1 (MQ=255)
                cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg                                                <  1:325083/90‑1 (MQ=255)
                  ttGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTa                                              <  1:14722/90‑1 (MQ=255)
                   tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                             <  1:12485/90‑1 (MQ=255)
                                       gAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAcgcgc                         <  2:347696/90‑1 (MQ=255)
                                              ccccATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                             <  1:371366/63‑1 (MQ=255)
                                              ccccATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                             >  2:371366/1‑63 (MQ=255)
                                                 cATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGc               <  1:169239/90‑1 (MQ=255)
                                                              cAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGa  >  2:222670/1‑90 (MQ=255)
                                                                       |                                                                                
GGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGA  >  NZ_CP009273/2818380‑2818531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑