Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,439,101 T→C K255E (AAG→GAG)  mepA ← penicillin‑insensitive murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,439,1010TC89.5% 41.6 / 0.1 19K255E (AAG→GAG) mepApenicillin‑insensitive murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (10/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATG  >  NZ_CP009273/2439016‑2439175
                                                                                     |                                                                          
aaaCAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                        <  2:608563/90‑1 (MQ=255)
aaaCAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                        <  1:456579/90‑1 (MQ=255)
     ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                   >  2:124920/1‑90 (MQ=255)
         aTTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTttgtt                                                               <  1:221887/90‑1 (MQ=255)
               cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTccc                                                         <  2:54114/90‑1 (MQ=255)
               cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                         >  2:552258/1‑90 (MQ=255)
               cGTTTCCAATAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                         <  1:816029/90‑1 (MQ=255)
                 tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                       >  1:144738/1‑90 (MQ=255)
                   tCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGtt                                                     >  1:177594/1‑90 (MQ=255)
                          gATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAgg                                              >  1:418106/1‑90 (MQ=255)
                              aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                          >  1:663746/1‑90 (MQ=255)
                              aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                          >  1:159843/1‑90 (MQ=255)
                                     cATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGc                                   >  1:25107/1‑90 (MQ=255)
                                       tCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCtt                                 >  2:480428/1‑90 (MQ=255)
                                          aGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGc                              <  2:767542/90‑1 (MQ=255)
                                            tAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAg                            >  2:265955/1‑90 (MQ=255)
                                             aGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGt                           >  1:513883/1‑90 (MQ=255)
                                                                      gAGGCTGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATg  <  2:508295/90‑1 (MQ=255)
                                                                      gAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATg  <  2:429755/90‑1 (MQ=255)
                                                                                     |                                                                          
AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATG  >  NZ_CP009273/2439016‑2439175

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