Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,851,374	A→T	6.9%	F153Y (TTC→TAC)	CJ019_04550 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,851,374	0	A	T	6.9%	168.8 / 3.1	58	F153Y (TTC→TAC)	CJ019_04550	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base T (2/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGTAC  >  ACB122/4851251‑4851506                                                                                                                            |                                                                                                                                     tgCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCcaca                                                                                                                           <  3:315590/135‑1 (MQ=255) tgCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCcaca                                                                                                                           <  3:233520/135‑1 (MQ=255)  gCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCCACTCCCTGACGTCACACTAAAGCAAGTAAACCCGCACTCAAGTACAACCCCacac                                                                                                                          >  7:46884/1‑135 (MQ=255)  gCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGTACAACCCCacac                                                                                                                          >  7:40567/1‑135 (MQ=255)  gCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCcaca                                                                                                                           <  5:204586/134‑1 (MQ=255)           ttACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGCCCCAACTCCCTGCCGTCACCCTAAAGCAAGTAAACCCGCACTCAAGAACAACCCAACACCAAACACaa                                                                                                                 >  7:359812/1‑135 (MQ=255)                   gcgcAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTgg                                                                                                        >  3:221838/1‑136 (MQ=255)                   gcgcAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTg                                                                                                         >  2:374883/1‑135 (MQ=255)                   gcgcAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCAACACCAAACACAAAATGCGTg                                                                                                         >  5:32957/1‑135 (MQ=255)                      cAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGGACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGcc                                                                                                      <  5:174355/135‑1 (MQ=255)                      cAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGcc                                                                                                      <  1:233437/135‑1 (MQ=255)                      cAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGcc                                                                                                      <  2:225020/135‑1 (MQ=255)                                aGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAg                                                                                           >  7:302028/1‑136 (MQ=255)                                 ggggTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAg                                                                                           >  1:296090/1‑135 (MQ=255)                                        tgctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAGCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAc                                                                                   >  8:239784/1‑136 (MQ=255)                                        tgctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGCACAACCCCACCCCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAc                                                                                   >  8:63807/1‑136 (MQ=255)                                        tgctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTa                                                                                    >  5:141546/1‑135 (MQ=255)                                        tgctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTa                                                                                    >  1:283291/1‑135 (MQ=255)                                          ctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACACCCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAcg                                                                                  >  1:2840/1‑135 (MQ=255)                                          ctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCcaca                                                                                                                           <  5:390698/94‑1 (MQ=255)                                          ctgTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCcaca                                                                                                                           >  6:390698/1‑94 (MQ=255)                                             ttAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCACCTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTAcgc                                                                                 >  2:296153/1‑133 (MQ=255)                                               aGCAGGCATGCTGATTTTGCGCTCGCTCTGCCCCAACTCCCTGCCGTCACACTAAAGCATGTAAACCCGCACTCAAGAACGACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGAAg                                                                             >  1:328427/1‑135 (MQ=255)                                               aGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAg                                                                             >  1:269964/1‑135 (MQ=255)                                                gCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAggc                                                                           >  3:223159/1‑136 (MQ=255)                                                  aGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAggcgt                                                                         >  7:352822/1‑136 (MQ=255)                                                       tGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAACCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGcc                                                                                                      <  1:291224/102‑1 (MQ=255)                                                       tGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGcc                                                                                                      >  2:291224/1‑102 (MQ=255)                                                       tGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAggcgtcggc                                                                     >  1:392751/1‑135 (MQ=255)                                                          tGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTcc                                                                 >  3:109668/1‑136 (MQ=255)                                                                   gctcgctcTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGa                                                          <  5:74581/134‑1 (MQ=255)                                                                        ctctGACACAACTCCCTGCCGTCACTCTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGGCGGCGTCCTTGACGAAAAGAAc                                                   <  8:359812/136‑1 (MQ=255)                                                                                    tCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGa                                                          <  3:334176/117‑1 (MQ=255)                                                                                    tCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGa                                                          >  4:334176/1‑117 (MQ=255)                                                                                      ccTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAAc                                                   <  6:321841/122‑1 (MQ=255)                                                                                      ccTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAAc                                                   >  5:321841/1‑122 (MQ=255)                                                                                      ccTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgc                                     >  8:24361/1‑136 (MQ=255)                                                                                       cTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgc                                     <  2:283291/135‑1 (MQ=255)                                                                                       cTGCAGTCACACTATATCAAGTAAACCCGCACTCAAGAAGAACCCCACACCAAACTCAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCTCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgc                                     <  6:374372/135‑1 (MQ=255)                                                                                          ccGTCACACTAAAGCAAGTACACCCGCACTCAAGAACAACCCCACCCCAAACACAAAATGCGTGGCCCGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGACCCCTGCCCCCAACgctgg                                  >  1:79473/1‑135 (MQ=255)                                                                                          ccGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAAc                                                   >  4:174293/1‑118 (MQ=255)                                                                                          ccGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAAc                                                   <  3:174293/118‑1 (MQ=255)                                                                                          ccGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgctgg                                  >  5:357966/1‑135 (MQ=255)                                                                                           cGTCAAACTAAAGCAAGGAAACCCGCACTCAACAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCgctggct                                <  4:223159/136‑1 (MQ=255)                                                                                                        cAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCa                                                                                                     <  1:223474/54‑1 (MQ=255)                                                                                                        cAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCa                                                                                                     >  2:223474/1‑54 (MQ=255)                                                                                                              aaCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGTCAGTTACGCGGAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGgcgc             >  8:383408/1‑136 (MQ=255)                                                                                                              aaCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCTCATCGCTGGCTGCATGAAGCACAGCGgcgc             >  3:359415/1‑136 (MQ=255)                                                                                                              aaCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAgcggcg              <  4:359415/135‑1 (MQ=255)                                                                                                              aaCCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAgcggcg              <  4:109668/135‑1 (MQ=255)                                                                                                                 ccGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCGGGCTGCATGAAGCACAGCGGCGCCa           <  2:255996/135‑1 (MQ=255)                                                                                                                  cGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCAc          >  6:346044/1‑135 (MQ=255)                                                                                                                      cTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCt      >  2:88509/1‑135 (MQ=255)                                                                                                                         aacaacaacCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGt    <  8:371704/131‑1 (MQ=255)                                                                                                                         aaGTACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGGCCTGGACGAAAAGAACCCTGCCCCCCTCGCGGGCTGCATGAAGCACAGCGGCGCCACGGCTGt    <  8:46884/134‑1 (MQ=255)                                                                                                                         aaGTACAACCCCACACCAAACACAAAATGCGGGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGt    <  8:40567/134‑1 (MQ=255)                                                                                                                          aGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGTac  <  1:374883/135‑1 (MQ=255)                                                                                                                          aGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGTa   <  5:346044/134‑1 (MQ=255)                                                                                                                          aGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGTa   <  2:44154/134‑1 (MQ=255)                                                                                                                            |                                                                                                                                     TGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAACCCCACACCAAACACAAAATGCGTGGCCAGGCTCTTCAGGCAGTTACGCGCAGGCGTCGGCGTCCTGGACGAAAAGAACCCTGCCCCCATCGCTGGCTGCATGAAGCACAGCGGCGCCACGGCTGTAC  >  ACB122/4851251‑4851506

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