Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,102,646	T→C	G140G (GGA→GGG)	wzx ←	O16 family O‑antigen flippase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,102,646	0	T	C	91.7%	60.8 / ‑0.7	24	G140G (GGA→GGG)	wzx	O16 family O‑antigen flippase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (11/11);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC  >  NZ_CP009273/2102557‑2102732                                                                                          |                                                                                       gACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTttttctttt                                                                                        >  1:780984/1‑90 (MQ=255)       gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                                  >  1:13688/1‑90 (MQ=255)       gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                                  >  2:290329/1‑90 (MQ=255)         gtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGaaa                                                                                <  2:382436/90‑1 (MQ=255)           tCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACa                                                                              >  1:601141/1‑90 (MQ=255)             acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:452012/90‑1 (MQ=255)             acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:601141/90‑1 (MQ=255)                  aataaatGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGAt                                                                       <  2:410360/90‑1 (MQ=255)                                          cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                               >  1:729033/1‑90 (MQ=255)                                           aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            >  1:391251/1‑60 (MQ=255)                                           aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:391251/60‑1 (MQ=255)                                                 gATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGa                                        >  1:411099/1‑90 (MQ=255)                                                        tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                                 <  2:729033/90‑1 (MQ=255)                                                          ctactTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGAc                               >  1:582412/1‑90 (MQ=255)                                                          ctactTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGAc                               <  1:11456/90‑1 (MQ=255)                                                             ctTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCat                            <  1:172202/90‑1 (MQ=255)                                                               tAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCatat                          <  1:547548/90‑1 (MQ=255)                                                               tAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCatat                          >  2:114662/1‑90 (MQ=255)                                                                      aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg                   >  1:80532/1‑90 (MQ=255)                                                                      aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg                   >  2:70390/1‑90 (MQ=255)                                                                       gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg                  >  1:120646/1‑90 (MQ=255)                                                                           tttcatttTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATg              <  1:153260/90‑1 (MQ=255)                                                                             tcatttTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt            <  2:745253/90‑1 (MQ=255)                                                                             tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt            <  2:512925/90‑1 (MQ=255)                                                                                ttttttCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGt         <  2:15678/90‑1 (MQ=255)                                                                                    ttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTtata     >  2:814075/1‑90 (MQ=255)                                                                                       tttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTc  <  1:140425/90‑1 (MQ=255)                                                                                          |                                                                                       GACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC  >  NZ_CP009273/2102557‑2102732

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