Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,986,621	A→C	33.2%	D31E (GAT→GAG) ‡	CJ019_03763 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,986,621	0	A	C	33.2%	55.5 / 3.8	54	D31E (GAT→GAG) ‡	CJ019_03763	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/24);  new base C (12/4);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGC  >  ACB122/3986496‑3986755                                                                                                                              |                                                                                                                                       gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGACCACCGCGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGcgcg                                                                                                                               >  4:294987/1‑135 (MQ=255) gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGGCGGGcgcg                                                                                                                               >  1:364992/1‑135 (MQ=255)      cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCCGCGGCGCCCGCGCGGGcgc                                                                                                                          >  1:54185/1‑133 (MQ=255)       acaaGCACCTGAGCCACCAACTGAGCAGGGAACAAGGCAAGGGGGCCCAGCATCGGGTCCAGTTGCTCCGGGGAGAACCCCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGACCCGGCGCGGGCACg                                                                                                                         <  2:364992/132‑1 (MQ=255)       aTCAGCCCCTGAGCCAGCAACTGATCAGGGTACAAGCCAATGGGGGCCAGCAGCTGGTCCAGTTCCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                         <  5:232957/135‑1 (MQ=255)       aTCAGCACCTGAGCCAGCAACTTATCTGGGTACAAGGCACTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGTACAACCCCGGGTCGTGGGCTACGAGGGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                         <  8:151392/135‑1 (MQ=255)            cACCTGAGCCCGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCCTCGGAAGGCGCCGGATCCGGCGCGCGCACGTCTGc                                                                                                                    >  6:386259/1‑135 (MQ=255)            cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTCCGAGCGTGGCAGCGGAAGGCGCCCGATCCGGCGCGGGCACGTCTGc                                                                                                                    >  4:278369/1‑135 (MQ=255)            cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCGGGGCCAGTTGCGCCTGGGACAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGc                                                                                                                    <  7:34767/135‑1 (MQ=255)            cACCTGAGCCAGCAACGGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGc                                                                                                                    <  6:307937/135‑1 (MQ=255)                  cgCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                              <  1:268688/134‑1 (MQ=255)                    ccAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCCCGTCTGCGGCCGGCa                                                                                                            >  2:286966/1‑135 (MQ=255)                                 aGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGCTCCGGCGCGGGGGCGGCGGCGGCCGGCCCCACGCCCCGGGc                                                                                               >  7:241471/1‑135 (MQ=255)                                     tACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGTCGGCGGCCGGCACCAGGCCCGGCGCcgcg                                                                                           >  3:125022/1‑135 (MQ=255)                                      aCAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCCGCGCGGGGCCGGCGGGGGGCAGCGCCAGCCCCGGGGCCgcgc                                                                                          >  3:382303/1‑135 (MQ=255)                                           gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGGGCGGGCCCGTCTGCGGCCGGCCCCCGGCCCGGGGCCgcgcggcg                                                                                      >  6:310644/1‑130 (MQ=255)                                           gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCACCTCTGCGGCCTGCACCAGGCCAGGGGCCGCGCGCAg                                                                                      >  8:458507/1‑134 (MQ=255)                                             aaTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCGGCGGCCTGCCCCCCGCCCGGGGCAGCGCGCAGGc                                                                                    >  2:119438/1‑134 (MQ=255)                                                gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGATCCGGCGCGGGCACCGCGGCGGCCTGCACCAGGGCAGGGGCAGCGCGCAGGCaaa                                                                                 >  6:244376/1‑134 (MQ=255)                                                  gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGCACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCCCCGGATCCGGAGCGGGCACGTCTGCGGCCGGCGACACGCCACCGGCAGCGCGCAGGCAAAgaa                                                                              >  8:431991/1‑133 (MQ=255)                                                  gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACCCCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGCTCCGGCGCGGGCCCGGCTGCGGCCGGCCCCAGGCCCGGGGCAGCCCGCAGGCACAGCa                                                                              >  2:146875/1‑135 (MQ=255)                                                  gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCGGCGGGGGCACGTCTGCGGCCTGCACCAGGCCCGGGGCAGCGCTCAGGCCAAGCa                                                                              >  3:66998/1‑135 (MQ=255)                                                  gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCa                                                                              >  4:85419/1‑135 (MQ=255)                                                       aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCGGCGCGGGGCCGGCAGCGGCCCCCCCCACCCCCCGGGCCGCGCGGAGGCCAAGCAGTGcc                                                                         >  7:264747/1‑135 (MQ=255)                                                       aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTGCGGCGGGGGCCCCTCTGCGGCCGGGGCCCCGGCCGGGGCCGCGCGCCGGCCACGCAGGGcc                                                                         >  5:438779/1‑135 (MQ=255)                                                       aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGGCTGCGGCCGGCCCCACGCCCCGGGCCGCGACCCGGCCAAGCAGTGcc                                                                         >  4:198039/1‑135 (MQ=255)                                                        gCATCGGGTCCAGTTGCGCCTGGGACAACACCGGGGCCTGGGCTCCGAGCGGGGCCGCGGACGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCg                                                                        <  6:413185/135‑1 (MQ=255)                                                           ccTGGCCCCGTTGCCCCGGGGACAACACGGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACa                                                                     <  5:76001/134‑1 (MQ=255)                                                                      tGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCaa                                                            <  6:5868/133‑1 (MQ=255)                                                                       gCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCCCGGGCCCCTCTGCGGCCAGCACCCGGCCACGGGCAGCGCGCACGCAAAGCAGAGCCGCCAGGAAAGCAAAc                                                          >  2:299040/1‑134 (MQ=255)                                                                       gCTCCTGGGAGAACACCGGGTCCGTGGCTAAGAGCGTGGCAGCGGACGGAGCCGGCTCCGGCGCCGGCACCGCATCGGCCATCCCCACCCCCTGGGCCTCGCTCATGCAAAGACCTGCCCACAGCAAAGCAAAc                                                          >  6:27956/1‑134 (MQ=255)                                                                            tggcacaaCCCCGGGTCCTTGGCACCGCGGGTGCCCGCTGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGc                                                    <  3:71109/129‑1 (MQ=255)                                                                              ggAGAACACGGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                                  <  2:146334/135‑1 (MQ=255)                                                                              ggAGAACACGGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                                  <  8:354674/135‑1 (MQ=255)                                                                                acaaCCCCGGGGCCTTGGCTACGAGCGGGGCAGCGGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                                <  1:41144/133‑1 (MQ=255)                                                                                       ggggTCCTTGGCTCCGGGCGTGGCCGCGGACGGCGCCGGCGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGc                                         <  5:98741/134‑1 (MQ=255)                                                                                        gggTCATTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCt                                        <  6:319582/135‑1 (MQ=255)                                                                                         tgTCCTTGGCTACGAGCGGTGCCGGGGAAGGGGCCGGAGCCGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                       <  3:267111/134‑1 (MQ=255)                                                                                         gggcattGGCTACGAGGGGGGCTGCGGACGGCGCGGTATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                       <  6:53627/130‑1 (MQ=255)                                                                                           tCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACCTCTGCGGCCTGCACCAGGCCAGGGGCCGCGCGCAGGCAAAGCACTGCCGCCACGAAAGCAAACCATGGCAAGCGCATCGCGTAt                                     >  4:66052/1‑135 (MQ=255)                                                                                                            ggccggggACGGCGGCGGACCGGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  2:2172/129‑1 (MQ=255)                                                                                                            ggCTGCGGCCGGCGCCGCCCCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  7:267449/135‑1 (MQ=255)                                                                                                            ggCCGCGGCCGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  7:51382/135‑1 (MQ=255)                                                                                                            ggCAGGTGACGGCGCCGGCCCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  3:294987/135‑1 (MQ=255)                                                                                                            ggCAGGGGCCGGCGCCGGCGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  1:182898/135‑1 (MQ=255)                                                                                                            ggCAGCGGAAGGCGCCGGTGCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  6:343006/135‑1 (MQ=255)                                                                                                            gcccggggACGGCTCCGGAGGCGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGGCAGAcc                    <  7:24644/129‑1 (MQ=255)                                                                                                              cAGCGGAAGGCGCCGGGTCCCGCGCGGGGCCGTCTGCGGCCTGCACCCGGCCACGGGCCGCGCTCCCGGAAACCAGGGCCGCCAGGCAATCAAACAATGGCATGCGCCACGCGTATACCCATGAGGTGAGAcccc                  >  8:267449/1‑133 (MQ=255)                                                                                                              cAGCGGAAGGCGCCGGGTCCCGCGCGGGCCCGTCTGCGGCCGGCACCCGGCCCGGGGGCGCGCGCAGGGCAACCAGTGCCGACACGAAAGCAAACAATGGCATGCGCATCGCGGATACCCCAGAGGTCAGACcgc                  >  2:19804/1‑135 (MQ=255)                                                                                                                  gtaaGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACAAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  8:202848/133‑1 (MQ=255)                                                                                                                  ggAGGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  6:418273/135‑1 (MQ=255)                                                                                                                  ggACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  8:241471/135‑1 (MQ=255)                                                                                                                      ggcggcGTAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTa          <  3:66052/135‑1 (MQ=255)                                                                                                                        cgcCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGTAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAg        <  4:409171/135‑1 (MQ=255)                                                                                                                              aTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGc  <  4:45608/135‑1 (MQ=255)                                                                                                                              |                                                                                                                                       GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGC  >  ACB122/3986496‑3986755

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