Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 98,192 T→C I314T (ATT→ACT)  murC → UDP‑N‑acetylmuramate‑‑L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927398,1920TC95.8% 68.7 / ‑3.5 24I314T (ATT→ACT) murCUDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (14/9);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CGCCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTATTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAAT  >  NZ_CP009273/98114‑98271
                                                                              |                                                                               
cgcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGgcgc                                                                      >  2:527198/1‑90 (MQ=255)
 gcCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCt                                                                     >  1:62246/1‑90 (MQ=255)
     ggTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                 >  2:276424/1‑90 (MQ=255)
     ggTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaa                                                                 >  2:531197/1‑90 (MQ=255)
      gtcgtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGaaa                                                                <  2:212649/90‑1 (MQ=255)
         gtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCt                                                             <  2:45710/90‑1 (MQ=255)
         gtcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCt                                                             >  2:330479/1‑90 (MQ=255)
          tcATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCtt                                                            >  2:386036/1‑90 (MQ=255)
                   gcTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGAc                                                   <  1:11306/90‑1 (MQ=255)
                   gcTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGAc                                                   >  1:580673/1‑90 (MQ=255)
                    cTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTCTGCGGGCGCTTGAAAGCTTCCAGGGGACt                                                  >  1:529481/1‑90 (MQ=255)
                         cgccgcAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCg                                             <  1:123651/90‑1 (MQ=255)
                                              tgcTACGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGa                        >  1:260549/1‑90 (MQ=255)
                                                  aCGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTc                    <  2:338692/90‑1 (MQ=255)
                                                   cGGAAGAGGGCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTcc                   <  1:402349/90‑1 (MQ=255)
                                                          gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCg                                             >  2:470775/1‑57 (MQ=255)
                                                          gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCg                                             <  1:470775/57‑1 (MQ=255)
                                                          gggCATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGa            >  2:567472/1‑90 (MQ=255)
                                                             cATTGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGcc         <  2:580673/90‑1 (MQ=255)
                                                               ttGACGACGAGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAg       >  2:250856/1‑90 (MQ=255)
                                                                   cgacgaGGCTATTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGaa   <  2:256573/90‑1 (MQ=255)
                                                                   cgacgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGaa   >  2:223149/1‑90 (MQ=255)
                                                                   cgacgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGaa   <  1:527198/90‑1 (MQ=255)
                                                                    gacgaGGCTACTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAAt  >  2:1981/1‑90 (MQ=255)
                                                                              |                                                                               
CGCCAGGTCGTCATAACGCGCTGAACGCCGCAGCTGCGGTTGCGGTTGCTACGGAAGAGGGCATTGACGACGAGGCTATTTTGCGGGCGCTTGAAAGCTTCCAGGGGACTGGTCGCCGTTTTGATTTCCTCGGTGAATTCCCGCTGGAGCCAGTGAAT  >  NZ_CP009273/98114‑98271

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