Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,684,652	G→C	22.2%	K124N (AAG→AAC)	CJ019_05336 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,684,652	0	G	C	22.2%	95.7 / 3.0	53	K124N (AAG→AAC)	CJ019_05336	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/27);  new base C (6/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGATCGCCAG  >  ACB122/5684521‑5684783                                                                                                                                    |                                                                                                                                    aGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAATATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCg                                                                                                                                  <  8:59067/135‑1 (MQ=255) aGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCg                                                                                                                                  <  1:72422/135‑1 (MQ=255)      cACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCCCCCCCCACAAAATCGGCAgt                                                                                                                             >  2:104393/1‑135 (MQ=255)              aCAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTgac                                                                                                                       <  4:175916/133‑1 (MQ=255)              aCAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTgac                                                                                                                       <  7:47887/133‑1 (MQ=255)              aCAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGGGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTgac                                                                                                                       <  3:194667/133‑1 (MQ=255)                      gCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGTACAGCAAGATCAACAAGATCACCCCACCCCCAGCGGCCAGGCCCGAGGTGGTGCCCCCCCCGCCCCCGCAAAAAGAGCCACTGGTTCCGACCGCGcc                                                                                                              >  4:124428/1‑134 (MQ=255)                      gCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCCCCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCCCCAAAGTCGGGAGTGGTGGCGCCGGCGccc                                                                                                             >  2:233797/1‑135 (MQ=255)                      gCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGccc                                                                                                             <  3:191892/135‑1 (MQ=255)                      gCGGGGCCCAAGACCACTGAGCGGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCCCCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGccc                                                                                                             <  1:104393/135‑1 (MQ=255)                          tGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGCCCCGAGGTGGTCCCCCCCCCTCCCCCGAAAAACGCGGCCGGGGGGGCGACGGCGCCCACGc                                                                                                         >  2:303566/1‑135 (MQ=255)                             ccAAGACCACTCCGCTGGCGCCGTTCCAGGACAGCAAGAGCACCAAGATCACCCCCCCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGc                                                                                                      <  1:289022/135‑1 (MQ=255)                               aaGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTGCCCCCCCCACCCCCGCAAAAATCGGCAGGGGTGGCGACGGCGCCCCACCCGcc                                                                                                     >  5:174238/1‑133 (MQ=255)                                  aCCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAg                                                                                                 <  4:188019/135‑1 (MQ=255)                                       tGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCCGGCCCGAGGTGGTCCCACCCCCACCCCCGCCAAAGACGGCAGTTGTGACGACGGCGCCCAAGCCCCCCAAGGTCgc                                                                                            >  6:135388/1‑134 (MQ=255)                                       tGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCCAAAAACTCGGCAGTGGTGACGACGGCGCCCCAGCCGCCCCAGGTCgc                                                                                            >  7:174168/1‑134 (MQ=255)                                         aGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCCAAAAGTCGGCAGGGGGGCCGGCGGCGCCCCAGCCCCCCCAGGTCgcgc                                                                                          >  6:34686/1‑135 (MQ=255)                                         aGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCGAAAAACTAGGCAGTGGTGCCGACGGCGCCCAACCCCCCCAACGTCggcg                                                                                          >  4:169897/1‑132 (MQ=255)                                         aGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCACGAGGTGGTCCCCCCCCCGACCCCGAAAAACTCGGGAGGGGTGGCGACCGCGGCACCGCCGCACAACGGCGAGc                                                                                          >  3:192664/1‑135 (MQ=255)                                              gcgcCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGAAAAAGTCGGCAGTGGTGTCGCCGGCGCCCCAGCCGCCCAAGGGCGAGCCCcaa                                                                                     >  2:719/1‑135 (MQ=255)                                              gcgcCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAAGTCGGCAGTGGGGCCGACGGCCCCCCCGCCGCCCAAGGTCGAGCCCaaa                                                                                     >  3:96693/1‑135 (MQ=255)                                                        aGGACAGCAAGATCAACAAGCTCACCCCACCGCCGGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGaaa                                                                           <  1:115852/135‑1 (MQ=255)                                                            cAGCAAGATCAACAAGATCACCCCCCCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCaa                                                                       <  7:80024/135‑1 (MQ=255)                                                             aGCAATATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                                                      <  7:113547/135‑1 (MQ=255)                                                             aGCAAGATCAACACGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                                                      <  1:43222/135‑1 (MQ=255)                                                             aGCAAGATCAACAAGAGCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAACTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                                                      <  6:108808/135‑1 (MQ=255)                                                             aGCAAGATCAAAAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                                                      <  7:194894/135‑1 (MQ=255)                                                                 aGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGAAAAACGACGCAGGGGTGACGCCGGCGCCCCAGCCGCACCACGGCCACCCCAAACCCAACGAAAGCCCGgcgc                                                                  >  5:203942/1‑135 (MQ=255)                                                                 aGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCGAAAAAGTCGGGAGGGGTGCCCCCCGCGCCCAAGCCGCCCAAGGGCGCGCCCAAACCCAAAGAAAGCCAGgc                                                                    >  8:78757/1‑133 (MQ=255)                                                                 aGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAACACGGCACTGGGGACGACGGCGCCCAAGCCGCACAATGTAGAGCCCAAACCCAATGAAAGccacgc                                                                    >  2:223796/1‑128 (MQ=255)                                                                  gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGCAAAAGTCGGCCGTGGTGCCGCCGGCGCCCCCGCCGCCCCAGGTCGCGCCCCAACCCCAGGAAAGCCACgcgcc                                                                 >  3:44813/1‑134 (MQ=255)                                                                  gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAAGTCCGCAGTGGTTGCGCCGGCGCCCCAGCCGCCCCAGGTCGGGCCCCCACCCAACCAAAGCAACgcgcc                                                                 >  5:274379/1‑134 (MQ=255)                                                                       acaaGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCGAAAAAATAGGGAGTGGTGACGACCGCGCCCCAGCCCCACAACGTGGCGCCCAAACCCAATGAAACCAAGGCGCAGCCCa                                                            >  4:318583/1‑135 (MQ=255)                                                                                ccccACCCAGTGCTGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAACTCTGCAGTGGAGACGACGGCGCCCAAGACGCACAAGTCCGAGCCCAAACCCAAGGAAACCAAGGCGCAGCCCAAGCCAGCcg                                                   <  7:261381/135‑1 (MQ=255)                                                                                      cGCCGGCGCCCAGCCCCGAGGGGGCCCCACCCCCGCCCCCGAAACAGTCGGCAGTGGGGACGACGGCGCCCAAGCCGCACAAGGTCGAGCACAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa                                             <  2:318887/135‑1 (MQ=255)                                                                                        ccggCGGCCAGCCCCCAGGGGGTCCCACCCCCTCCCCCGAAAACGGCGGCAGGGGTGACGACGGCGCCCAAGCCGCACAAGGGCGCGGCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACg                                           <  7:75817/132‑1 (MQ=255)                                                                                        ccTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACg                                           <  8:91821/135‑1 (MQ=255)                                                                                        ccTGCGGCCAGGCCCGAGGTGGGCCCACCCCCTCCCCCGAAAACGTCGGCAGTGGTTACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACg                                           <  7:350136/135‑1 (MQ=255)                                                                                         cTGCGGCCAGGCCCGAGGTGGTCCCACCCCCACCCCCCCAAAAATAGGCAGTGGGGACCCCGGCGCCCCACCCGCACAAGGTCGAGACCCAACCCCAGGAAAGCCAGGCGCAGCCCAAACCCCCCGCGccccccc                                          >  8:229608/1‑130 (MQ=255)                                                                                                cAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACtt                                   <  6:174238/135‑1 (MQ=255)                                                                                                         ggtggtCCCACCCCCTCCCCCGAAAACGGCGGCAGGGGGGACGACGGCGCCCAAGCCGCACACGGGCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTc                          <  7:350991/135‑1 (MQ=255)                                                                                                               ccccccccCTCCCCCTAAAACGTCGGCAGTGGGGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGACAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCt                    <  5:294300/135‑1 (MQ=255)                                                                                                               cccACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCt                    <  7:229608/135‑1 (MQ=255)                                                                                                                ccccccccTCCCCCGAAAAAGTGGGCAGGGGTGACGACGGCGCCCAAGCCCCACACCGTCGCGCCCCAACCCAAGGAAACCAAGGGCCCGCCCAACCCCGCCCCCCCCCCCCCGGGCTTAGACCGCGCGCCGCTg                   >  6:87653/4‑135 (MQ=255)                                                                                                                 cACCCCCTCCCCCGAAAAAGTCGGCAGGGGGGCCGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAACCCCACGGAACGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGt                  <  5:62023/135‑1 (MQ=255)                                                                                                                          ccccGAAAAAGTCGGGAGTGGTGACGACGGCGCCCAACCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAAGCGCAGCCCAAACCACCCGCGCCGACCCCGGACATAGACAGCTCGCAGCTGTACAGCCAGa         >  6:135639/1‑135 (MQ=255)                                                                                                                          ccccGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGa         <  4:115631/135‑1 (MQ=255)                                                                                                                          ccccCAAAAAGTCGGCAGTGGTGACGCCGGCCCCCAACCCGCCCCAGGTCGCGCCCAAACCCAACGAAAGCCACGCGCAGCCCCAGCCCGCCCCGCCCCCGCCCGACTTCGACAGCTCGCAGCGGTcccgcccgc         >  4:344534/1‑132 (MQ=255)                                                                                                                           cccGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGAt        <  4:192664/135‑1 (MQ=255)                                                                                                                                aaaaGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGATCGCCa   <  7:232186/135‑1 (MQ=255)                                                                                                                                 aaaGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGATCGCCAg  >  1:35854/1‑135 (MQ=255)                                                                                                                                   aGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCa                                                                 <  2:66724/70‑1 (MQ=255)                                                                                                                                   aGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCa                                                                 >  1:66724/1‑70 (MQ=255)                                                                                                                                    |                                                                                                                                    AGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGATCGCCAG  >  ACB122/5684521‑5684783

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