Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,003,043	G→A	7.1%	R5K (AGA→AAA)	CJ019_04696 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,003,043	0	G	A	7.1%	170.8 / 2.8	56	R5K (AGA→AAA)	CJ019_04696	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base A (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GGAGTCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTCGCTGT  >  ACB122/5002913‑5003175                                                                                                                                   |                                                                                                                                     ggAGTCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTCCTATGGCGAAAATGACATGGAAGCCTCCCCCCACGACAGCCTCAGCTTTCTCCCTGCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTc                                                                                                                                  >  3:125350/1‑135 (MQ=255)     tccaCCGCTCTGGGCAGCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGt                                                                                                                              <  8:454254/132‑1 (MQ=255)     tCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAATACTCCACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGt                                                                                                                              <  8:346176/135‑1 (MQ=255)     tCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGt                                                                                                                              <  7:255226/135‑1 (MQ=255)     tCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGt                                                                                                                              <  8:448914/135‑1 (MQ=255)     tCAACCGCGCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCTAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGt                                                                                                                              <  3:336372/135‑1 (MQ=255)      catccGCTCTGGGCATCAGTACGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGta                                                                                                                             <  8:192031/132‑1 (MQ=255)            ctctGGGCAGCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGa                                                                                                                       <  8:43649/135‑1 (MQ=255)                  gCATCAGTCCGTATTAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTGTCTCCCTTCTATCGATGAGGCCCAGCCATGAGCTACTCCAGATTCCTGTATACGGAGTAc                                                                                                                   <  7:240125/133‑1 (MQ=255)                  gCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTAc                                                                                                                   <  7:399603/133‑1 (MQ=255)                  gCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACg                                                                                                                  <  7:77798/134‑1 (MQ=255)                  gCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGa                                                                                                                 <  5:96667/135‑1 (MQ=255)                  gCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTCTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTAc                                                                                                                   <  3:231138/133‑1 (MQ=255)                  gCATCAGTCCGAATGACGGCCCCTGCTATGGCGAAAATTCCCTGGAACACTACACCCACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGa                                                                                                                 <  4:42430/135‑1 (MQ=255)                    aTCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGa                                                                                                                 >  3:127249/1‑133 (MQ=255)                    aTCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGa                                                                                                                 >  8:20313/1‑133 (MQ=255)                    aTCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGa                                                                                                                 >  6:241333/1‑133 (MQ=255)                        gTCCGTATGAGCTCCCCTGCTATGGCGAAACTGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGagcagc                                                                                                            <  6:142825/134‑1 (MQ=255)                         tCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGcc                                                                                                           >  8:328440/1‑134 (MQ=255)                         tCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAAATTCCTGTATACGGAGTACGAGCAGcc                                                                                                           >  5:422957/1‑134 (MQ=255)                         tCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAAATTCCTGTATACGGAGTACGAGCAGcc                                                                                                           >  5:429776/1‑134 (MQ=255)                                  gCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAg                                                                                                 <  7:449275/135‑1 (MQ=255)                                              tGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAggc                                                                                     >  1:461919/1‑135 (MQ=255)                                                 cGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAggcggc                                                                                  >  1:38761/1‑135 (MQ=255)                                                   aaaaTGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCt                                                                                 >  1:113774/1‑134 (MQ=255)                                                   aaaaTGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCt                                                                                 >  2:72119/1‑134 (MQ=255)                                                           cTGCAAGCCGACACCGACGACAGGCTGAGCTTTCTCGCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACTGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGAttt                                                                        <  5:354270/135‑1 (MQ=255)                                                            tGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGAtttt                                                                       <  3:276894/135‑1 (MQ=255)                                                                   cTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACt                                                                <  7:328440/135‑1 (MQ=255)                                                                   cTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAAATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACt                                                                <  6:422957/135‑1 (MQ=255)                                                                   cTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAAATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACt                                                                <  6:429776/135‑1 (MQ=255)                                                                        ccGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCa                                                           >  8:268589/1‑135 (MQ=255)                                                                              acAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCAtt                                                     >  4:123663/1‑135 (MQ=255)                                                                              acAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCAtt                                                     >  3:256262/1‑135 (MQ=255)                                                                                 gCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTggg                                                  <  7:287851/135‑1 (MQ=255)                                                                                         tttCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGAcaac                                           >  6:26776/1‑134 (MQ=255)                                                                                               ccTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGTCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGAcaacgcaaccg                                    >  5:34716/1‑133 (MQ=255)                                                                                                 ttCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGAc                                  <  8:95622/135‑1 (MQ=255)                                                                                                     aTCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTGATTGGGTGACAACGCAACGGACGa                                >  7:22464/1‑133 (MQ=255)                                                                                                         aTGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAgcg                          >  5:66678/1‑135 (MQ=255)                                                                                                           gAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGcgcg                        <  3:123663/135‑1 (MQ=255)                                                                                                           gAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGcgcg                        <  4:127249/135‑1 (MQ=255)                                                                                                             gggCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGt                       <  5:219875/134‑1 (MQ=255)                                                                                                                ccagccaTGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGc                   <  7:20313/135‑1 (MQ=255)                                                                                                                     caTGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCt              <  8:185335/135‑1 (MQ=255)                                                                                                                      aTGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCt              >  1:228841/1‑134 (MQ=255)                                                                                                                       tGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTa            >  6:237917/1‑135 (MQ=255)                                                                                                                           cTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACtt        >  6:399754/1‑135 (MQ=255)                                                                                                                           cTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACtt        >  6:451342/1‑135 (MQ=255)                                                                                                                            tACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTc       <  8:343711/135‑1 (MQ=255)                                                                                                                               tCCAGATTCGTGGATACGGAGGACGATCAGCCTCATGACGAGGTCGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCACGGGTCGAACAGCGCGTGGGCTTTCTTATCt         <  7:71509/130‑1 (MQ=255)                                                                                                                               tCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTc       <  2:461919/132‑1 (MQ=255)                                                                                                                               tCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGGGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACt         <  2:113774/130‑1 (MQ=255)                                                                                                                                 cAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTCGCTGt  >  7:445642/1‑135 (MQ=255)                                                                                                                                 cAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTCGCTGt  >  3:8802/1‑135 (MQ=255)                                                                                                                                 cAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGCGCTTTCTTAACTTCGCTg   >  5:91896/1‑134 (MQ=255)                                                                                                                                   |                                                                                                                                     GGAGTCAACCGCTCTGGGCATCAGTCCGTATGAGCTCCCCTGCTATGGCGAAAATGACCTGGAAGACTACACCGACGACAGCCTGAGCTTTCTCCCTTCTATCGATGAGGGCCAGCCATGATCTACTCCAGATTCCTGTATACGGAGTACGAGCAGCCTCATGACGAGGACGGTGAAGGCGGCTTTACGATTTTTAGCACTCCGCAGTCATTGGGTGACAACGCAACGGACGAACAGCGCGTGGGCTTTCTTAACTTCGCTGT  >  ACB122/5002913‑5003175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑