Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,249,801	C→T	7.1%	D251N (GAC→AAC)	CJ019_03985 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,249,801	0	C	T	7.1%	174.9 / 2.5	56	D251N (GAC→AAC)	CJ019_03985	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/20);  new base T (2/2);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AATACAGCTGCACGCTACTGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCAGGTAG  >  ACB122/4249670‑4249934                                                                                                                                    |                                                                                                                                      aaTACAGCTGCACGCTACTGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTcgcg                                                                                                                                    <  5:239545/135‑1 (MQ=255)       gCTGCACGCTACTGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGt                                                                                                                              <  2:254638/135‑1 (MQ=255)             cGCTACTGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTc                                                                                                                        <  2:55805/135‑1 (MQ=255)                   tGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTTAGTGAGCTCACTTTg                                                                                                                  >  5:216647/1‑135 (MQ=255)                        aGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCaa                                                                                                             >  3:35966/1‑135 (MQ=255)                           aCGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCCCTTTGCTCAAGcc                                                                                                          >  8:264836/1‑135 (MQ=255)                           aCGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGc                                                                                                           >  1:49996/1‑134 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCGATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  4:102303/1‑135 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  2:22195/1‑135 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  5:104951/1‑135 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  2:218417/1‑135 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  5:364021/1‑135 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAAc                                                                                                       >  8:71761/1‑136 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAAc                                                                                                       >  7:321754/1‑136 (MQ=255)                             ggTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGATGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCaa                                                                                                        >  5:316857/1‑135 (MQ=255)                                  aGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAAcgcgcg                                                                                                  >  7:174009/1‑136 (MQ=255)                                  aGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAAcgcg                                                                                                    >  4:84588/1‑134 (MQ=255)                                      caGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGAtc                                                                                               <  1:211247/135‑1 (MQ=255)                                       aGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTACGTGCCCTCACTTTGCTCCAGCCAACGCGCGAtct                                                                                              >  2:180194/1‑135 (MQ=255)                                       aGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCATGCCAACGCGCGAtct                                                                                              >  4:316060/1‑135 (MQ=255)                                       aGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGAtct                                                                                              >  6:158642/1‑135 (MQ=255)                                       aGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGAtct                                                                                              >  6:248117/1‑135 (MQ=255)                                            gAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTcgcccg                                                                                        >  8:310664/1‑136 (MQ=255)                                                aacaGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCt                                                                                     >  3:285659/1‑135 (MQ=255)                                                 acaGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTgg                                                                                   >  8:231945/1‑136 (MQ=255)                                                   aGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGt                                                                                  <  6:382645/135‑1 (MQ=255)                                                        gAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTTGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATcc                                                                            >  8:61123/1‑136 (MQ=255)                                                        gAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTTGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATc                                                                             >  2:376212/1‑135 (MQ=255)                                                          aCCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCg                                                                           >  5:154523/1‑135 (MQ=255)                                                          aCCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGc                                                                          >  8:246734/1‑136 (MQ=255)                                                                      ggcgCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTgcgc                                                               >  3:104709/1‑135 (MQ=255)                                                                      ggcgCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTgcgc                                                               <  3:172526/135‑1 (MQ=255)                                                                           agcagGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGa                                                          >  6:151034/1‑135 (MQ=255)                                                                               ggCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTc                                                      <  3:158893/135‑1 (MQ=255)                                                                               ggCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCa                                                     <  7:139967/136‑1 (MQ=255)                                                                                gCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCa                                                     <  1:314024/135‑1 (MQ=255)                                                                                gCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCa                                                     <  6:154523/135‑1 (MQ=255)                                                                                      ggTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTTGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCAtg                                               <  7:61123/135‑1 (MQ=255)                                                                                       gtgtGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTTGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCAtg                                               <  1:376212/134‑1 (MQ=255)                                                                                           gCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTg                                          <  4:324020/135‑1 (MQ=255)                                                                                             aGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGc                                        >  2:13508/1‑135 (MQ=255)                                                                                                  aaTCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAg                                                    >  8:94162/1‑118 (MQ=255)                                                                                                  aaTCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAg                                                    <  7:94162/118‑1 (MQ=255)                                                                                                     catcGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAggg                                 >  3:191066/1‑134 (MQ=255)                                                                                                     catcGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGt                                >  1:99946/1‑135 (MQ=255)                                                                                                     catcGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGGCGATTTCAGCCATGTGCTGGCTCAAGGGt                                <  2:188906/135‑1 (MQ=255)                                                                                                       tcGGGCATCTGCCGGCATGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTTCGCGTCGATTTAAGCAATGTGCTGGCTCAATGTTGa                              >  6:10928/1‑135 (MQ=255)                                                                                                       tcGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTg                               >  5:395400/1‑134 (MQ=255)                                                                                                        cGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGAc                             <  5:248117/135‑1 (MQ=255)                                                                                                                  ccGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGAGTCGATTTCAGCCATGTGCTGGCTCAAGGGGGACTGTCGCGATAc                  <  7:35352/136‑1 (MQ=255)                                                                                                                   cGGCAGGGTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATAc                  <  1:202152/135‑1 (MQ=255)                                                                                                                        ggTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCAtg                                               >  2:179365/1‑101 (MQ=255)                                                                                                                        ggTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCAtg                                               <  1:179365/101‑1 (MQ=255)                                                                                                                              tCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGc        <  4:35966/134‑1 (MQ=255)                                                                                                                              tCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCa       <  2:229962/135‑1 (MQ=255)                                                                                                                              tCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCa       <  1:13508/135‑1 (MQ=255)                                                                                                                               cATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCa       <  6:316857/134‑1 (MQ=255)                                                                                                                                    cgcgGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCAGGTAg  <  3:81584/134‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      AATACAGCTGCACGCTACTGGCCAGTACGGTAAAGCACAGGCTGAACAACAGGATGAACCGCAGCAGGCGGCGCAGCAGGCCATCGGTGTGCAGGCGAATCATCGGGCATCTGCCGGCAGGTTGATCATGTCGCGGTAAGTGACCTCACTTTGCTCAAGCCAACGCGCGATCTCGCCCGCCTGGTTCATCCGCTGCAACTGCGCGTCGATTTCAGCCATGTGCTGGCTCAAGGGTGACTGTCGCGATACCGCCACCCGCAGGTAG  >  ACB122/4249670‑4249934

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