Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,986,622	T→G	29.8%	D31A (GAT→GCT)	CJ019_03763 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,986,622	0	T	G	29.8%	40.8 / 2.8	47	D31A (GAT→GCT)	CJ019_03763	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/13);  new base G (4/10);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  ACB122/3986496‑3986752                                                                                                                               |                                                                                                                                   gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCgg                                                                                                                           >  7:121934/1‑136 (MQ=255)  cccccATCCGCCCCTAACCCACAAACTAATCCGGGGACACGCCAAGGGGGGCCAGCAGCTGGTCCCGATGCGCCGGGGAGAACACCGGGTCCTTGGCTACGCGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCgg                                                                                                                           <  5:296340/132‑1 (MQ=255)      cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGcgc                                                                                                                       >  3:286442/1‑133 (MQ=255)       agcaccACCTGAGCCAGCAACTGATCAGGGTACAAGGCATTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGCGAACACGGGGTCCTTGGCTACGTGCGTGGCAGCGTAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                      <  6:99100/130‑1 (MQ=255)                 tacccAGCAAATGCTCAGGTTACAACGCACTTGTGGCCTGCATCGGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGACCGGGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                           <  8:222146/133‑1 (MQ=255)                   gCCTGCATCTGATCCGGGTACAAGGCAATGGGGCCCAGCATCTGGTCCAGTGGCTCCTGGGCGACCACGGGGTCCTTGGCTACGAGCGTGGCAGCGTAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGc                                                                                                          <  2:103574/135‑1 (MQ=255)                   gCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGACCACCGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGc                                                                                                          <  2:109044/135‑1 (MQ=255)                     cAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGATCCGGCGCGGGCACGTCTGCGGCCGGCAcc                                                                                                       >  8:95679/1‑136 (MQ=255)                         aaCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATGCGGCGCGGGCACGTCAGCGGCCGGCCcccgg                                                                                                    >  6:70879/1‑135 (MQ=255)                         aaCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGCAGTACAACGGGTCCTTGGCTACTATCGTGGCAGCGGAAGTTTACGGATCCGGCGCGTGCACGTCTGCGGCCTGCAccagg                                                                                                    >  2:365801/1‑135 (MQ=255)                          aCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCCccaggcc                                                                                                  >  8:191638/1‑136 (MQ=255)                          aCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCAccaggcc                                                                                                  >  8:282624/1‑136 (MQ=255)                                  gggTACACGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAg                                                                                          <  8:295199/136‑1 (MQ=255)                                  gggTACAAGGCTATGGGGGCCATCATCTGGTGCAGTTGCGCCTGGGAGACCACGGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAg                                                                                          <  8:367508/136‑1 (MQ=255)                                   ggTACAAGGCTCGGGGGGCCAGCAGCTGGGCCAGGGGCTCCGGGGAGCACCCGGGGGCGGTGCGTACGCGGGTGGCCGCGGACGGCGCGGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAg                                                                                          <  6:271927/135‑1 (MQ=255)                                   ggTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCCGCGCGGGCCCGTCTGCGGCCAGCACCACGCCCCGGGCAg                                                                                          >  4:302678/1‑135 (MQ=255)                                      aCAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGGGGCGCCGGCGCGGGGCCGGCGGCGGCCGGCCCCCGGCCCCGGGCAgcgcg                                                                                      >  7:252224/1‑135 (MQ=255)                                         aGGCAAGGGGGGCCAGCCTCTGGTCCAGTGGCCCCTGGGAGACCCCCGGGTCCTGGGCTACGCGCGGGGCCGCGGACGGCGCCGGCCCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAg                                                                                   <  8:252224/136‑1 (MQ=255)                                                gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGCGCCGGCGCGGGCACGGCGGCGGCCGGCCCCAGGCCAGGGGCCGCGCGCAGGCAAAGc                                                                            >  8:339759/1‑136 (MQ=255)                                                gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCCGATCCGGCGCGGGCACGTCATCGGCCGGCACCCCGCCCGGGGCCGCGCGCAGGCaaa                                                                              >  4:300421/1‑134 (MQ=255)                                                gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCCGCGCGGGGCCGGCTGCCGCCAGCCCCACCCCAGCGGCCGCCCGCCGGCCAAg                                                                             >  8:29148/1‑135 (MQ=255)                                                gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCGCAGGCAAAg                                                                             >  4:248293/1‑135 (MQ=255)                                                  gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGGCCGTCGGCGGCCGGCCCCCGGCCAGGGGCGGCGCGGAGGcccacca                                                                           >  5:236531/1‑129 (MQ=255)                                                       aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCGGCGCGGGCCCGGGTGCGGCCTGCCCCCGGCCCGGGGGCGCCCGCAGGCCAACCAGTggc                                                                      >  1:272892/1‑133 (MQ=255)                                                        gCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATACGGCGCGGGCACGTCTTCGGCCAGCACCACGCCCGGGGCAGCGCTCAGGCAAAGCAGTGCCg                                                                     >  4:70915/1‑135 (MQ=255)                                                          aTCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCCGCGCGGGCCCGGCGGCGGCCTGCCCCCGGCCCGGGGGCGCGCTCAGGCCAAGCAGTGGCg                                                                     >  1:70552/1‑133 (MQ=255)                                                          aTCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTACGGCCCGCGCACGTCTGCCGCCCGCACCAGGCCAGGGGCCGCCCTCAGGCAAAGCAGTGCCgcc                                                                   >  2:149259/1‑133 (MQ=255)                                                                   aggtGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCa                                                          <  2:293610/132‑1 (MQ=255)                                                                      tGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGGCCGGCTGCGGCCGGCCCCACGCCCGGGGGCGCCCGCCGGCAAAGCCGGGCCGGCCCgcaagcaac                                                        >  2:236830/1‑133 (MQ=255)                                                                            tGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATg                                                   <  6:313137/133‑1 (MQ=255)                                                                               gagaACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  7:282624/136‑1 (MQ=255)                                                                                agaACCCCGGGTCCTTGGCTACGAGGGGGGCTGCGGACGGCGCCGGAGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  3:302678/135‑1 (MQ=255)                                                                                    cacGGGGTCGTTGGCTAGGCGCGTGGGAGGGGACGGCGCCGGTCGGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCAt                                         <  6:370359/135‑1 (MQ=255)                                                                                        ggggCCTGGGCTACGCGGGGGGGCGGGGACGGCGCGGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCt                                     <  1:278045/135‑1 (MQ=255)                                                                                          gTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCACGTCTGCGGCCTGCCCCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTa                                   >  8:89572/1‑135 (MQ=255)                                                                                             cTTGGCACCGCGCGGGGCCGCGGACGGCGCCGGAGCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATcc                                <  5:70879/135‑1 (MQ=255)                                                                                                     cGAGCGTGGCAGCGGATGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAgg                        <  2:70552/135‑1 (MQ=255)                                                                                                       aGCGTGGCAGCGGAAGGCGCCGGCTCCCGCGCGGGCCCGGCGGCGGCCGGCCCCAGGCCCGGGGCCGCGCGCACGCAAACCAGTGCCGACACGAAAGCAAACAAATGCATGCGCCTCGCGTTTCCCCAAGAgg                        >  4:178014/1‑133 (MQ=255)                                                                                                           gggCCGCGGACGGGGCCGGATCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  7:191638/135‑1 (MQ=255)                                                                                                            ggCAGCGGAATGTGCCGGATCCGGCGAGGGAACGTCTGCGGCCTGCACCAGGCCAGGGGCAGGGCTCAGGCAAAGCTGTGCCGACAGGAAATCAAACAATGGCATGCGCATCGCTTATCCCCTTGAGGTCAGAcc                 <  1:365801/135‑1 (MQ=255)                                                                                                            ggCAGCGGAAGGCGCCGGCTCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  3:371964/135‑1 (MQ=255)                                                                                                            ggCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  4:286442/135‑1 (MQ=255)                                                                                                            ggCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  2:272892/135‑1 (MQ=255)                                                                                                                  ggACGGCGCCGGTTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGGAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg           <  1:149259/135‑1 (MQ=255)                                                                                                                      ggCGCCGGCTCCGGCGCGGGCACGTCTGCCGCCTGCCCCACGCCAGGGGCAGAGCTCAGGCAAAGCCGGGCCGCCAGGAAAGCAAACAATGGCATGCGCAGCGCATATACCCATGAGGTAAGCCcgcggtggt         >  6:230607/1‑128 (MQ=255)                                                                                                                           cGGCTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGt  <  3:248293/135‑1 (MQ=255)                                                                                                                           cGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGt  <  4:375621/135‑1 (MQ=255)                                                                                                                               |                                                                                                                                   GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  ACB122/3986496‑3986752

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