Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,630,540	A→C	31.0%	L419R (CTG→CGG)	czcS ←	Sensor protein CzcS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,630,540	0	A	C	31.0%	44.9 / 2.8	44	L419R (CTG→CGG)	czcS	Sensor protein CzcS
Reads supporting (aligned to +/- strand):  ref base A (16/13);  new base C (11/4);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGCC  >  ACB122/1630408‑1630671                                                                                                                                     |                                                                                                                                    gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCCGTCccgcccgcccgc                                                                                                                                   >  3:433354/1‑135 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGc                                                                                                                                   >  7:86587/1‑135 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGc                                                                                                                                   >  2:24481/1‑135 (MQ=255)   cGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcg                                                                                                                                 <  1:198375/135‑1 (MQ=255)    gATCACTCGGCAGGAAACACCAGGGTAATCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgc                                                                                                                                <  7:140305/135‑1 (MQ=255)    gATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgc                                                                                                                                <  2:94271/135‑1 (MQ=255)         cTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                           <  4:209187/135‑1 (MQ=255)         cTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                           <  6:252701/135‑1 (MQ=255)               aGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGa                                                                                                                     >  2:29687/1‑135 (MQ=255)                 gAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                   <  7:96975/135‑1 (MQ=255)                 gAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                   <  1:298329/135‑1 (MQ=255)                 gAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                   <  7:130469/135‑1 (MQ=255)                           ggTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGCCCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTACCGCCCACCCGCCGCGGCCCCCCCCCACCAACCCGGt                                                                                                          >  7:339054/1‑134 (MQ=255)                           ggTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGt                                                                                                          >  8:254156/1‑134 (MQ=255)                           ggTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGt                                                                                                          >  4:137959/1‑134 (MQ=255)                           ggTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCGCCCGCCGCGAACGCCCCGACCCTACCCGGt                                                                                                          >  5:85788/1‑134 (MQ=255)                                   gggTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGCTCGCCAACCCCCACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGCACCTACCCGGTAAAACCGCt                                                                                                 >  2:306793/1‑135 (MQ=255)                                   gggTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGAACCAACCCGGGAAAACCCCt                                                                                                 >  6:137864/1‑135 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCCCGCCCGCCAGCCCCGCACGCCCCCCACCAACCCGCTAAAACCCCTCGac                                                                                             >  6:91890/1‑133 (MQ=255)                                         gCCCCAGGGGGAAGCGGCTGACGTCGACCCGCCCGCCATGCAAGTGCATGATCGACTTCACGTGCCCCAACCCCAACCCGGTCCCGCCCTCCCCCCGCGAACGCCCCTAACCTACCCGTTAAAACCGCTCTAACa                                                                                           <  8:443665/135‑1 (MQ=255)                                                 gggCAGCTGCTGACCTCGCCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCAGGTCCCGCCCTCCAGCCGCGCCCGCCCCGAACCTACCCCGTAAAGCCGCTCGGACAAGCGCGAc                                                                                   >  7:169924/1‑135 (MQ=255)                                                 gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACCTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGCCCGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGcgcgcc                                                                                   >  1:263251/1‑133 (MQ=255)                                                      gctgctGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTACCGCCCACCCGCCGCGCACCCCCCGCACCGACCCCGTAAAACCGCTCCAAAAAGCCCGACAGGTg                                                                              >  1:136455/1‑135 (MQ=255)                                                      gctgctGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCCGTCCCGCCCTCCCGCCGCGCACGCCCCGACCCTACCCGGTTAAACCGCTCGCACAAGCGCGACGGGTg                                                                              >  4:432185/1‑135 (MQ=255)                                                       ctgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCCGTCCCCCCCTCACGCCGCCGCCCCCCCCCACCCACCCCGTAAAAACCCTCGAACCAGCCGGACGGGTGc                                                                             >  6:273814/1‑135 (MQ=255)                                                        tgctgcCCACGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCACCCCCACCCCGGGCCCGCCCTCCAGCCGCGCACGCCCCGAACCGACCCCGTAAAACCGCTCGAACACGCGGGCCTGGTGct                                                                            >  2:261890/1‑135 (MQ=255)                                                        tgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCCCGCCCTCCAGCCGCCACCGCCCCGCACCTCCCCGGGAAAACCCCTCCACCAAGCGCCACAGGTGcg                                                                            >  1:130607/1‑134 (MQ=255)                                                        tgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACCCCCCGAACCTACCCCGGAAAACCCCTCGAACAACCGCGCCAGGTGct                                                                            >  8:386514/1‑135 (MQ=255)                                                                         ccgccATGCAACTGCATGATCTACTTCACGATCGCCAACCCCACCCCGGACCCGCCCTCCCGCCGCGACCGCCCCGAACATACCCGTTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGa                                                           <  1:29687/135‑1 (MQ=255)                                                                             cATGCAACTGCATGATCGACTTCACGATCGCCACCCCCAACCCGGGCCCGCCCTCCAGCCGCGAACGCCCCGACCCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCt                                                       <  7:386514/135‑1 (MQ=255)                                                                               tGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGGACCGCCCGCCCGCCGCGGACCCCCCCCACCGACCCGGTAAAACCGCGCGACCAATCGCGACAGGTGCTCTGCAGGCACAACCGGAcccgtg                                                     >  6:55080/1‑132 (MQ=255)                                                                                   acgggaTGATCGACTTCACGAGCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCTAACGCCCCGAACCTACCCGGTAAAACCGCGCGAACAACCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                                 <  6:308567/130‑1 (MQ=255)                                                                                   aCTGCATGATCGACTTCACGATCGCCACCCCCAACCCGGGCCCGCCCGCCACCCGCGAACGCCCCGAACCTCCCCGGTAAAACCGCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                                 <  7:41309/135‑1 (MQ=255)                                                                                   aCTGCATGATCGACTTCACGAGCGCCAACCCCAACCCGGCCCCGCCCTCCAGCCGCGACCGCCCCGAACCTACCCGGAAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                                 <  6:102733/135‑1 (MQ=255)                                                                                        aTGATCGACTGCACGATCGCCAACCCCAACCGGGGCCGGCCCTCCCGCCGCGAACGCCCGGAAACTACCCGTAAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCAc                                             <  4:440084/134‑1 (MQ=255)                                                                                          gATCGACTTCACGATCGCCAACCCCAACCCGGGCCCGCCCTCCAGCCGCGAACGCCCCGCACCTACCCGGGAAAACCCCTCGAACAATCGCGCCTGGTGCTATGCCGGCACACCCCGCCCCTGGGTCTCCACgg                                           >  5:108138/1‑134 (MQ=255)                                                                                            tCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCACCAGCCGCCAACGCCCCGCACCTACCAGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCAGCAGGCACACCCCGACCCAGGTTATCCAccg                                           >  8:271018/1‑130 (MQ=255)                                                                                                           ccacccccaaccCGGACCCGCCCTCCAGCCGCGAACGCCCCTAACCTCCCCGGTAAAACCGCTCTAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGc                         <  5:55080/135‑1 (MQ=255)                                                                                                                   aacccGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGTAAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGt                 <  2:91065/135‑1 (MQ=255)                                                                                                                     cccGGTCCCGCCCTCCAGCCGCGCACGCCCCGAACCTACCCCGTAAAACCGCTCGCACAAGCCCGCCTGGTGCGCTGCAGGCACACCCGGACCCAGGTTGTCCACGGCTAGAATGACATAGCTGCCCGCGCGGTc               >  6:108572/1‑135 (MQ=255)                                                                                                                             cGCCCTCCCGCCGCGCACGCCCCGAACCAACCCGGTAAAACCGCTCGAACAAGCGCGACAGGTGCTCTGCAGGCCCCCCCGGCCCCAGGTTCTCCACGGATAGAAAGACTTGGCTGCCCGAGCGTTCGCTTCGCa       >  4:81335/1‑135 (MQ=255)                                                                                                                             cGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgc          >  7:43339/1‑132 (MQ=255)                                                                                                                             cGCCCTCCAGCCGCCCACGCCCCGAACCAACCCGCTAAAACCGCTCGAACAAGCGCGACTGGTGGTCTGCAGGCACACCCGGGCCCTGGGTCTACCAGGATAGAAATTCGTAGCTTCCCGTGAGTTCGATgggcc       >  8:444875/1‑134 (MQ=255)                                                                                                                              gCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgcgc        >  8:53244/1‑133 (MQ=255)                                                                                                                                  tCCAGCCGCGAACCCCCCGAACCTACCCGGGAAAACCGCTCGACCAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATACAAAGCCGTCGCTGCCCGAGCGGTCGAGGCGCAGAGcc  >  7:428047/1‑135 (MQ=255)                                                                                                                                     |                                                                                                                                    GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGCC  >  ACB122/1630408‑1630671

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