Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,944,425	T→G	28.3%	E86A (GAG→GCG)	CJ019_03729 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,944,425	0	T	G	28.3%	63.8 / 4.2	50	E86A (GAG→GCG)	CJ019_03729	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base G (4/11);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGGTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGCGGGG  >  ACB122/3944292‑3944547                                                                                                                                      |                                                                                                                             gcgcTCGTCGGCCACTGTCAGTTCGAAGCGTGGCCTGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCGCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGGGTAGGCCGGATAGCGGCGGCCGTTGGCCTc                                                                                                                             <  7:162471/135‑1 (MQ=255) gcgcGGGTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTc                                                                                                                             <  1:365180/135‑1 (MQ=255)  cgcTGGTCGGCCTAGCTGAGTGCGAAGCGTGGCATGGCCGGGGCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGATGCTGTCGGGGTAGGCCGGATAGCGGCGGCCGTTGGCGTc                                                                                                                             <  1:139363/134‑1 (MQ=255)  cgcTGGTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTc                                                                                                                             <  1:321835/134‑1 (MQ=255)    cTGGTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTcgcg                                                                                                                          >  6:138241/1‑135 (MQ=255)       gTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCgg                                                                                                                       >  5:75652/1‑135 (MQ=255)             cAAGGTCAGTTCGAAGCGTGGCATGGCCGGGGCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTg                                                                                                                 <  5:155081/135‑1 (MQ=255)               aGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGAc                                                                                                               >  7:361020/1‑135 (MQ=255)                    aGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCAGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGt                                                                                                          >  7:213900/1‑135 (MQ=255)                        cGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGGCCGAGTGCCa                                                                                                      >  6:44138/1‑135 (MQ=255)                         gAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCa                                                                                                      <  7:208761/134‑1 (MQ=255)                         gAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAg                                                                                                     >  8:163014/1‑135 (MQ=255)                         gAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAg                                                                                                     >  7:269847/1‑135 (MQ=255)                         gAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCGCGCGCGGGCCCGGCCCGAGTGGCAg                                                                                                     >  3:303/1‑135 (MQ=255)                         gAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCCTTGGCCACGCGCGGGCCCTGACCGAGTGCCAg                                                                                                     >  1:16659/1‑135 (MQ=255)                             cggggCATGGCCGGGGCGAGGCGAGGCCCCGCCGGGGCGAGGCCCTGCGGTAGGGCCCGCGCCTGGCGGCTGTCGGTGTAGGCCGGATCGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGt                                                                                                 <  2:400029/132‑1 (MQ=255)                             cGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGt                                                                                                 <  4:213264/135‑1 (MQ=255)                              gTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCGCGCGCGGGCCCTGGCCGCGTGCCACGCGtt                                                                                                >  1:203167/1‑135 (MQ=255)                                  cagggCCGGGGGGAGGCGATGGCCTGCCGGGTCGAGGCCCTGCTGGATGGCCCGGGCCAGGCTGCTGTCGGGGTAGGCCGGATAGCGGCGGCGGTTGGCCTGGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCa                                                                                            <  3:40518/132‑1 (MQ=255)                                  cATGGCCGGGTCTAGGCTATTGCCTGCCGGGGCTATGCCGTGCTGGATGGCACGCGCGAGGCTGATGTCGGGGTAGGCCGGATTGCGGCGGCCGTTTGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCa                                                                                            <  2:344802/135‑1 (MQ=255)                                              gAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑ggg                                                                              >  2:51307/1‑135 (MQ=255)                                                     ttGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCGGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGGCG‑AGGCC‑GGGGGGgcgc                                                                       <  7:148569/135‑1 (MQ=255)                                                         cAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGGCGCGCGCCGGCCCGGCCCCCGTGGCCGGCGGTGGCAGGGGCCGGCC‑ggggggggccccc                                                                    >  1:202755/1‑129 (MQ=255)                                                          aGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCGGCGGGCGGTGTAGGCCGGATAGCGGCGGCCGTGGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACAcc                                                                  <  7:265182/135‑1 (MQ=255)                                                                ggcGATGCCCTGCTGGAGGGCCCGCGCCAGGCGGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTGGGCCTCGCGCGGGCCCGGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTc                                                            <  3:227323/133‑1 (MQ=255)                                                                gTCGATGCCCGGCGGGAGGGCCCGCGCCAGGGGGCGGGCGGTGGAGGCCGGCTAGCGGCGGCCGTTGGCCGCGCGCGGGCCCGGACCGAGTGCCAGGCGTTGCCAGTCG‑CGGCC‑GGGGGGGCGCACACCGCCTTc                                                            <  7:265722/135‑1 (MQ=255)                                                                        ccTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCTGATAGCGGCGTCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTTCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGcaccgc                                                    >  3:391300/1‑135 (MQ=255)                                                                              ggAGGGCCCGCGCCAGGCTGCGGTCGGTGTAGGCCGGATCGCGGCGGCCGTTGGCCGCGCGCGGGCCCTGACCGAGTGCCAGGCGTGGCCAGGCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGcacggccacggc                                              <  8:214541/135‑1 (MQ=255)                                                                                aTGGCCCGCGCCAGGCGGCGGGCGGGGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCcg                                            <  4:296961/135‑1 (MQ=255)                                                                                 gggCCCGCGCCAGGCTGCTGGCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCGCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCcgt                                           <  3:453247/134‑1 (MQ=255)                                                                                  ggCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCGCGCGCGGGCCCGGACCGAGTGCCCGGCGGTGCCAGTCG‑AGGCC‑GGGGGGGCGCACCCCGCCTTCGGCACGGCCCCGGCcgtc                                          >  1:2922/1‑135 (MQ=255)                                                                                    cccGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCcgtcgt                                        >  3:149227/1‑135 (MQ=255)                                                                                       gcgcCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGc                                     >  4:61517/1‑135 (MQ=255)                                                                                            aGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCGGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCCCGGCCCCGGCCGTCGTTGCCAggg                                >  5:328854/1‑135 (MQ=255)                                                                                             ggcggcGGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCGCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGGCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGc                               <  6:220489/135‑1 (MQ=255)                                                                                                  ctgGCGGGGTAGGCGGGATCGGGGCGGCCGGTGGCCGCGCGCGGGCCCTGCCCGAGTGCCAGGCGTTGCCCGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTg                          <  6:75652/135‑1 (MQ=255)                                                                                                     tCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCCCGGCCGTCGTTGCCATGGCCCCTGgcg                       >  8:107737/1‑135 (MQ=255)                                                                                                      cGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCCGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACCCCCCCGTCGGCACCGCCCCGGCCGTCCTTGCAAAGGCAAAAGgcgc                      >  4:119371/1‑135 (MQ=255)                                                                                                              gcgggATCGGGGCGCCCGTGGGCCGCGCGGGGGCCCGGCCCGAGGGCCCGGCGTGGCCCGCCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGgcagc              <  3:29395/132‑1 (MQ=255)                                                                                                              gCCGGTTCGGGGCGGCCGTTGGCCCGGCGCGGGCCGTGACCGTGGGCCAGGCGTGGCCCGGCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGgcagc              <  4:132928/135‑1 (MQ=255)                                                                                                                    ccgcggcgGCCGTGGGCCGCGCGGGGCCCCGGACCGCGGGCCAGGCGTGCCCAGGCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTg        <  5:320184/133‑1 (MQ=255)                                                                                                                      gggggggcGGTTGCCCGCGCGGGGGCCCTGACCGTGGGCCCGGGGTGCCCAGTCG‑AGGCC‑GGGGGGGCGAACACGGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  8:269847/130‑1 (MQ=255)                                                                                                                      gcggcgGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCCGGGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCCACTGGCGCCGGGCCGCAcccggg       >  5:214074/1‑130 (MQ=255)                                                                                                                       gggcggcGGTGGGCCTGGCGGGGGCCCGGACCGAGGGCCCGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  6:200404/133‑1 (MQ=255)                                                                                                                       gggcggcCGTTGGCCTCGCGCGGGCCCTGACCGAGGGCCGGGGGGGGCCTGCCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  3:119371/133‑1 (MQ=255)                                                                                                                       cggcggcggTTGCCCGCGGGCGGGCCCTGCCCGATTGCCAGGGGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  3:462726/134‑1 (MQ=255)                                                                                                                       cggcggCCGTTGGCCGCGCGCGGGCCCGGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  7:114222/134‑1 (MQ=255)                                                                                                                       cggcggCCGTGGCCCGCGGGGGGCCCCGTCCCGAGGCCCAGGCGTGCCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGc      <  2:137369/134‑1 (MQ=255)                                                                                                                        ggcggcCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGCgg    <  7:163014/135‑1 (MQ=255)                                                                                                                          cggcCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGCgggg  <  8:361020/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                             GCGCTGGTCGGCCAAGGTCAGTTCGAAGCGTGGCATGGCCGGGTCGAGGCGATTGCCAGCCGGGTCGATGCCCTGCTGGATGGCCCGCGCCAGGCTGCTGTCGGTGTAGGCCGGATAGCGGCGGCCGTTGGCCTCGCGCGGGCCCTGACCGAGTGCCAGGCGTTGCCAGTCG‑AGGCC‑GGGGGGGCGCACACCGCCTTCGGCACGGCCACGGCCGTCGTTGCCATGGCAACTGGCGCAGGGCAGCACGCTGGCGGGG  >  ACB122/3944292‑3944547

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