Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,502,076 A→G N347D (AAC→GAC)  rlhA → 23S rRNA 5‑hydroxycytidine C2501 synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,502,0760AG94.1% 45.2 / ‑4.0 17N347D (AAC→GAC) rlhA23S rRNA 5‑hydroxycytidine C2501 synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (5/11);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGG  >  NZ_CP009273/1501989‑1502156
                                                                                       |                                                                                
ggCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGAc                                                                                <  1:388096/90‑1 (MQ=255)
    tGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACgtgt                                                                            <  2:123748/90‑1 (MQ=255)
         gTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTGGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGAtt                                                                       <  2:826617/90‑1 (MQ=255)
                     gtGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAgtc                                                           >  1:782361/1‑90 (MQ=255)
                                        atcTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATa                                        >  1:790632/1‑90 (MQ=255)
                                         tctcGATGTTGCCGTTCCCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATAc                                       <  2:233503/90‑1 (MQ=255)
                                         tctcGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATAc                                       <  1:266145/90‑1 (MQ=255)
                                             gATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTc                                   >  1:284570/1‑90 (MQ=255)
                                                 ttGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCgaga                               >  2:166696/1‑90 (MQ=255)
                                                           cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  2:652587/90‑1 (MQ=255)
                                                           cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  1:416714/90‑1 (MQ=255)
                                                           cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  2:19309/90‑1 (MQ=255)
                                                                      cAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTAcc          <  1:505927/90‑1 (MQ=255)
                                                                       aaaTGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACcg         <  1:183716/90‑1 (MQ=255)
                                                                            gCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCt    <  1:166696/90‑1 (MQ=255)
                                                                              gATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTgg  >  1:10127/1‑90 (MQ=255)
                                                                              gATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTgg  >  1:431051/1‑90 (MQ=255)
                                                                                       |                                                                                
GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGG  >  NZ_CP009273/1501989‑1502156

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