Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,917,312	C→T	7.1%	E283E (GAG→GAA)	CJ019_04613 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,917,312	0	C	T	7.1%	179.2 / 3.7	56	E283E (GAG→GAA)	CJ019_04613	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/32);  new base T (2/2);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCAGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCAGCGAGCGGCTGACCTG  >  ACB122/4917180‑4917443                                                                                                                                     |                                                                                                                                    gCAGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTTACAGCACCGCTGGCGATCAGCTc                                                                                                                                   >  8:357507/1‑135 (MQ=255)   aGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTc                                                                                                                                   <  8:85915/133‑1 (MQ=255)   aGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCgg                                                                                                                                 >  6:389793/1‑135 (MQ=255)   aGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCgg                                                                                                                                 >  4:319673/1‑135 (MQ=255)        cggatcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGt                                                                                                                            >  3:136901/1‑135 (MQ=255)          gatcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTg                                                                                                                           <  1:233017/134‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  1:12565/135‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  6:143213/135‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  5:423446/135‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  4:246982/135‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  1:345405/135‑1 (MQ=255)           atcTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGGGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcg                                                                                                                         <  8:305198/135‑1 (MQ=255)             cTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCCTCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcggc                                                                                                                       >  2:133353/1‑135 (MQ=255)             cTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGTTCGGCCAGTgcgg                                                                                                                        >  7:364933/1‑134 (MQ=255)             cTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGTTCGGCCAGTgcgg                                                                                                                        >  7:371226/1‑134 (MQ=255)             cTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcggc                                                                                                                       >  3:287318/1‑135 (MQ=255)             cTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTgcgg                                                                                                                        >  4:329499/1‑134 (MQ=255)                     ccTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGa                                                                                                               <  3:322353/135‑1 (MQ=255)                           gCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGt                                                                                                          <  7:458684/134‑1 (MQ=255)                           gCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCAGCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGt                                                                                                          <  2:102033/134‑1 (MQ=255)                               acgacTGTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGca                                                                                                     <  8:457334/135‑1 (MQ=255)                               aagacTTTTTCATTCCCCTACCCTCAAAATCAGCCGTTGGCGGTATGTTCTGGGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGca                                                                                                     <  6:198837/133‑1 (MQ=255)                                  acTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACAt                                                                                                  <  4:402722/135‑1 (MQ=255)                                           ttCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGa                                                                                          <  4:463893/134‑1 (MQ=255)                                           ttCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGa                                                                                          <  3:166128/134‑1 (MQ=255)                                           ttCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGa                                                                                          <  5:442541/134‑1 (MQ=255)                                           ttCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGa                                                                                          <  6:138749/134‑1 (MQ=255)                                           ttCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCAGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGa                                                                                          <  6:70935/134‑1 (MQ=255)                                                    ccTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGc                                                                                 <  5:186987/134‑1 (MQ=255)                                                          taTCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGgccag                                                                          <  8:429492/134‑1 (MQ=255)                                                          aaTCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGgccag                                                                          <  1:125977/135‑1 (MQ=255)                                                                 ggTTGGCGGAATGTTATGTGTCATTTCTGATTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCTGCCAGTGCGGCGACGTCGCGCAGTGCGCACATGAGTTCGATTACGGTGCCGGCCAGCCACggg                                                                   >  8:260655/2‑135 (MQ=255)                                                                                    gtCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTgc                                                >  6:324285/1‑135 (MQ=255)                                                                                     tCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTgcg                                               >  8:205791/1‑135 (MQ=255)                                                                                      cATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTgc                                                                                                        <  7:242546/77‑1 (MQ=255)                                                                                      cATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTgc                                                                                                        >  8:242546/1‑77 (MQ=255)                                                                                      cATTTCTCTTTGGCGTGCATCTGCTTGAAAGCACCGCTGGCGATCAGCTCGCCCACTGCGGCGACCTCGAGCAGTGCGCACATGTCTTCGATTACTGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTgcg                                               >  1:85892/1‑134 (MQ=255)                                                                                        tttCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGTTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCa                                            <  8:364933/135‑1 (MQ=255)                                                                                        tttCTGTTTGGCGGGCATCTGCTTGACCGCACCGCTGGCGATCCGCGCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGGGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGcc                                             <  4:460885/134‑1 (MQ=255)                                                                                        tttCTGTTTGGAGTGCATCTGCTTGACCGCACCGCTGGCGATCAGTTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCa                                            <  8:371226/135‑1 (MQ=255)                                                                                        ttgcTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGGCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGcc                                             <  4:130968/131‑1 (MQ=255)                                                                                         ttCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGa                                                                                                               >  2:351287/1‑67 (MQ=255)                                                                                         ttCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGa                                                                                                               <  1:351287/67‑1 (MQ=255)                                                                                           gtGTTTGGCGGGCATCTGCTTGACCGCCCCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCCCGGGCGCTGGCCACGCTGGCTGCGCCACt                                          <  2:174462/133‑1 (MQ=255)                                                                                             gTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCTATCAGCGCGGCCAGGGCGGCGACGGCTAGCAGTGCGCACATGTGTTCGATTACGGGGCCGGCCCGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGa                                       <  7:209317/135‑1 (MQ=255)                                                                                             gTTTGGCGTGCAGCTGCTTGACCGCACCGCGGGCGAGCCGCTCGGCCAGTGCGGCGACGGCGCGCAGTGCGCACTTGTGTTGGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGa                                       <  1:124045/135‑1 (MQ=255)                                                                                              tttGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGAt                                      >  2:459842/1‑135 (MQ=255)                                                                                                      gCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGCCGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACg                              >  2:449801/1‑135 (MQ=255)                                                                                                      gCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACg                              >  3:291133/1‑135 (MQ=255)                                                                                                                gACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGCGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCCGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGgcg                    >  5:4526/1‑135 (MQ=255)                                                                                                                gACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGgcg                    >  7:410081/1‑135 (MQ=255)                                                                                                                     caccgcTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCag                 >  6:81355/1‑133 (MQ=255)                                                                                                                        cgcTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCagcgagc            >  1:430372/1‑135 (MQ=255)                                                                                                                            ggCGATCAGCGCGGCCCGGGGGGCGACGTCGCGCAGGGCGCACATGTTTTCGATTCGGGTCCCCGCCAGCCACGGGCGCGGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCAGCGAGCGGCt        <  2:409334/135‑1 (MQ=255)                                                                                                                                  ccgCTCGGCCAGGGGGGCGACGTGGAGCAGTGCGAACAGGTGTTCGATTACGGGGCCGCCCAGCCCCGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCAGCGAGCGGCTGACCTg  <  3:347949/133‑1 (MQ=255)                                                                                                                                    gCTCGGCCAGGGCGGCGCCGTCGAGCAGTGCGCACATGTGGTCGATGACGGGGCCGGCCCGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCAGCGAGCGGCTGACCTg  <  3:33630/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    GCAGGATCGGATCTTCGGAACCTTGTGCAGACGACTTTTTCATTCCCCTACCCTCAAAATCCGCCGTTGGCGGTATGTTCTGTGTCATTTCTGTTTGGCGTGCATCTGCTTGACCGCACCGCTGGCGATCAGCTCGGCCAGTGCGGCGACGTCGAGCAGTGCGCACATGTGTTCGATTACGGTGCCGGCCAGCCACGGGCGCTGGCCACGCTGGCTGCGCCACTTGATCTCCGACGGGTCCAGGCGCAGCGAGCGGCTGACCTG  >  ACB122/4917180‑4917443

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