Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	3,440,326	A→C	49.7%	S18R (AGC→CGC)	CJ019_03266 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	3,440,326	0	A	C	49.7%	‑5.8 / 13.4	36	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/5);  new base C (13/5);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB143/3440193‑3440456                                                                                                                                      |                                                                                                                                   cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                                  >  7:229403/1‑136 (MQ=255) cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGCCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGCCCACCCGCCTGGTCGCCGTGGCGCTGCTCCTGGCCCCCAgc                                                                                                                                  >  5:28805/1‑136 (MQ=255)  gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCcgcg                                                                                                                                 >  6:281740/1‑136 (MQ=255)       ccAAACTTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGggc                                                                                                                            >  2:310291/1‑136 (MQ=255)                  aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCAGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCAGGCGGCGCCCCGGc                                                                                                                 >  6:62071/1‑136 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGcc                                                                                                                >  8:180732/1‑135 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCAGGCGGCGCCCGGGCa                                                                                                                >  3:272812/1‑136 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCGGGcc                                                                                                                >  2:84770/1‑135 (MQ=255)                     ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAgg                                                                                                              >  7:274466/1‑135 (MQ=255)                            gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGgtg                                                                                                       >  5:159620/1‑136 (MQ=255)                                   tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGGCGCCTGGGTGTGGCAgg                                                                                                >  6:280034/1‑136 (MQ=255)                                      ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                              >  6:155118/1‑135 (MQ=255)                                      ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCa                                                                                             >  7:51038/1‑136 (MQ=255)                                                aaGGCACTCTACAAGCGAGCGCTCGAGGTGATGCAGTGAGCGCTTGGGCCACCCGCCAGGTCGCCGTGGCGCTGCTGCTGGCCGCCAACGCTGGCCGCCCCCGCGCCCACTCGGAGTGGCAGGCCCACTCCtat                                                                                     >  7:296934/1‑134 (MQ=255)                                                       tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAg                                                                            >  7:315608/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                               >  5:340532/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGGCGGCCGGGCCCGGGCCGCCGGGGTGTGGCAGGGCCACGCCGAAAATGGGCAGCTCGCTGATCAgg                                                               >  6:116700/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGCGGGCCACCCCGATAAGGCGGCGCTCGCGGCTCAgg                                                               >  2:197861/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGGGGCGGCCGGGCAGCCTGGGTGTGGGAGGCCCACGACAATAATGGGCCGCACGCAGATAAgg                                                               >  2:133660/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGCGCCGCCTGGGGGTGGCAGGCCACCCCCGATAAGGCGCCGCGCGGTGATCAgg                                                               >  5:138979/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGCGCCGCCGGGGTGTGGCAGGCCACCGCCGCAAAGGCGCAGCTCGCGGATCAggac                                                             >  5:156083/1‑134 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACCCCGATAAGGCGCCGCTCGCGGGTCAGGcc                                                             >  8:375355/1‑136 (MQ=255)                                                                                 cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTTTGGCAGGGCAACGCCTATAAGGCGCCGCGCGCTGATAATGCCGATTTCTTCgg                                                  >  2:114333/1‑136 (MQ=255)                                                                                         gggggggCCCCCCGCGGGGGCGGGGGGGCGCGGCGGGGGGCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                          <  2:256244/132‑1 (MQ=255)                                                                                         gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGTCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                          >  6:132193/1‑136 (MQ=255)                                                                                           ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGAAGACTTCGGCAAGCAGTTg                                        >  7:214264/1‑136 (MQ=255)                                                                                                    ccGCCTGGGCGCGGGGGCGCGGCGGCGGGCCCCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAg                               <  1:270476/136‑1 (MQ=255)                                                                                                      gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCCACCCCTATAAGGCGCCGCTCGGGGATCAGGCCCATGTCGACGGCAAGCAGATGGCGGCGAAGGa                             >  4:73899/1‑136 (MQ=255)                                                                                                           ggggcggTGGCGCTGCTGCTGGCCGCCTGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                        <  1:332966/130‑1 (MQ=255)                                                                                                           ggcggggTGGCGCTGCTGGGGGCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                        <  8:315608/130‑1 (MQ=255)                                                                                                            tggCCGTGGGGGTGGGGCCGGCCGCCCGCGCGCCCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:202067/134‑1 (MQ=255)                                                                                                            tCGCCGTGGCGGTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  7:239753/136‑1 (MQ=255)                                                                                                            gcggggTGGGGCGGCTGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:66969/131‑1 (MQ=255)                                                                                                            gcgcgGTGGCGGGGCGGCGGGCCGCCAGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  8:259868/135‑1 (MQ=255)                                                                                                            gcgcgGTGGCGCTGCGGGGGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:44413/135‑1 (MQ=255)                                                                                                                 gggcggctgcggctGTCCTCCTGTGGTGGCGGCGCCCGGGCAGCCTGGGTGTTTCAGGCCAACTCCTATAAGGCGCAGCTAGCTGATCATGCCGATGAATTCGGAAAGCAGTTGGATGAGAATGATCGCGCGCACg                  <  8:335305/131‑1 (MQ=255)                                                                                                                              gggCCGCGCGGGCGGGCGCCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga     <  7:90022/135‑1 (MQ=255)                                                                                                                                 ccgccAGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  5:309495/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB143/3440193‑3440456

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