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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,630,540	0	A	C	27.6%	54.8 / 6.2	47	L419R (CTG→CGG)	czcS	Sensor protein CzcS
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base C (8/6);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAG  >  ACB143/1630408‑1630669                                                                                                                                     |                                                                                                                                  gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCCCGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGcc                                                                                                                                >  6:283072/1‑136 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGcc                                                                                                                                >  5:81884/1‑136 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGcc                                                                                                                                >  3:157701/1‑136 (MQ=255) gcCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGcc                                                                                                                                >  2:38918/1‑136 (MQ=255)   cGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgc                                                                                                                              <  3:293998/136‑1 (MQ=255)   cGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgc                                                                                                                              <  5:311246/136‑1 (MQ=255)      tCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTCCCGCCCTCCCGCcgcgga                                                                                                                           >  1:289789/1‑134 (MQ=255)      tCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGaa                                                                                                                           >  5:167222/1‑136 (MQ=255)      tCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCCGCCTCCCGCcgcgcc                                                                                                                           >  6:837/1‑134 (MQ=255)        aCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                         <  2:64882/136‑1 (MQ=255)        aCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                         <  7:175553/136‑1 (MQ=255)        aCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                         <  1:42747/136‑1 (MQ=255)                ggAAACACCAGGGTAAACCTGGTCTGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGAGCGCCAACCCCAACCCGGGCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                 <  6:306879/136‑1 (MQ=255)                ggAAACACCAGGGTAAACCGGGGAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                 <  4:11762/136‑1 (MQ=255)                ggAAACACCAGGGTAAACCGGGGAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                 <  7:148941/136‑1 (MQ=255)                         aGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCCGTCCCGCCCTCCCGCCGCGCACGCCCCGAACCTAccccg                                                                                                         >  2:83438/1‑133 (MQ=255)                         aGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCCGGCCCGCCCTCCCGCCCCGCACGCCCCCCACCAACCCgg                                                                                                         >  5:256676/1‑135 (MQ=255)                              aaaCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCACCCGCCGCGAACGCCCCCACCCTACCCGGTAAAAc                                                                                                   >  3:310520/1‑136 (MQ=255)                                  cGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCt                                                                                               >  6:165267/1‑136 (MQ=255)                                  cGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCGCCAGCCGCGACCGCCCCGAACCTACCCGGTAAAACCGCt                                                                                               >  2:36418/1‑136 (MQ=255)                                   gggTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCAGCGAACGCCCCGAACCTACCCGGTAAAAACGCTc                                                                                              >  8:234884/1‑136 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGCCCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGAGCGCCAACCCCAAACCGGTCCCGCCCTCCACCCGCGTACGCCCCGAACCTACCCGGAAAAACCCCACGAAc                                                                                          >  1:308131/1‑135 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACCCCCCGCACCGACCCGGTAAAACCGCTCGAAc                                                                                          >  2:254242/1‑135 (MQ=255)                                        ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGACCCTACCCGGTAAAACCGCTCGAAc                                                                                          >  3:204730/1‑135 (MQ=255)                                             cAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAgcgc                                                                                    <  1:24592/136‑1 (MQ=255)                                              aGGGGGCAGCTGCTGCCCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCACCCCGGTCCCGCCCTCCAGCCGCGCACGCCCCCCACCAACCCGGTAAAACCGCGCGAACCACcgcg                                                                                   >  2:62134/1‑136 (MQ=255)                                                 gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGCCt                                                                                >  7:2739/1‑136 (MQ=255)                                                 gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACt                                                                                >  3:322382/1‑136 (MQ=255)                                                   gCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTgg                                                                              <  4:157701/136‑1 (MQ=255)                                                       ctgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGAACCCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGct                                                                          >  3:120836/1‑136 (MQ=255)                                                          ctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCACCCCGGTCCCGCCCTCCAGCCCCCGACCCCCCCCACCAACCCGGTAAAACCGCTCGAACCAGCGCGACTGGTGCTCTg                                                                       >  1:371574/1‑136 (MQ=255)                                                             aCCTCGCCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCGCGAACAAGCGCGACTGGTGCTCTGc                                                                      <  2:126524/134‑1 (MQ=255)                                                                 cgccccgcccgCCATGCAACTGATTTATCCACTTCACGATCGCCAACCCCAACCCGTGCCCGCCCTCCCCCCGCGAACACCCCGAACCTACCCGTTAAAACCGCTCCAACACGCGCGACTGGTGCTCTGCAGGcac                                                                <  1:116931/133‑1 (MQ=255)                                                                   aCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCacac                                                              <  8:62258/136‑1 (MQ=255)                                                                     ccgcccgcCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCGCCCGCCGCGTACGCCCCGAACCTACCCGGTAAAACCGCGCGAACCAGCGCGCCTGGGGCTCTGCCGGCACAccc                                                            >  3:87897/1‑136 (MQ=255)                                                                       gcccgccATGCAACTGAATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCGGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCgg                                                          <  7:234884/136‑1 (MQ=255)                                                                               tGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTAACCGGTAAAACCCCTCGAACAAGCGCCACTGGTGCTCTGCAGGCACACCCGGACCCTGGt                                                  >  8:385608/1‑136 (MQ=255)                                                                                  aaCTTCATGATCGACTTCCCGAACGCCAACCCCAACCCGGTCCCGCCCCCCAGCCCCGAACGCCCCCAACCTACCCGGTAAACCCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                               <  6:256676/136‑1 (MQ=255)                                                                                  aaCTGCATTTTCGACTTCACGATCGCCAACCCCACCCCGGGCCCCCCCTCCCGCCGCTAACGCCCCGAACCTACCCGTTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                               <  3:69421/136‑1 (MQ=255)                                                                                  aaCTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCCAACGCCCCGAACCTCCCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                               <  2:9043/136‑1 (MQ=255)                                                                                              gACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGa                                    <  6:375362/135‑1 (MQ=255)                                                                                                                       cGGTCCCGCCCTCCAGCCGCGCACGCCCCGAACCAACCCCGTAAAACCGCTCGAACAAGCGCGACAGGTGCTCTGCAGGCACACCCGGCCCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGAt          >  3:372949/1‑136 (MQ=255)                                                                                                                         gcccccccCACCCGCCCCCACCGCCCCGCACATACCCGGAAAAACCGCTCCAACAACCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgc        <  2:395293/134‑1 (MQ=255)                                                                                                                         gTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCTAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGGTAGAATGACTTCGCTGCCCGAGCGTTCGATgc        <  8:198814/136‑1 (MQ=255)                                                                                                                             cGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCCCGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCag    >  8:233719/1‑136 (MQ=255)                                                                                                                             cGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgcg       >  4:391383/1‑133 (MQ=255)                                                                                                                               cccTCCCCCCGCGACCGCCCCTAACCTCCCCGTTAAACCCGGTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCagag  <  5:38971/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  GCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCTACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAG  >  ACB143/1630408‑1630669

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