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	82.6%	37.9 / ‑1.0	47	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (31/7);  total (35/12)
Fisher's exact test for biased strand distribution p-value = 7.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  ACB122/3440210‑3440468                                                                                                                                      |                                                                                                                              cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAg                                                                                                                              >  2:424871/1‑135 (MQ=255)     cACCAAACGTATCAATGGGGGTAGCACCGGGCTGGCGGATCGCACGGCCCTCTACACGCGAGCGCGCGCGGGGCGGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGGGGCGCGGCTGCGGGCCGCCAgcgcg                                                                                                                          <  8:78552/135‑2 (MQ=255)               ctccttgggggTACCCAGGGTCTGGCGCAGCGCCACTCCCTCTACAAGCGAGCGCGCGCGGGGCTGCCGGGCGGGCGTGGCCCCCCCGCGGGGGCGGCGGGGCGCGGGGGCTGGCCGCCAGCGCGGGCGGCGccc                                                                                                                <  4:447420/130‑1 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCAGCGCGGGAGGCGCCCGGGc                                                                                                            >  6:7800/1‑135 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCgcgc                                                                                                            >  1:217384/1‑135 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGc                                                                                                            >  2:17983/1‑135 (MQ=255)                   aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAACGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGc                                                                                                            >  7:424920/1‑135 (MQ=255)                            gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCggggg                                                                                                   >  4:459883/1‑134 (MQ=255)                                cGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTgg                                                                                               >  3:513905/1‑135 (MQ=255)                                   tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAg                                                                                            >  8:440628/1‑135 (MQ=255)                                       gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGccc                                                                                        >  3:177855/1‑134 (MQ=255)                                        aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCaa                                                                                       >  4:529901/1‑135 (MQ=255)                                               cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCTCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCCGGCGCCCGGGGGGTGGGGGGCCCCCGccga                                                                                 >  5:15236/1‑132 (MQ=255)                                               cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTCGCCGCGCCCGGGCCGCCCGGGGGTGGGAGGGCCACCCCtc                                                                                 >  6:164115/1‑133 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCAATGGCGGCGCCCGGGCAGCCGCGGGGTGTGAGGCCACCGCCaata                                                                               >  7:453102/1‑135 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCACCGccgc                                                                                 >  3:193923/1‑131 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCGGGGGGTGGGGGGCCCACGccga                                                                                 >  5:262843/1‑131 (MQ=255)                                                aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCGGCCGGGGTGGGGGAGGCCCACGccga                                                                                 >  7:372088/1‑131 (MQ=255)                                                    cccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGgc                                                                           >  3:469562/1‑135 (MQ=255)                                                       tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCTGGGGGGGGCAGGCCACCGCCGATAAGGCGCa                                                                        >  5:239147/1‑135 (MQ=255)                                                         tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCCCCCCGATAAGGcgccgc                                                                      >  2:184406/1‑135 (MQ=255)                                                          aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCAGGCCGCGCCCGGGCCGCCTGGGGGTGGCAGGCCCACGCCGAAAAGGcgccgcg                                                                     >  1:380902/1‑134 (MQ=255)                                                          aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCCGGCGGCGCCCCGGCAGCCGGGGGGTGGGCGGCCCACCCCGATAACGCGCAGCt                                                                     >  6:104661/1‑135 (MQ=255)                                                              gcgagcgCGCGCGGGGCTGCCGTGCGCGCGGGGGCCCCCCGGCTGGGCGCCGGGGGGCTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                 <  3:208603/135‑1 (MQ=255)                                                                 agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCCCCCCCCATAAGTCGCCGCTCGCGGa                                                               >  6:238570/1‑134 (MQ=255)                                                                 agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCGGGCCCACGCCCATAAGGCGCACCgcgcggc                                                               >  3:225383/1‑131 (MQ=255)                                                                 agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCAGGCCCACCCCGATAAGGCGCCGCGCGCGGa                                                               >  8:103909/1‑134 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCCCGCCGCTAAGGCGCAGCGCGCTGCTCAGGc                                                         >  7:331411/1‑135 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGCTCAGGc                                                         >  7:362066/1‑135 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCCGGCCGCCTCGGTGGGGGCGGCCAACCCCGATAAGGCGCCGCGCGCGGATccgcg                                                         >  4:153969/1‑131 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCCCTGGCGGCGCCCCGCCGGCCCGGGTGGGGCCGGCCCACGCCCATAAACCGCACCACCCGGATCacaa                                                         >  6:310377/1‑132 (MQ=255)                                                                      tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGGGGGCGGCCCACGCCGAAAAGGCGCCGCGCGCGGATCCGGc                                                         >  7:289145/1‑135 (MQ=255)                                                                        gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGGCCGGGCCGCCGGGGTGTGGGGGGCCCACGCCGATAAGGCGCCGCGCGCGGCTCAGGCCg                                                       >  1:444003/1‑135 (MQ=255)                                                                                 cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCCCCCGGGCCGCCGGGGGGTGGGCGGCCACCCCCGAAAACGCGCACCGCGCTGCTCAGGCCGATGACATCg                                              >  8:313013/1‑135 (MQ=255)                                                                                         gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCGGGCGGCGCCCGGGGCGCCGGGGGGTGGGCGGCCCCCCCCGCTAAAGCGCCGCGCGCGGGTCAGGCCGCGGACTTCGGGcagcag                                      >  6:160476/1‑135 (MQ=255)                                                                                         gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCCCGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCTGACGGCGGCGAGCAg                                      >  7:411501/1‑135 (MQ=255)                                                                                          cgtGGGCCCCCCGCCGGGGCGCCGGGGCGCTGCGGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                     <  7:69764/133‑1 (MQ=255)                                                                                                             ggcggtggcggtggcgCCGCCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:7800/131‑1 (MQ=255)                                                                                                             cgcgggggCGCGGCGGCGGGCCGCGCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:200640/129‑1 (MQ=255)                                                                                                              gCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGGGGCGCCCGGGCCGCCGGGGGGGGGCCGGGCCCCGCCGCTAAGGCGCCGCTCGCGGATCAGGCCGCTGACTTCGGCGCGCCGGTGGGGGCGAAGGCTcgcgcg                 >  8:285610/1‑135 (MQ=255)                                                                                                               gggTGGGGCTGCGGCGGGCCGCCAGCGCGGCCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCgcgcgc                <  1:76311/133‑1 (MQ=255)                                                                                                               ccGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGGCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCCGCTCGCGGGTCAGGCCGATGCCGTGGGCGAGCAGTTGGGGGAGAAGGCTCgcgcgc                >  6:482279/1‑135 (MQ=255)                                                                                                                       ctgctgctgGCCGCCAGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg        <  8:98568/135‑1 (MQ=255)                                                                                                                            gctgGCCGCCCGCGCCGCCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg   <  1:204726/135‑1 (MQ=255)                                                                                                                             gtgGCCGCCCGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  5:104661/134‑1 (MQ=255)                                                                                                                             ggggCCGCCGGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  2:89471/133‑1 (MQ=255)                                                                                                                             cgggCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  7:103909/133‑1 (MQ=255)                                                                                                                                      |                                                                                                                              CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  ACB122/3440210‑3440468

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