Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,593,447	T→C	S52P (TCC→CCC)	yjjQ →	response regulator transcription factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,593,447	0	T	C	88.5%	60.8 / 2.9	26	S52P (TCC→CCC)	yjjQ	response regulator transcription factor
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (16/7);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAA  >  NZ_CP009273/4593364‑4593530                                                                                    |                                                                                    agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                               >  1:450612/1‑90 (MQ=255)  ggTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAtt                                                                              >  1:317977/1‑90 (MQ=255)   gTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTa                                                                             >  2:398051/1‑90 (MQ=255)                        gAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGt                                                        <  2:450090/90‑1 (MQ=255)                         aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTg                                                       >  2:66306/1‑90 (MQ=255)                         aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGACTAGTCATTGCTGATTTAGCCGGTg                                                       >  1:295866/1‑90 (MQ=255)                           tatGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGaa                                                     <  2:226605/90‑1 (MQ=255)                                           gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                     >  2:320824/1‑90 (MQ=255)                                           gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                     >  1:230802/1‑90 (MQ=255)                                            cgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTc                                    <  1:169543/90‑1 (MQ=255)                                            cgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTc                                    <  1:141193/90‑1 (MQ=255)                                             gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   <  1:289140/90‑1 (MQ=255)                                             gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   >  2:297860/1‑90 (MQ=255)                                                tctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGtt                                >  2:359104/1‑90 (MQ=255)                                                        ggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACa                        >  2:107388/1‑90 (MQ=255)                                                           ccttaccttaTTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   >  2:224056/1‑76 (MQ=255)                                                           ccttaccttaTTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   <  1:224056/76‑1 (MQ=255)                                                                  ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt              >  1:295496/1‑90 (MQ=255)                                                                  ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt              >  2:514749/1‑90 (MQ=255)                                                                   tattaCAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCtt             >  1:233673/1‑90 (MQ=255)                                                                    attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCttt            >  1:49085/1‑90 (MQ=255)                                                                         aaTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtct       >  1:200724/1‑90 (MQ=255)                                                                          aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc      <  1:276918/90‑1 (MQ=255)                                                                              cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  >  2:137427/1‑90 (MQ=255)                                                                              cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  >  2:204917/1‑90 (MQ=255)                                                                              cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  <  1:221464/90‑1 (MQ=255)                                                                                    |                                                                                    AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAA  >  NZ_CP009273/4593364‑4593530

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