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,4510TC95.5% 58.9 / ‑3.4 22E162G (GAA→GGA) mltBlytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (5/16);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

GCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAG  >  NZ_CP009273/2818375‑2818537
                                                                            |                                                                                      
gcgTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAAt                                                                           >  2:215602/1‑90 (MQ=255)
               tGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCTAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTc                                                            >  1:469710/1‑90 (MQ=255)
                gACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTcc                                                           <  1:671889/90‑1 (MQ=255)
                   aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                                        >  2:687263/1‑90 (MQ=255)
                     cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCcggcgg                                                      <  1:50458/90‑1 (MQ=255)
                        tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCCCGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                                   <  2:33074/90‑1 (MQ=255)
                        tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                                   <  2:130629/90‑1 (MQ=255)
                        tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                                   <  2:18013/90‑1 (MQ=255)
                        tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                                   <  2:155481/90‑1 (MQ=255)
                        tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                                   <  1:421501/90‑1 (MQ=255)
                                            gAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAcgcgc                               <  2:472977/90‑1 (MQ=255)
                                                      cATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGc                     <  1:591188/90‑1 (MQ=255)
                                                      cATCACCGGCCCCCAGCGGGTTCCAACGCCGATAATCCCCCCGATAATTTCCGGCGGTACTCCAAACACCTGCCACGCGCGATTCAACGc                     <  1:687263/90‑1 (MQ=255)
                                                        tCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAt                   <  1:531340/90‑1 (MQ=255)
                                                        tCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAt                   >  1:424142/1‑90 (MQ=255)
                                                        tCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAt                   <  2:424142/90‑1 (MQ=255)
                                                        tCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAt                   <  2:491959/90‑1 (MQ=255)
                                                                    aGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGAt       <  1:198492/90‑1 (MQ=255)
                                                                         gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAAc                       <  1:387050/69‑1 (MQ=255)
                                                                         gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAg  <  1:591388/90‑1 (MQ=255)
                                                                         gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAg  <  2:249899/90‑1 (MQ=255)
                                                                         gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAg  >  2:591388/1‑90 (MQ=255)
                                                                            |                                                                                      
GCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAG  >  NZ_CP009273/2818375‑2818537

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