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_CP009273	1,502,076	0	A	G	96.4%	80.3 / ‑3.4	28	N347D (AAC→GAC)	rlhA	23S rRNA 5‑hydroxycytidine C2501 synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (9/18);  total (9/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGG  >  NZ_CP009273/1501993‑1502156                                                                                    |                                                                                 tGCCGGTTGGCGAAGTAGTGAAAGTGACAAAAGATCACCGCCTTGTTGCGGTTCCCGAGCCACTGGCAAATGGCGATGGCCTGGACgtgt                                                                            <  1:1342985/90‑1 (MQ=255)       tAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTa                                                                      <  2:707936/90‑1 (MQ=255)            gAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGt                                                                 >  2:360135/1‑90 (MQ=255)              agtagtGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGa                                                               <  1:272985/90‑1 (MQ=255)                     aaaGTGGCGAAAGATCATCTCGATGTTGCCGTTCCCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAgtcgtc                                                        <  2:315871/90‑1 (MQ=255)                         tGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGtt                                                    >  1:561415/1‑90 (MQ=255)                            cGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTc                                                 <  2:388478/90‑1 (MQ=255)                            cGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTc                                                 <  1:266486/90‑1 (MQ=255)                            cGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTc                                                 <  2:41943/90‑1 (MQ=255)                                      tctcGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATAc                                       <  1:185607/90‑1 (MQ=255)                                      tctcGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATAc                                       <  1:678618/90‑1 (MQ=255)                                      tctcGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATAc                                       <  2:1046812/90‑1 (MQ=255)                                           aTGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCg                                  >  1:859434/1‑90 (MQ=255)                                               tGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGaa                              >  1:1508552/1‑90 (MQ=255)                                                gCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGaaa                             >  2:773446/1‑90 (MQ=255)                                                   gTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAAcc                          <  1:567975/90‑1 (MQ=255)                                                     tACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCgg                        >  2:1469828/1‑90 (MQ=255)                                                     tACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCgg                        >  2:1571544/1‑90 (MQ=255)                                                     tACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCgg                        >  2:899305/1‑90 (MQ=255)                                                        cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  2:405815/90‑1 (MQ=255)                                                        cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  1:899305/90‑1 (MQ=255)                                                        cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  1:323159/90‑1 (MQ=255)                                                        cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  1:298085/90‑1 (MQ=255)                                                        cGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga                     <  1:141350/90‑1 (MQ=255)                                                          aGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGaaa                   <  2:561415/90‑1 (MQ=255)                                                                    aaaTGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACcg         <  1:1533723/90‑1 (MQ=255)                                                                         gCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCt    <  1:770397/90‑1 (MQ=255)                                                                           gATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTgg  >  2:477815/1‑90 (MQ=255)                                                                                    |                                                                                 TGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGG  >  NZ_CP009273/1501993‑1502156

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