Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,174,381	T→G	C270G (TGC→GGC)	yegU →	ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,174,381	0	T	G	84.6%	68.2 / 5.4	26	C270G (TGC→GGC)	yegU	ADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (13/9);  total (15/11)
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.16e-01

GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG  >  NZ_CP009273/2174297‑2174462                                                                                     |                                                                                  gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                              >  2:614479/1‑90 (MQ=255) gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                              >  1:767250/1‑90 (MQ=255)              gTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTa                                                                 <  2:36310/90‑1 (MQ=255)                        tATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                        <  1:812587/73‑1 (MQ=255)                        tATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                        >  2:812587/1‑73 (MQ=255)                         aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                      >  1:223109/1‑90 (MQ=255)                         aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                      >  2:377245/1‑90 (MQ=255)                              gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                        <  1:462182/67‑1 (MQ=255)                              gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                        >  2:462182/1‑67 (MQ=255)                              gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                 >  2:544077/1‑90 (MQ=255)                                 cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                              >  1:154254/1‑90 (MQ=255)                                 cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                              >  2:62949/1‑90 (MQ=255)                                 cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                              >  1:700445/1‑90 (MQ=255)                                   ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                            >  2:280765/1‑90 (MQ=255)                                   ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                            <  1:372293/90‑1 (MQ=255)                                       gtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGTTAACCTTGGCGGCGACACAGACACCAt                                        >  1:408340/1‑90 (MQ=255)                                                        ttGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGaca                                                  >  1:122421/1‑63 (MQ=255)                                                        ttGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGaca                                                  <  2:122421/63‑1 (MQ=255)                                                                gCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc               >  1:537334/1‑90 (MQ=255)                                                                gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       >  1:783340/1‑66 (MQ=255)                                                                gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                       <  2:783340/66‑1 (MQ=255)                                                                  acaAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcgg             <  1:618445/90‑1 (MQ=255)                                                                            ccGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt   <  2:704463/90‑1 (MQ=255)                                                                             cGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg  <  1:645499/90‑1 (MQ=255)                                                                             cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg  <  1:377245/90‑1 (MQ=255)                                                                             cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg  <  2:46205/90‑1 (MQ=255)                                                                                     |                                                                                  GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATG  >  NZ_CP009273/2174297‑2174462

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