Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB122	3,440,343	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
*	ACB122	3,440,343	0	A	C	76.0%	14.7 / ‑1.6	26	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base C (16/3);  total (22/4)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGA  >  ACB122/3440227‑3440470                                                                                                                     |                                                                                                                                aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGcc                                                                                                             >  1:30989/1‑136 (MQ=255)           gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGgtgt                                                                                                   >  5:248898/1‑137 (MQ=255)           gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtgt                                                                                                   >  6:449134/1‑137 (MQ=255)                     ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCaa                                                                                         >  3:378758/1‑137 (MQ=255)                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCACCCGGGGTGTGGCAGGCCACCGCCTATaa                                                                                >  8:8440/1‑137 (MQ=255)                              cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCAGCCGGGGGGTGGCAGGCCAACGCCTATaa                                                                                >  4:125807/1‑137 (MQ=255)                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCAGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCCCAGCCGATAAg                                                                               >  2:123152/1‑137 (MQ=255)                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCCACCCCGATaaa                                                                               >  2:123148/1‑136 (MQ=255)                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACCCCGATaaa                                                                               >  3:242415/1‑136 (MQ=255)                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGGCACCGCCGATAAg                                                                               >  5:232086/1‑137 (MQ=255)                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAATGCGCAGCTc                                                                      >  6:72931/1‑137 (MQ=255)                                         aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCTGGCCAACGCCGAAAAGGCGCAGCTCg                                                                     >  5:30583/1‑137 (MQ=255)                                           aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGCGCCGCCGGGGTGTGGCACGGCAACGCCGATAATGCGCAgcgcgcg                                                                   >  1:329454/1‑136 (MQ=255)                                            cccgcgcgcTCGTGGTGCGGCAGTTCCGGCTTGGGCCCCCCGGGGGGGCGCGGGGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg                                                                  <  8:194619/132‑1 (MQ=255)                                                   gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                           >  8:136680/1‑137 (MQ=255)                                                   gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGACGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                           >  1:127766/1‑137 (MQ=255)                                                        aGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACGCCGATAAGGCGCCGCTCGCTGATCAGGCCGCTg                                                      >  6:398540/1‑137 (MQ=255)                                                           tgctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCAGCGCGGGGATCAGGCCGCTGAcg                                                   >  4:191546/1‑136 (MQ=255)                                                                        gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCAACCCCGATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGtt                                      >  8:108020/1‑137 (MQ=255)                                                                        gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCCACCCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACGTCGGCAAGCAGtt                                      >  6:225777/1‑137 (MQ=255)                                                                                  cccGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCTGGCGGGGCCCCGGCCCCCGGGGGGTGGGGCGCCCACGCCTATAACGCGCAGCACGCCGCTCAGGCCGCGGTCATCGGCCAGCACTTGGGGGGGaaa                            >  4:349608/1‑136 (MQ=255)                                                                                     gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACTTCGGCAAGCAGTTGGGGGGGAAGGAt                         >  6:182620/1‑137 (MQ=255)                                                                                          ggcgcgGTGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                    <  8:232619/135‑1 (MQ=255)                                                                                                gTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGGCGGCCCACGCCGCTAAGGCGCCGCGCGGGGGTCAGGCCGAGGGCTTGGGCGAGCGGTTGGGGGGGAAGGGTcgcgcgggcgg              >  2:50988/1‑132 (MQ=255)                                                                                                     gccgcggctgGCCGGGGGGGCCGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTgg         <  1:92086/131‑1 (MQ=255)                                                                                                            ggggCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga  <  6:373163/135‑1 (MQ=255)                                                                                                                     |                                                                                                                                AATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGA  >  ACB122/3440227‑3440470

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