Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	3,440,326	A→C	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	80.0%	21.9 / ‑2.6	31	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base C (21/3);  total (24/7)
Fisher's exact test for biased strand distribution p-value = 7.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG  >  ACB143/3440193‑3440458                                                                                                                                      |                                                                                                                                     cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                                    >  4:119732/1‑136 (MQ=255)                  aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCgcgc                                                                                                                   >  7:294052/1‑136 (MQ=255)                            gCAACGGTCTGGCGGCTCGCCAGGCCCTCTACACGCCTGCGCTCGCGGTGCTGCTGTGCGCGCGGGGGCCCCCCGCCGGGGCGCCGGGGCGCGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                         <  8:220117/136‑1 (MQ=255)                            gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                         >  3:368715/1‑136 (MQ=255)                            gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                         >  8:108816/1‑136 (MQ=255)                            gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGgtg                                                                                                         >  4:237919/1‑136 (MQ=255)                                   tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAgg                                                                                                  >  8:101852/1‑136 (MQ=255)                                    cgggaggcGCGCAAGGCCCTCTACACGCGCGCGCGCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGGGGCGCTGCTGCTGGCCGGCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                 <  1:109158/134‑1 (MQ=255)                                     gggaggaGCGCAAGGCCCTCTACCCGCGCGCGCGCGTGGTGCTGCAGGGCGCGCTGGGGCCACCCGCCGGGGCGCCGTGGCGCGGCGGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                <  8:43002/135‑1 (MQ=255)                                        aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCAGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCAc                                                                                             >  7:303069/1‑136 (MQ=255)                                               cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCTGGCGGCGCCCCGGCCGCCTCGGGGTTGGAGGGCACCGccct                                                                                        >  2:365452/1‑132 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCACCGCCGATaa                                                                                     >  1:271519/1‑136 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACCCCTATaa                                                                                     >  1:273858/1‑136 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCCGGCGGCGCCCCGGCCGCCTAGGTGTGGCAGGCCCACCCCCATaa                                                                                     >  6:15393/1‑136 (MQ=255)                                                              gcgcgcgcTCGAGGGGCTGCAGGGCGCGCGTGGGCCACCCGCCGGGGCGCCGGGGCGCTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg                                                                       <  8:48649/136‑1 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCCCCGGGGTGTGGCCGGCCCCCGCCGATAATGCGCCGCTCGCTGTTAAgg                                                                 >  2:43999/1‑136 (MQ=255)                                                                    gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCCCCCCGGCCGCCGGGGTGTGGCAGGCCACCGCCGATAAGTCGCAGCGCGCGGATCAgg                                                                 >  1:164783/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCCGGCCCCCGCCTATAAGGCGCAGCGCGCGGATCAGGcc                                                               >  3:298449/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGGAGGCCAACGCCGAAAATGCGCAGCTCGCGGCTCAGGcc                                                               >  3:283915/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGAGGCCACCGCCGATAAGTCGCCGCTCGCTGATCAGGcc                                                               >  7:130206/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTCGCGGATCAGGcc                                                               >  3:227235/1‑136 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTGGCGGCGCCCGGGCCCCCGGGGGGTGGGAGGCCCCCGCCTATAAAGCGCCGCTCGCGGATCAGGcc                                                               >  3:262133/1‑136 (MQ=255)                                                                                         gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCTATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                            >  7:291367/1‑136 (MQ=255)                                                                                         gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCGATAAGGCGCCGCTCGCGGCTCAGGCCGATGACGTCGGCCAGCAGt                                            >  3:31444/1‑136 (MQ=255)                                                                                           ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCGCCCGCCGATAAGGCGCAGCGCGCGGATCTGGCCGAGGCCTTCGGGAGGCAGTTg                                          >  5:26988/1‑136 (MQ=255)                                                                                                ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGGGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                       >  7:189168/1‑134 (MQ=255)                                                                                                            ggcccGTGGCGGGGGTGCTCGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                         <  8:265412/133‑1 (MQ=255)                                                                                                            gcgcgGTGGCGCTGGTGGTCGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                         <  4:177707/135‑1 (MQ=255)                                                                                                             cGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgcg                        <  7:16844/136‑1 (MQ=255)                                                                                                                   ggCGCTGCTGCTGGCCGACAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGt                  >  7:89984/1‑136 (MQ=255)                                                                                                                                 ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa    <  6:121672/136‑1 (MQ=255)                                                                                                                                  cgccCGCGCTGGCGGCGCCCGGGCACCCGGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTCGCTGATCAGGCCGATGACGTAGGGAAGCCGTTGGCGGCGAAGGATAGCGCGCACCGGCGGGAGCGGGAGAAg   >  2:84876/1‑136 (MQ=255)                                                                                                                                   gccACCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCTATAATGCGCAGCTCGCGGATCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTGGTGCGGGAGAAgg  >  4:160643/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                                     CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG  >  ACB143/3440193‑3440458

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