Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,630,540	A→C	42.8%	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	42.8%	6.4 / 4.9	27	L419R (CTG→CGG)	czcS	Sensor protein CzcS
Reads supporting (aligned to +/- strand):  ref base A (11/3);  new base C (6/6);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAG  >  ACB122/1630408‑1630669                                                                                                                                     |                                                                                                                                    gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGCCCTCGACCCGCCCGCCATGCACCTGCATGATCGCCTTCACGATCGCCAACCCCAACCCCGTCCCGCCCTCCAGc                                                                                                                                   >  5:104743/1‑135 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGCCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCccgcccgcccgc                                                                                                                                   >  6:214072/1‑135 (MQ=255)   cGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcg                                                                                                                                 <  1:163659/135‑1 (MQ=255)    gATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTAACGATCGCCCACCCCAACCCCCTCCCGCCCTCCAGcccc                                                                                                                                >  2:285459/1‑133 (MQ=255)      tCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGGCCCGCCCACCCGccccca                                                                                                                              >  1:129006/1‑131 (MQ=255)               aGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGa                                                                                                                     >  5:149017/1‑135 (MQ=255)                 gAAACACCAGGGTAAACCGGGTAGGCCCCAGGCGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                   <  3:272784/135‑1 (MQ=255)                         aGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCgg                                                                                                           >  7:57270/1‑135 (MQ=255)                              aaaCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGAACCTACCCGGG‑aaaa                                                                                                     >  1:43651/1‑135 (MQ=255)                                       aGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCCCGCCCTCCCGCCGCGCCCGCCCCGCCCCGACCCGGT‑AAAACC‑GCTCGaa                                                                                           >  4:43531/1‑135 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGCTCGACTTCACGATCGCCAACCCCCCCCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GGGCGaaa                                                                                          >  7:235924/1‑134 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGCCTTCACGATCGCCAACCCCCACCCGGTCCCGCCCTCCCGCCCCGAACGCCCCGAACCTACCCGGT‑AAGACC‑GCTCgga                                                                                           >  5:196419/1‑132 (MQ=255)                                                 gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGAACCGACCCGGT‑AAAACC‑GCTCGAACAATcgcgccg                                                                                >  8:398693/1‑133 (MQ=255)                                                        tgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCCCCCACCAGCCCCGACCCCCCCCACCCCACCCGGT‑AAAAAC‑CCTCGCACAAGCGCGACATGTGct                                                                          >  1:338598/1‑135 (MQ=255)                                                        tgctcACCTCGCCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCCCCACCCCCACCCCGGCCCCGCCCACCCGCCGCGACCGCCCCCAACCTACCCGGT‑AAAACC‑GCTCGAACACGCGCGACTGGTGct                                                                          <  1:53382/135‑1 (MQ=255)                                                                                   aCTGCATGATCGCCTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                               <  6:51642/135‑1 (MQ=255)                                                                                            tCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACgg                                         >  4:334564/1‑132 (MQ=255)                                                                                            tCGACTTCACGATCGCCAACCCCAACCCCGTCCCGCCCTCCAGCCGCGAACGCCCCGCACCAACCCCGT‑AAAACC‑CCTCGCACAAGCGCGCCTGGTGCTCTGGAGGCCCCCCCGGCCCCGGGGTCTCCCCgg                                         >  4:40758/1‑132 (MQ=255)                                                                                                  tCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGCACCTACCCCGT‑AAAACC‑GCTCGAACAACCGCGACTGGTGCTCTGCAGGCCCACCCGGACCCTGGTTCTGCACGGATAGAATGa                                >  5:341794/1‑135 (MQ=255)                                                                                                  tCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGa                                >  5:164898/1‑135 (MQ=255)                                                                                                   cccGAGCGCCACCCCCCCCCCGG‑CCCCCCCCCCCCCCCCGCACCCCCCCAACCATCCCGGTAAAAACC‑GCTCCACACCCCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGa                                <  8:141819/132‑1 (MQ=255)                                                                                                              accccaacccGGTCCCGCCCTCCAGCCGCGCACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGCCTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGAGAGAATGACTTCGCTGccc                     >  7:300118/1‑134 (MQ=255)                                                                                                                  caccccGGCCCCCCCCCCCCGCCGCCAACGCCCCCACCCAACCCGTT‑AAAACC‑GCCCCAACACCCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGt               <  7:360006/133‑1 (MQ=255)                                                                                                                   cccccGGCCCCGCCCCCCCGCCGC‑TACGGCCCCAACACTCCCCGT‑AAAACCTGCTCCACCACCCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGTATAGAATGACTTCGCTGCCCGAGCGt               <  1:341618/133‑1 (MQ=255)                                                                                                                         gTCCCGCCCTCCCCCCCCCAACCCCCCGAACCTCCCCGGA‑AAACCC‑GCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATg         <  3:40758/135‑1 (MQ=255)                                                                                                                         gTCCCCCCCTCCCCCCCCCAACGCCCCGACCCAACCTGCA‑AAAACC‑CCTCCAACACCCGCGACTGGAGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGacg         <  5:361053/135‑3 (MQ=255)                                                                                                                             cGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCCCGGATAGAATGACTTCGCTGCCCGCGCGTTCGATGCGCag    >  7:94274/1‑136 (MQ=255)                                                                                                                               ccgtCCTGCCGGGAACGCCCCGAAGCTCCCCTTT‑AAAACG‑GGTCGACCTCGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCagag  <  8:73996/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGT‑AAAACC‑GCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAG  >  ACB122/1630408‑1630669

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: ‑