Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	841,643	G→T	22.5%	H609N (CAT→AAT)	CJ019_00805 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	841,643	0	G	T	22.5%	74.7 / 2.6	49	H609N (CAT→AAT)	CJ019_00805	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/16);  new base T (2/10);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 9.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCCGACGATCACACCGG  >  ACB143/841510‑841774                                                                                                                                      |                                                                                                                                    aaCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaa                                                                                                                                   <  7:85266/136‑1 (MQ=255) aaCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaa                                                                                                                                   <  2:91341/136‑1 (MQ=255)    gAGGCGGCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGt                                                                                                                                <  6:88330/136‑1 (MQ=255)         gTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTaaa                                                                                                                             >  2:183714/1‑134 (MQ=255)         gTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAACCCAGCAGGCCGGTTACGGTGGTGGCGCCGAAAACCCCGTGTTTCATTGGAACCACGGATTACAAGGGATGAAAGTaaaaa                                                                                                                           >  7:283984/1‑134 (MQ=255)            gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGt                                                                                                                         <  7:39880/135‑1 (MQ=255)               cTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTg                                                                                                                     <  6:30520/136‑1 (MQ=255)                     aatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGccc                                                                                                                <  2:201598/135‑1 (MQ=255)                     aatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGccc                                                                                                                <  6:250183/135‑1 (MQ=255)                     aatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGccc                                                                                                                <  5:25709/135‑1 (MQ=255)                          aGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCCGCAGGCCGCTTTCGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAACAGGGATGAAAGTAAATTGTGCTGAGCCCGAAccg                                                                                                          >  8:9205/1‑134 (MQ=255)                          aGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGAAGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATATTGCTGAGCCCGaa                                                                                                             >  5:101570/1‑133 (MQ=255)                            gTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCa                                                                                                        <  1:253014/136‑1 (MQ=255)                               aCGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTTAGGTGCTGTCGCCGAAAACCCCGTGTTTGATTGGAACCGCGGATTACCAGGGCTTAAAGTAAATAGTTCTGAGCACGCACAGCaatt                                                                                                     >  8:114730/1‑133 (MQ=255)                                     ccgCCATGCCGTTTGGCTGCCTCAAAAGCCAGAAGCCCGCTAACGGGGCGGCCGACTAAAACCCGGTGTTTCATGGGAACCTCGGATTCCCAGGGATCAAAGAAAATAGTGCTGAGCCCGAACAGCTTTGCGaacc                                                                                               <  2:261033/134‑3 (MQ=255)                                     cAGCCAAGCCGTTTGGCTGCCTCAAAAGCCCGCAGGCCGCTTACGGTGGTGTCGCCCAAAACCCCGTGTTTGATTGGAACCTCGGATTTCCAGGGATTGAAGTAAATATTGCTGCGCCCGCAACGCATTGCGaagc                                                                                               >  4:5697/1‑136 (MQ=255)                                     cAGCCAAGCCGGTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCGGTCGCCGAAAACCCCGTGGTTCATTGGAACCTCGGATTACCAGGGTTGAAAGTAAATATTGCTGAGCCCGGACAGCATTGCGaatc                                                                                               >  4:1421/1‑134 (MQ=255)                                       gCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACTAAAACCACGTGTTTCATTGGAACCTCGTATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCTAACAGCATTGCGaagcaa                                                                                             <  2:91291/136‑1 (MQ=255)                                                  gggCTGCCTCAAAACCCCGCAGGCCGCTTACGGTGCTGTCGACTAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTCAGCCCGAACACCATTGCTAAGAAAGCGACCGTCg                                                                                   <  1:81668/134‑1 (MQ=255)                                                        ccTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCGGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTgc                                                                            <  2:170001/136‑1 (MQ=255)                                                         ctaaaacGCCCGCAGGCCGGTTCGGGTGCTGTCGACTAAAACCAGGTTTTACATTGCAACCTCGTTTTACCAGGGATGAAATTAAATAGTGCTCCCCCCCAACAGCATTGCCAACCAAGCGACCGTCGAAGGTgcg                                                                           <  2:276914/133‑1 (MQ=255)                                                         cTCAAAAGCCAGCAGGCCGATTACGGTGCTGTCGACTAAAACCACGTTTTTCATTGGAACCTCGGATAACCAGGGATTAAAGTAAATAGTGCTGAGCCCTAACACCATGGCGAAGCAAGCGACCGTCGAAGGTgcg                                                                           <  8:62645/136‑1 (MQ=255)                                                         cTCAAAACCCAGCAGGCCGCTTAGGGTGCTGTCGACTAAAACCACGTGTTTCATTGTAACCTCGTATAACCTGGGATCAAAGTAAATAGTGGTGAGCCCTAACAGAATTGCTACGAAAGCGACCGTCGAAGGTgcg                                                                           <  2:71297/136‑1 (MQ=255)                                                           cAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACCTGTTTCATTGGAACCTCCGATTAACAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCAATGCGAAGCAAGCGACCGTCGAAGGTgcgc                                                                          >  6:212484/1‑135 (MQ=255)                                                                      agGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGaa                                                              >  4:261877/1‑136 (MQ=255)                                                                      agGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGaa                                                              >  6:91641/1‑136 (MQ=255)                                                                         ccGCTTACGGTGCTGTCGACCAAAACCACGTTTTTCATTGCACCCTCGTATTACCAGGTATTAAATAAAATAGTGCTCCGCCCGAACACCTTGGCTAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATtc                                                           <  6:82081/136‑1 (MQ=255)                                                                               cgggggCTGTCGACTAAACCCAGTTGTTTAATTGCAACCTCGTATTACCAGGTATTAATTAAAATTGTGCTGAGCCCTAACAGCATTGCCAAGAAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAg                                                      <  8:283984/130‑1 (MQ=255)                                                                               agggTGCTGTCTACTAAAACAACGTGTTTATTTGCAACCTGGTTTAACCTGGGATTAAAGTAAATAGTGCTGACCCCGAACACATTTGCTAAGAAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAg                                                      <  6:123013/133‑1 (MQ=255)                                                                               agggTGCGTTCTACTAAAACCAGGTGTTTCATTGCAACCTCTTATAACCTGGGATTAATTTAAATTGTGGTGAGCCCTACCTGATTTGCTATCCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAg                                                      <  6:238571/133‑1 (MQ=255)                                                                               aCGGGGCGGGCGACGAAACCCACGTGTTTCATTGTAACCTCGTATAACCAGGTATCAAAGAAAATAGGGCTGAGCCCGAACACATTTCCTACGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAg                                                      <  6:103785/135‑1 (MQ=255)                                                                                cGGGGCTGTCGACGAAAACCACGTTTTTCATGGGAACCTCGTATTACCAGGCATGAATGAAAATCGCGCTGCCCCCGAACAGCATTGCTACGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCt                                                    <  4:120646/136‑1 (MQ=255)                                                                                   tGCTGTCTACGAAAACCACGTGTTTCATTGTAACCTCGTATTACCAGGGATTAAAGTAAATAGTGCTGAGCCCGAACAGAATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGt                                                 <  5:310493/136‑1 (MQ=255)                                                                                       gTCCACCAAAACCACGTGTTACATTGTAACCGCGTATACCCCGGGATCAAAGTAAATAGGGCTGCCCCCCAACCCCATGGCGACCCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGAcc                                             <  4:205327/136‑1 (MQ=255)                                                                                        tCGACGAAACCCCCGTTTTTAATGGGAACCTGGGATTCCCTGGTATTAATGTAAATCTTCGCCCGCCCCAACTCCTTTGGGAAGAAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGAcc                                             <  2:100062/135‑1 (MQ=255)                                                                                        tCGACGAAAACCACGTGTTTCATTGGACCCTCGGATAACCAGGGATTAAAGTAAATAGTGCTGAGCCCGAACAGCTTTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCt                                            <  1:51144/136‑1 (MQ=255)                                                                                         cgacgaAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTg                                           >  8:260037/1‑136 (MQ=255)                                                                                             gaAAACCACGTGTTTCATTGTAACCTCGGATTACCAGGTATTAAAGTAAATTGTGCTGAGCCCGAACAGAATTGCTAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTggagg                                       <  2:81434/136‑1 (MQ=255)                                                                                              aaaaCCACGTGTTTCTTTGTAACCTCGGATTACACGGTATTAAATTAATTAGTGCTGAGCCCTAACAGCTTTGCGAAGAAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTggagg                                       <  6:101570/135‑1 (MQ=255)                                                                                                  ccACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGcc                                  >  3:193247/1‑136 (MQ=255)                                                                                                     cGTGTTTCATTGGAACCTCGGATAACCAGGTATCAAAGTAAAACGTGCTCCGCCCTACCAGCATGGCCAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGAc                               <  7:26199/136‑1 (MQ=255)                                                                                                        gtTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAg                            >  8:23579/1‑136 (MQ=255)                                                                                                                gcacccTCGCATAACCAGGGATCAAATAAAATAGTGCTGCCCCCCAACACAATGGCCAAGCACGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCAtt                     <  1:63339/131‑1 (MQ=255)                                                                                                                 gAACCTCGGATTACCAGGTATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGAGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAg                   <  6:57882/136‑1 (MQ=255)                                                                                                                 caaCCTCGCATTCCCCGGTATCAAATAAAATCGTGCTCCGCCCGAACACCATTGCCACCAACGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAg                   <  1:227657/135‑1 (MQ=255)                                                                                                                        ggATTACCAGGGAGGAAAGTAAATAGTGCTGCGCCCGAACCGCTTGGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCcgac             <  7:142656/135‑1 (MQ=255)                                                                                                                          aTTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCCGACGAt          >  7:24317/1‑136 (MQ=255)                                                                                                                          aTTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCCGACGAt          >  8:3030/1‑136 (MQ=255)                                                                                                                                  ggATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCCGACGATCACACCgg  >  4:106443/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                                    AACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATACTTGGGAATTCTCAAGCTGGTGACCTGGAGGATGCCGACCAGGGCCATTAGGCCGACGATCACACCGG  >  ACB143/841510‑841774

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