Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,718,476	T→G	27.1%	T61P (ACC→CCC)	CJ019_01614 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,718,476	0	T	G	27.1%	50.8 / 5.0	47	T61P (ACC→CCC)	CJ019_01614	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/19);  new base G (11/3);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 2.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGC  >  ACB122/1718343‑1718610                                                                                                                                      |                                                                                                                                       aGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTc                                                                                                                                       >  2:22712/1‑135 (MQ=255)   cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGCACg                                                                                                                                    >  7:147271/1‑136 (MQ=255)   cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGCGGCAc                                                                                                                                     >  1:318268/1‑135 (MQ=255)      cTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCAccg                                                                                                                                   >  8:134978/1‑132 (MQ=255)       ttGCCATTCAGCTCAATGCGCGCAGTTCCCGCATAGCGGTGGCCAGGCTTGACCGGCGCGGGGTCATAGGTGAGATGGAACGCCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCa                                                                                                                                 <  2:267360/135‑1 (MQ=255)                aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTcc                                                                                                                        <  4:299944/135‑1 (MQ=255)                aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCAGAGGTGAGATGGAACGGCAGCGGCCCGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCCCGGCAGGGGCGTc                                                                                                                         >  3:291567/1‑134 (MQ=255)                       tGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCGGGGCGCGCCAGGGGCCCGGCAGGGGCGGCCCGCAGCgc                                                                                                                >  3:283083/1‑135 (MQ=255)                        gcgcgcACTACCCGCATAGCGGGGGCCAGGCTTGCCGGGCGGGGGGGCATAGGTGAGCGGGACCGGCCGCGGCCCGTGGCCTTTGCCCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                                               <  3:354001/136‑1 (MQ=255)                              aCTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                                          <  1:47899/135‑1 (MQ=255)                               cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGTAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCTGGCGCCCCCGGCTCACGGGCCGGGCGTGCAGCACCGCCACAt                                                                                                           >  6:168144/1‑133 (MQ=255)                                   ccGCATAGCGGTGGCCAGGCTTGACCGTCGCGGGGGCATAGGTGAGAGGGAACGGCAGCGGCACGGTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAg                                                                                                    <  8:276440/136‑1 (MQ=255)                                    cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGCCCGGCCGGGGGGTCCATGCGCGCCACAACCAGCacg                                                                                                   >  3:198414/1‑134 (MQ=255)                                         aGCGGTGGCCAGGCTTGACCGGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                                              <  7:166329/136‑1 (MQ=255)                                           cGGTGGCCAGGCTTGACCTGCGCGGGGTCTTAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGGGCGGCACGGGGACGGCCGGGGCGGCCCTCCGCGCCCCAACCTGCAGGTGCACGc                                                                                            >  3:29247/1‑136 (MQ=255)                                           cGGTGGCCAGGCTTGACCTGCGCGGGGTCTTAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCGCGGCAGGGGGGGCCCGCCGCGCCACCTCCGGGAGTTCCCCg                                                                                             >  6:10900/1‑135 (MQ=255)                                             gTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCCCGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCt                                                                                           >  8:192483/1‑135 (MQ=255)                                                 ccAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGCCCGGGGCCCGGCCGGGGCGGCGCGCCGCGCCCCCTCCGGCAGGTGCACGCGCAggg                                                                                      >  7:10010/1‑136 (MQ=255)                                                 ccAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCACGGCCGGGGCGTCCATCAGCGCCCCCTCCGGCAGTTCCACGCTCAgg                                                                                       >  2:214374/1‑135 (MQ=255)                                                   aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTGGCGGGCGGGGCGTCCCTGAGCGACCGCTCCGGCCGGTCCCCGCGCAGGGt                                                                                     >  2:22363/1‑135 (MQ=255)                                                       ttGACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                                 <  2:37897/135‑1 (MQ=255)                                                        tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggc                                                                               >  8:329839/1‑136 (MQ=255)                                                        tGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTGCCGGCCGGGGCGTCCATAAGCGCCCCCACCAGCCGGTCCCCGCTCAGGGTggcgg                                                                                >  2:90937/1‑135 (MQ=255)                                                         gACCTGCGCGGGGTCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggc                                                                               <  8:238463/135‑1 (MQ=255)                                                           ccTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                                            <  3:320569/136‑1 (MQ=255)                                                           ccTGCGCGGGGGCATAGGGGAGAGGGAACGGCCCCGGCACGTTGCCTTGGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                                            <  3:147589/136‑1 (MQ=255)                                                             agcgcGGGGGCATCGGGGAGATGGACCGGCAGCGGCACGTGGCCTTTGACCGGCCCGGCCGGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcaggcg                                                                          <  7:320097/135‑1 (MQ=255)                                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGCGCGGCCGGGGCCCGGGCGGGGCGGGCCTCCGCGCCACCTCCTGCAGGTCCCCGCGCAGGGTggcggccggggggcgg                                                                     >  7:227728/1‑132 (MQ=255)                                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGCACGGGCGGGGCGGCGCTCAGCGCCACATCCCGCCGTTGCACGCTCCGGGgggcggcgggcggccg                                                                      >  2:187165/1‑133 (MQ=255)                                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGGGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGcggcaggcggcagg                                                                     >  4:95219/1‑136 (MQ=255)                                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGcggcaggcggcagg                                                                     >  8:181382/1‑136 (MQ=255)                                                                  ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGGGCATCAGCGCCACATCCTGGACTTCCACCCGCACGGTGgcggcaggcggcag                                                                      >  2:299937/1‑135 (MQ=255)                                                                   gggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGCGTCACGGCAGGGGCGTCCAACAGCCCCCCATCCTGCACATCCACGCACAGCGTGGcggcacgcggcagg                                                                     >  5:83414/1‑135 (MQ=255)                                                                         gaGGTGAGAGGGAACGGCAGCGGCACGTTGCCTGTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAAt                                                                <  2:88809/133‑1 (MQ=255)                                                                          aGGTGAGATGGACCGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgc                                                              <  4:29247/135‑1 (MQ=255)                                                                                 aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                                                 >  4:125310/1‑77 (MQ=255)                                                                                 aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                                                 <  3:125310/77‑1 (MQ=255)                                                                                 aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAg                                                       <  5:383829/135‑1 (MQ=255)                                                                                   ggAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCCCGGCCGGGGCGGCCATCAGCGCCACCACCGGCAGGTCCACGCGCCCGGGGGCGGGCGGCGGCCGGGGCATGCGCTGCAg                                                       >  1:309880/1‑133 (MQ=255)                                                                                          cAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                                             <  4:283083/136‑1 (MQ=255)                                                                                                        cTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCt                               <  7:181382/136‑1 (MQ=255)                                                                                                                    cggccGGGCGCGCCCGGGGCACGGCAGGGGCGGCCCTCAGCGCCCCATCCTGCAGTTCCACGCGCAGGGTGGCGGCCGGCGGCCGGGCCATGCGCGGCCGGGAGAAGACTTCACCCGCCAGGCGGGGCGGGGtgg                    >  2:287711/1‑135 (MQ=255)                                                                                                                       ccTGGCGCGCCTGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtgg                <  3:218534/136‑1 (MQ=255)                                                                                                                        cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtgg                <  2:159892/135‑1 (MQ=255)                                                                                                                            ggcgcCAGTGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGtt            <  1:115487/134‑1 (MQ=255)                                                                                                                                   gggTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTgg     <  2:340259/135‑1 (MQ=255)                                                                                                                                     gTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGc  <  4:176067/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       AGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGC  >  ACB122/1718343‑1718610

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