Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,812,541	A→T	6.9%	D85V (GAC→GTC)	CJ019_02665 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,812,541	0	A	T	6.9%	177.4 / 3.4	58	D85V (GAC→GTC)	CJ019_02665	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

CTCAGGCTATGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGACACGAGAGCAGTGGTATC  >  ACB122/2812407‑2812675                                                                                                                                       |                                                                                                                                       cTCAGGCTAAGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTga                                                                                                                                        <  1:139142/135‑1 (MQ=255)    cggCTATGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCCCCAACGGCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGCCTACAGCGCGCACCTCTTAGCGTTCCGGGAAGCCTGTgacga                                                                                                                                     >  6:101381/2‑135 (MQ=255)     ggCTATGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTgacg                                                                                                                                      <  7:255012/133‑1 (MQ=255)      gCTATGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACa                                                                                                                                   <  8:141687/135‑1 (MQ=255)         aTGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGCAGCCTGTGACGACATCg                                                                                                                                >  7:170943/1‑135 (MQ=255)             gCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCt                                                                                                                            >  8:37014/1‑135 (MQ=255)                   cGAGCAGCCCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  2:223572/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGTCGACATCGGTCTTCAGCg                                                                                                                      >  1:244313/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGTCGACATCGGTCTTCAGCg                                                                                                                      >  1:240342/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCcgcg                                                                                                                      >  2:138742/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  6:140591/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  7:56802/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  7:50081/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  7:272288/1‑135 (MQ=255)                   cGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCg                                                                                                                      >  4:65054/1‑135 (MQ=255)                             agacTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAAt                                                                                                            <  4:265454/135‑1 (MQ=255)                              gacTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATg                                                                                                           <  1:138742/135‑1 (MQ=255)                              gacTCTCCCGGTACTTCACCAACGTCCAGCAGAGGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTGCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATg                                                                                                           <  7:197859/135‑1 (MQ=255)                              gacTCTCCCGGTACTTCACCAAAGTCCAGCAGATGTACAACCTGATCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATg                                                                                                           <  7:114716/135‑1 (MQ=255)                                ctctCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGt                                                                                                         <  3:142067/135‑1 (MQ=255)                                  ctcCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTcc                                                                                                       >  3:286494/1‑135 (MQ=255)                                       ggTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCt                                                                                                  >  6:299310/1‑135 (MQ=255)                                                ccAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGTGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCg                                                                                         <  3:264651/135‑1 (MQ=255)                                                ccAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAg                                                                                                                        <  7:24479/104‑1 (MQ=255)                                                ccAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAg                                                                                                                        >  8:24479/1‑104 (MQ=255)                                                ccAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAg                                                                                           >  7:19403/1‑133 (MQ=255)                                                     gTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagacg                                                                                    <  8:19403/135‑1 (MQ=255)                                                     gTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagacg                                                                                    <  5:299310/135‑1 (MQ=255)                                                     gTCCAGCAGATGTACAACCTGTTCGAGGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagacg                                                                                    <  3:155823/135‑1 (MQ=255)                                                      tCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagac                                                                                     >  5:103780/1‑133 (MQ=255)                                                      tCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagac                                                                                     >  4:147534/1‑133 (MQ=255)                                                      tCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagac                                                                                     >  6:189437/1‑133 (MQ=255)                                                          gcagATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCgagacggagac                                                                               <  5:189437/135‑1 (MQ=255)                                                                gTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGGTCCGGGAAGCCTGTGACGACATCGGTCGTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTg                                                                         <  8:179943/135‑1 (MQ=255)                                                                  aCAACCTGTTCGATGATCGAGAGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGAGACGACATCGGTCATCAGCGTTCACCCAATGGACCTGTCTGCATGAGACAGAAGGAGACGAACTTGc                                                                       >  6:169125/1‑135 (MQ=255)                                                                    aaCCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTc                                                                     <  8:56802/135‑1 (MQ=255)                                                                           tCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGTCGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcat                                                               <  2:240342/134‑1 (MQ=255)                                                                           tCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGTCGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcat                                                               <  2:244313/134‑1 (MQ=255)                                                                           tCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcatc                                                              >  8:84924/1‑135 (MQ=255)                                                                           tCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcat                                                               <  6:244730/134‑1 (MQ=255)                                                                            cGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGCCGGAGACGTACTGGCTCAGtcatca                                                             >  7:170581/1‑135 (MQ=255)                                                                             gatgatCGAGTGCCGTACGCCTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGAAGACATCGGTCTTCAGCGTTCACCCAATGGTCATGTCTGCATGAGCGAGACGGAGACGTACTGCCTCAGtcatca                                                             >  6:183895/1‑134 (MQ=255)                                                                             gatgatCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTg                                                                         >  1:111522/1‑122 (MQ=255)                                                                             gatgatCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTg                                                                         <  2:111522/122‑1 (MQ=255)                                                                             gatgatCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcatca                                                             >  7:81658/1‑134 (MQ=255)                                                                             gatgatCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcatca                                                             <  1:223572/134‑1 (MQ=255)                                                                             gatgatCGAGTGCCGTAAGACTACAGCGAGCACCTGTTAGCGTTCCGTGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcatca                                                             <  8:125808/134‑1 (MQ=255)                                                                             gatgatAGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTg                                                                         <  4:148539/122‑1 (MQ=255)                                                                             gatgatAGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCCGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTg                                                                         >  3:148539/1‑122 (MQ=255)                                                                              atgatCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGGGACGACATCGGGCTTCAGCGTTCACCCAATGGTCATGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGtcctct                                                             >  6:164156/1‑132 (MQ=255)                                                                                   cGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGc                                                      <  1:261150/135‑1 (MQ=255)                                                                                           gTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGt                                              <  8:272288/135‑1 (MQ=255)                                                                                             aCGACTACAGCGAGCACGTGTTAGCTTTCCGGGAAGCGTGTTACGAAATCGGTCTTCAGCGTTCACCAAATGGTCCTGTCTGCATGAGCGAGACGTAGACGTAGTGGCTCTGTAATCACCATAGAATGAACGTcc                                            <  8:81658/135‑1 (MQ=255)                                                                                               gACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCAc                                                                                                                 >  2:95150/1‑64 (MQ=255)                                                                                               gACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCAc                                                                                                                 <  1:95150/64‑1 (MQ=255)                                                                                               gACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTg                                          >  7:245470/1‑135 (MQ=255)                                                                                                     aGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGc                                     <  4:286494/134‑1 (MQ=255)                                                                                                              cTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGt                           <  3:147534/135‑1 (MQ=255)                                                                                                                 ttAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGaa                        >  4:232445/1‑135 (MQ=255)                                                                                                                   aGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAAc                       >  5:274092/1‑134 (MQ=255)                                                                                                                     cGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGa                    <  7:37014/135‑1 (MQ=255)                                                                                                                                      gacgacATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGACACGAGAGCAGTGGTAt   >  4:298705/1‑135 (MQ=255)                                                                                                                                      gacgacATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGACACGAGAGCAGTGGTAt   <  4:207168/135‑1 (MQ=255)                                                                                                                                       acgacATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGACACGAGAGCAGTGGTATc  >  5:215132/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                       CTCAGGCTATGAGCAGGTCGAGCAGACCAGACTCTCCCGGTACTTCACCAACGTCCAGCAGATGTACAACCTGTTCGATGATCGAGTGCCGTACGACTACAGCGAGCACCTGTTAGCGTTCCGGGAAGCCTGTGACGACATCGGTCTTCAGCGTTCACCCAATGGTCCTGTCTGCATGAGCGAGACGGAGACGTACTGGCTCAGTCATCACCAGAGCATGAACGTCCTGACCGCACGGATCCGTGAACTGACACGAGAGCAGTGGTATC  >  ACB122/2812407‑2812675

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