Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,494,123	T→C	E195G (GAA→GGA)	ydcK ←	YdcK family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,494,123	0	T	C	93.9%	92.2 / ‑0.4	33	E195G (GAA→GGA)	ydcK	YdcK family protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (22/9);  total (24/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

CATCTTCTTCTATGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTTCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTTGTGCCTGG  >  NZ_CP009273/1494037‑1494206                                                                                       |                                                                                    cATCTTCTTCTATGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCaa                                                                                  >  2:682541/1‑90 (MQ=255)     ttcttctATGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACt                                                                              >  1:627927/1‑90 (MQ=255)         tctATGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCa                                                                          >  2:109247/1‑90 (MQ=255)            aTGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCaaaa                                                                       >  2:108816/1‑90 (MQ=255)            aTGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCaaaa                                                                       >  1:725737/1‑90 (MQ=255)                 aGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCaaa                                                                  >  2:173899/1‑90 (MQ=255)                            gcgcGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGgcgc                                                       <  2:1256629/90‑1 (MQ=255)                            gcgcGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGgcgc                                                       >  2:353538/1‑90 (MQ=255)                                 tGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGt                                                  >  1:760569/1‑90 (MQ=255)                                      gTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAAt                                             >  1:54189/1‑90 (MQ=255)                                      gTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAAt                                             >  1:999350/1‑90 (MQ=255)                                      gTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAAt                                             <  1:567350/90‑1 (MQ=255)                                             ttAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGa                                                     >  2:1547187/1‑75 (MQ=255)                                             ttAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGa                                                     <  1:1547187/75‑1 (MQ=255)                                             ttAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACATCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCa                                      >  1:1205553/1‑90 (MQ=255)                                                     aTCGCATAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGgac                              >  2:1510606/1‑90 (MQ=255)                                                          aTAGCCAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACgg                         >  1:371927/1‑90 (MQ=255)                                                              ccAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATc                     <  1:1178150/90‑1 (MQ=255)                                                              ccAGACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATc                     <  2:1319225/90‑1 (MQ=255)                                                                 gACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGcc                  >  1:1516113/1‑90 (MQ=255)                                                                 gACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGcc                  >  1:1247171/1‑90 (MQ=255)                                                                 gACATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGcc                  >  2:446889/1‑90 (MQ=255)                                                                   cATTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCAt                <  2:996469/90‑1 (MQ=255)                                                                    aTTGTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATa               <  2:1247171/90‑1 (MQ=255)                                                                       gTTCTCTTCATTACCTTCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAAt            >  2:934351/1‑90 (MQ=255)                                                                       gTTCTCTTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAAt            >  2:346331/1‑90 (MQ=255)                                                                            ctTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTtgtg       >  1:456765/1‑90 (MQ=255)                                                                            ctTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTtgtg       <  2:300106/90‑1 (MQ=255)                                                                            ctTCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTtgtg       <  1:108816/90‑1 (MQ=255)                                                                              tCATTACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTTGTGcc     >  1:1025862/1‑90 (MQ=255)                                                                               cATTACCTTCAATACTGGGAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTTGTGCCt    >  1:710713/1‑90 (MQ=255)                                                                                 ttACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTTGTGCCtgg  >  1:1137616/1‑90 (MQ=255)                                                                                 ttACCTCCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCACACCGGACGACGGCATCGCCATAAATTTGTGCCtgg  >  1:1249767/1‑90 (MQ=255)                                                                                       |                                                                                    CATCTTCTTCTATGCCAGCTTTCACCTGCGCGTGACCGTAGACTTTAGCGCAATCGCATAGCCAGACATTGTTCTCTTCATTACCTTCAATACTGGCAAAATCAAAAACTTCGGCGCGATGTTCAATAAAAGCATACCGGACGACGGCATCGCCATAAATTTGTGCCTGG  >  NZ_CP009273/1494037‑1494206

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