Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 92,505 C→G A6G (GCC→GGC)  mraY → phospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927392,5050CG85.0% 117.0 / 14.2 40A6G (GCC→GGC) mraYphospho‑N‑acetylmuramoyl‑pentapeptide‑ transferase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base G (20/14);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

CGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCT  >  NZ_CP009273/92419‑92589
                                                                                      |                                                                                    
cGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGa                                                                                   <  2:569770/90‑1 (MQ=255)
cGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGa                                                                                   <  1:196398/90‑1 (MQ=255)
     ttAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAtt                                                                              <  1:901344/90‑1 (MQ=255)
      tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt                                                                             <  1:819652/90‑1 (MQ=255)
      tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt                                                                             <  1:302702/90‑1 (MQ=255)
      tAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACAttt                                                                             <  2:390733/90‑1 (MQ=255)
           aaGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCa                                                                        <  2:117380/90‑1 (MQ=255)
                      tAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTcc                                                             <  2:705500/90‑1 (MQ=255)
                      tAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTcc                                                             <  1:948899/90‑1 (MQ=255)
                       aGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCg                                                            >  1:264720/1‑90 (MQ=255)
                       aGTGCCGCCATGGAAGAGGCAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCg                                                            >  2:453718/1‑90 (MQ=255)
                        gTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTCAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCgg                                                           >  2:811342/1‑90 (MQ=255)
                         tGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGc                                                          >  1:676404/1‑90 (MQ=255)
                           ccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCtt                                                        <  2:664132/90‑1 (MQ=255)
                              ccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTaa                                                     <  2:382628/90‑1 (MQ=255)
                              ccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTaa                                                     <  1:707442/90‑1 (MQ=255)
                                          tagtaCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTa                                         <  1:941512/90‑1 (MQ=255)
                                              aCGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTg                                     >  2:987590/1‑90 (MQ=255)
                                                gcgcTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGAc                                   >  2:574695/1‑90 (MQ=255)
                                                     ttACAGGAGAATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTc                              >  2:617745/1‑90 (MQ=255)
                                                             gaATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATc                      >  2:48231/1‑90 (MQ=255)
                                                             gaATGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATc                      >  2:452345/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  2:631976/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  1:437551/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  1:369574/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  2:91883/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCCGACGTTTCGCGCCATCGt                    >  1:670094/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  2:760404/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  1:95898/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  1:665274/1‑90 (MQ=255)
                                                               aTGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGt                    >  1:646298/1‑90 (MQ=255)
                                                                 gggACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCa                  <  2:439090/90‑1 (MQ=255)
                                                                  ggACATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAg                 >  1:556662/1‑90 (MQ=255)
                                                                     cATGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCct              >  1:225593/1‑90 (MQ=255)
                                                                       tGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc            >  1:991218/1‑90 (MQ=255)
                                                                       tGTTAGTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCctgc            >  2:704137/1‑90 (MQ=255)
                                                                            gTTTGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGAcc       >  2:450639/1‑90 (MQ=255)
                                                                              ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc     <  2:556662/90‑1 (MQ=255)
                                                                              ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc     <  1:574695/90‑1 (MQ=255)
                                                                              ttGGCTGGGCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgc     >  1:424158/1‑90 (MQ=255)
                                                                               tggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACcgcg    >  1:343223/1‑90 (MQ=255)
                                                                                ggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCgcgc   <  1:720373/90‑1 (MQ=255)
                                                                                ggctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCgcgc   <  1:631976/90‑1 (MQ=255)
                                                                                 gctgggCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCt  >  1:760222/1‑90 (MQ=255)
                                                                                      |                                                                                    
CGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCT  >  NZ_CP009273/92419‑92589

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