Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,930,927	A→G	9.8%	F396L (TTC→CTC)	CJ019_01847 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,930,927	0	A	G	9.8%	111.3 / 5.0	41	F396L (TTC→CTC)	CJ019_01847	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base G (2/2);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

CGAAGAGGAGTTTGTTAACCACGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCAACTGCAGG  >  ACB122/1930800‑1931045                                                                                                                                |                                                                                                                       cGAAGAGGAGTTTGTTAACCACGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATg                                                                                                                 >  5:288470/1‑135 (MQ=255)   aaGAGGAGTTTGTTAACCACGGAAGCACTCCAGCATGGCGTCAGTGTTGACGCGGGAGTTGAGGTTGCCGCAATGGCCGCCAAGGGGGTAGAGGGTCCGCCGTTCGCCGAACACCTTCCCCAGGAAGCCCAAGTc                                                                                                               >  6:36878/1‑135 (MQ=255)      aggagTTTGTTAACCACGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGcc                                                                                                            >  8:21081/1‑135 (MQ=255)         agTTTGTTAACCACGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGGTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCgg                                                                                                         <  3:260007/135‑1 (MQ=255)                      cGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGAcg                                                                                            >  4:110538/1‑135 (MQ=255)                        gaagaaCTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAGGCCAATGTCGCCCGGGCCGAGAATGAtgtc                                                                                          >  1:213474/1‑135 (MQ=255)                        gaagaaCTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAGGCCAATGTCGCCCGGGCCGAGAATGAtgtc                                                                                          >  1:217561/1‑135 (MQ=255)                         aagaaCTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGAcgtcg                                                                                         <  6:288470/135‑1 (MQ=255)                                    cATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCg                                                                                                          >  6:147566/1‑107 (MQ=255)                                    cATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCg                                                                                                          <  5:147566/107‑1 (MQ=255)                                    cATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTtgt                                                                               >  5:270897/1‑134 (MQ=255)                                      tGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCa                                                                            <  6:260458/135‑1 (MQ=255)                                           tCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATCGCCGCCATGGGGGTAGAGTGGCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGGCGTCGGCGTTGTGCATCACg                                                                       <  5:96950/135‑1 (MQ=255)                                              cTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGc                                                                     <  8:110376/134‑1 (MQ=255)                                                  tGACGCGGTAGTTGAGGGTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGGCGTCGGCGTTGTGCATCACGGCGATCt                                                                <  7:50694/135‑1 (MQ=255)                                                          tAGTTGAGGTTGCCGCAATGGCCGCCATGGTGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCActg                                                        <  1:175731/135‑1 (MQ=255)                                                            gTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCActgct                                                      >  4:217531/1‑135 (MQ=255)                                                              tGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACtgctgc                                                    >  4:276396/1‑135 (MQ=255)                                                                   ttGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAgg                                               <  3:276396/135‑1 (MQ=255)                                                                   ttGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAg                                                <  5:74052/134‑1 (MQ=255)                                                                         cAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAGGCCAATGTCGCCCGGGCCGAGAATGATGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCt                                         <  2:213474/135‑1 (MQ=255)                                                                         cAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAGGCCAATGTCGCCCGGGCCGAGAATGATGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCt                                         <  2:217561/135‑1 (MQ=255)                                                                           aTGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTc                                       <  8:231147/135‑1 (MQ=255)                                                                           aTGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTc                                       <  7:114849/135‑1 (MQ=255)                                                                           aTGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTc                                       <  5:36878/135‑1 (MQ=255)                                                                                      gggggTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACa                            <  1:67566/135‑1 (MQ=255)                                                                                           tAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCCCGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTg                       >  1:65750/1‑135 (MQ=255)                                                                                                ggTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGATAATGACGTCGTCGGCGTTGTTCAACACGGAGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGCGCATACAGGCAGACCTg                  >  6:283351/1‑135 (MQ=255)                                                                                                ggTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACtgctgc                                                    <  4:278287/101‑1 (MQ=255)                                                                                                ggTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACtgctgc                                                    >  3:278287/1‑101 (MQ=255)                                                                                                 ggcAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGACGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTgg                 <  7:263734/133‑1 (MQ=255)                                                                                                     gCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGa             >  4:81177/1‑135 (MQ=255)                                                                                                     gCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGa             >  1:257390/1‑135 (MQ=255)                                                                                                      ccGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCCGGGCATACAGGCTGACCTGGTTGAc            >  2:61282/1‑135 (MQ=255)                                                                                                      ccGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGa             >  3:62870/1‑134 (MQ=255)                                                                                                      ccGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGAc            <  1:142748/135‑1 (MQ=255)                                                                                                       cGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCCCGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTggc            >  5:10221/1‑132 (MQ=255)                                                                                                         ttCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCaa         <  6:270897/135‑1 (MQ=255)                                                                                                              cGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCAACTGCa    >  8:50386/1‑135 (MQ=255)                                                                                                              cGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGGCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCAACTGCa    <  7:280218/135‑1 (MQ=255)                                                                                                                aaCACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCAACTGCAgg  <  8:95389/135‑1 (MQ=255)                                                                                                                                |                                                                                                                       CGAAGAGGAGTTTGTTAACCACGGAAGAACTCCAGCATGGCGTCACTGTTGACGCGGTAGTTGAGGTTGCCGCAATGGCCGCCATGGGGGTAGAGGGTCAGCCGTTCGCCGAACACCTTGCGCAGGAAGCCAATGTCGCCCGGGCCGAGAATGACGTCGTCGGCGTTGTGCATCACGGCGATCTTGTCACTGCTGCGCAGGTAGTCTTCCAGGGCATACAGGCTGACCTGGTTGACCAACTGCAGG  >  ACB122/1930800‑1931045

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