Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,003,336	C→A	7.2%	W385C (TGG→TGT)	CJ019_01904 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,003,336	0	C	A	7.2%	164.3 / 4.6	56	W385C (TGG→TGT)	CJ019_01904	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base A (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

ATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCT  >  ACB122/2003205‑2003471                                                                                                                                     |                                                                                                                                        agtAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGCACAGGTCAATGTCCCaa                                                                                                                                      >  2:363006/3‑135 (MQ=255) aTTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCaaa                                                                                                                                     >  7:47902/1‑136 (MQ=255) aTTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCaa                                                                                                                                      >  4:57703/1‑135 (MQ=255) aTTACATCGTGGCCGTTAAAGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTTCGCCAGTAGAGAAAATGTCTGGGCCAATGAACAGGTCAATGTCCCaa                                                                                                                                      >  5:37967/1‑135 (MQ=255)      aTCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCg                                                                                                                                  <  2:69362/134‑1 (MQ=255)              ggTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAAcc                                                                                                                         >  4:104695/1‑135 (MQ=255)              ggTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCa                                                                                                                        <  8:211240/136‑1 (MQ=255)                    tGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCAtggc                                                                                                                    <  4:207058/134‑1 (MQ=255)                    tGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCAtggc                                                                                                                    <  3:18234/134‑1 (MQ=255)                         gACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTa                                                                                                              <  4:132422/135‑1 (MQ=255)                         gACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTa                                                                                                              <  4:78332/135‑1 (MQ=255)                            tGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGt                                                                                                           <  6:35722/135‑1 (MQ=255)                             gCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGtt                                                                                                          >  2:259431/1‑135 (MQ=255)                             gCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGtt                                                                                                          >  2:88404/1‑135 (MQ=255)                             gCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGCCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTg                                                                                                         >  8:238066/1‑136 (MQ=255)                                 tttGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGcc                                                                                                      >  5:306644/1‑135 (MQ=255)                                         cTTTGAGCAATGCCAAAGCCTCTTTGACAAGGGCGGCAGAG‑AAAATTGCGCCAGTAGGGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGcc                                                                                              <  1:355516/135‑1 (MQ=255)                                           ttGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCg                                                                                            >  1:81887/1‑135 (MQ=255)                                               gCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTcgc                                                                                        <  3:57703/135‑1 (MQ=255)                                                 aaatccAAAGCCTCTTTGTCAAATGCCGCAGAGAAAAATCGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTcgc                                                                                        <  8:142631/130‑1 (MQ=255)                                                 aaTGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCgg                                                                                      >  4:72763/1‑135 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcaca                                                                                 >  6:250276/1‑135 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcaca                                                                                 >  2:203569/1‑135 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  <  3:104925/134‑1 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  <  1:25108/134‑1 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  <  6:274525/134‑1 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAg                                                                                >  7:121678/1‑136 (MQ=255)                                                      cAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAg                                                                                >  8:29259/1‑136 (MQ=255)                                                           ccTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCgg                                                                            <  5:118513/135‑1 (MQ=255)                                                           ccTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCgg                                                                            <  2:219355/135‑1 (MQ=255)                                                           ccTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGc                                                                           >  7:142631/1‑136 (MQ=255)                                                                  gACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGATAAAATGTCTGCGCGAATGAAGAGGTCAATGTCCCAAACCGGCTTGAAGGATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAg                                                                      <  1:391480/134‑1 (MQ=255)                                                                        ggcggcAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAgc                                                               >  6:353237/1‑135 (MQ=255)                                                                            gcAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCaa                                                           >  6:165442/1‑135 (MQ=255)                                                                                              aGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCCACCGTCCTCCAGCAAAtt                                          >  5:229353/1‑135 (MQ=255)                                                                                                        aTGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATaa                                                                       >  7:159572/1‑96 (MQ=255)                                                                                                        aTGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATaa                                                                       <  8:159572/96‑1 (MQ=255)                                                                                                             tGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCACCCTTCCTCCAGCAAATTGCGAAGCTGCGCCAt                           >  5:200223/1‑135 (MQ=255)                                                                                                             tGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATg                          >  7:54425/1‑136 (MQ=255)                                                                                                             tGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATg                          >  7:60855/1‑136 (MQ=255)                                                                                                                gcCAATGAACAGGGCAATGTCCCAAACCGGCGTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGGCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATg                          <  1:357404/133‑1 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCg                                                                                            <  2:159290/65‑1 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCg                                                                                            >  1:159290/1‑65 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAg                                                                >  2:297982/1‑93 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAg                                                                <  1:297982/93‑1 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCACAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  >  4:192977/1‑75 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCACAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  >  4:186271/1‑75 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCACAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  <  3:192977/75‑1 (MQ=255)                                                                                                                  cAATGAACAGGTCAATGTCACAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTcac                                                                                  <  3:186271/75‑1 (MQ=255)                                                                                                                       aaCAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGgc                 >  1:151724/1‑135 (MQ=255)                                                                                                                        aCAGGTCACTGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGGCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGgc                 <  1:135773/134‑1 (MQ=255)                                                                                                                        aCAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTGGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGgc                 <  4:286837/134‑1 (MQ=255)                                                                                                                            gTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAg            <  2:81887/135‑1 (MQ=255)                                                                                                                               aaTGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGttt         >  2:338731/1‑135 (MQ=255)                                                                                                                               aaTGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGttt         >  1:188421/1‑135 (MQ=255)                                                                                                                                aTGTCCCAAACCGTCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGt       <  8:288558/136‑1 (MQ=255)                                                                                                                                     ccAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCt  <  7:37150/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAG‑AAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAATGAACAGGTCAATGTCCCAAACCGGCTTGAACCATGGCTGGCTACGTTGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCT  >  ACB122/2003205‑2003471

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