Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,847,475	A→T	6.4%	Q62L (CAG→CTG)	CJ019_02703 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,847,475	0	A	T	6.4%	200.8 / 3.3	62	Q62L (CAG→CTG)	CJ019_02703	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GGCTCGGCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCGATCACCAGGCTGGTCGAG  >  ACB122/2847346‑2847605                                                                                                                                  |                                                                                                                                   ggCTCGGCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATTCCGACCTTGCCCAGCCCACGCCGGCTGTAGCACAGATcg                                                                                                                               >  2:155250/1‑135 (MQ=255) ggCTCGGCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATcg                                                                                                                               >  7:253866/1‑135 (MQ=255)       gCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACTGATCGCGGATc                                                                                                                         >  1:268195/1‑135 (MQ=255)       gCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACTGATCGCGGATc                                                                                                                         >  1:260993/1‑135 (MQ=255)          aGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCggg                                                                                                                      >  1:276018/1‑135 (MQ=255)          aGCCAGGGGACACGCTGAGCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCggg                                                                                                                      <  3:387904/135‑1 (MQ=255)             cAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCAc                                                                                                                   <  4:78563/135‑1 (MQ=255)             cAGGGGACACGCCGATCGGCGTATCACCCAAACCGCCGTCACCCTTGGCCCTTGCCCCGGGCATCGACCTGCAGGGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCAc                                                                                                                   <  5:350172/135‑1 (MQ=255)                 ggACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGc                                                                                                               <  2:276018/135‑1 (MQ=255)                   acacGCTCAGCGGCGTACCACGAAAACAGCCGTCAGACTGGCCCCTTGCCCCGGGCATCGCCCGCCCGTGGGCCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCg                                                                                                             <  3:463210/135‑1 (MQ=255)                        cTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGCCCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTg                                                                                                        >  7:122897/1‑135 (MQ=255)                        cTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGCCCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTg                                                                                                        >  6:169469/1‑135 (MQ=255)                        cTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGCCCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTg                                                                                                        >  2:334293/1‑135 (MQ=255)                        cTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCt                                                                                                         >  7:347136/1‑134 (MQ=255)                         tGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGCCCTGCAGTGGGTCGATGCCGCCCGTGCCAAGCCACCGGCGGCGGGAGCACAGAGCGCGGCTCGGGCCCCGGCCgcgcggg                                                                                                       >  1:455521/1‑129 (MQ=255)                         tGAGCGGCGTATCACCCAAACAGCCGGCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGGGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTgg                                                                                                       <  3:399251/135‑1 (MQ=255)                           agcGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCAGCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGAt                                                                                                     <  6:198691/133‑1 (MQ=255)                           aTCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGGGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGAt                                                                                                     <  2:285758/135‑1 (MQ=255)                              ggCGTATCACGAAAACAGCCTTCATACTTGGCCCTTGCCCCGGGCATCGACCTGCTGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTc                                                                                                  <  7:217079/135‑1 (MQ=255)                                   aTCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGAt                                                                                             <  3:20195/135‑1 (MQ=255)                                        gCAAACCGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACa                                                                                        <  7:198293/135‑1 (MQ=255)                                        gCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGCCCTTGCCAAGCCCACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGCTCTGGATGTCCTGATCTACa                                                                                        >  7:78369/1‑135 (MQ=255)                                          aaaCAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAAc                                                                                      <  2:69057/135‑1 (MQ=255)                                             cAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGtt                                                                                   >  7:454619/1‑135 (MQ=255)                                               gCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGCCCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACCACGTTg                                                                                  >  6:242872/1‑134 (MQ=255)                                                    cAGACTGGGCCCTTGCCCCGGGCATCGACGTGCAGGGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCt                                                                            <  7:291426/135‑1 (MQ=255)                                                     aGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGCTCGCGGATCGGGCCCCTGCCGAGCTGGATGTGCTGATCTACAACGTTGGCATcag                                                                           >  1:96348/1‑133 (MQ=255)                                                             cccTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAac                                                                   >  8:238648/1‑135 (MQ=255)                                                               cTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAacac                                                                 >  4:435843/1‑135 (MQ=255)                                                                ttGCCCCGGGCATCGTCCTGCAGGGGGGGGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGAATCTGGGAAGAACACg                                                                <  1:338928/135‑1 (MQ=255)                                                                ttGCCCCGGGCATCGCCCTGCAGTGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACg                                                                <  5:184238/135‑1 (MQ=255)                                                                tggCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACg                                                                <  6:68835/133‑1 (MQ=255)                                                                    cccGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCt                                                             >  7:109100/1‑134 (MQ=255)                                                                     ccGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTc                                                           <  4:216508/135‑1 (MQ=255)                                                                      cGGCCAGCGACCGGCCGTGGGGCTAGGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCa                                                          <  5:71833/135‑1 (MQ=255)                                                                         gCATCGACCTGCAGTGGGTCGATGACGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCg                                                       <  2:323310/135‑1 (MQ=255)                                                                          cATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGa                                                      <  2:48913/135‑1 (MQ=255)                                                                          cATCGACCTGCAGGGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGa                                                      <  7:390936/135‑1 (MQ=255)                                                                            tCGCCCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACTGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGc                                                    <  2:268195/135‑1 (MQ=255)                                                                            tCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACTGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGc                                                    <  2:260993/135‑1 (MQ=255)                                                                               aCCTGCAGTGGGGCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCa                                                   <  6:456352/133‑1 (MQ=255)                                                                                  tGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATg                                              <  8:263552/135‑1 (MQ=255)                                                                                     aGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGCTCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCAt                                           >  1:411540/1‑135 (MQ=255)                                                                                       tGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACg                                                                >  4:255153/1‑112 (MQ=255)                                                                                       tGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACg                                                                <  3:255153/112‑1 (MQ=255)                                                                                        ggggCGATGCCGACGTTGCCAAGCCAACGGGGGCTGTAGCACAGTTCGCGGTTGGGGCACCTGCCTATCTGTATGTCGTTTTCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTc                                         <  4:172716/134‑1 (MQ=255)                                                                                            cGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACCACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTga                                    >  8:276955/1‑135 (MQ=255)                                                                                            cGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTga                                    >  4:397383/1‑135 (MQ=255)                                                                                            cGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTga                                    >  3:172716/1‑135 (MQ=255)                                                                                            cGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTga                                    >  3:445086/1‑135 (MQ=255)                                                                                            cGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTga                                    >  5:69369/1‑135 (MQ=255)                                                                                                 ccGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACa                               >  1:305520/1‑135 (MQ=255)                                                                                                 ccGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACa                               >  3:352890/1‑135 (MQ=255)                                                                                                          ccatgccacGGGGGGCTGTTGCACAGACGGCGGATCGGGCACGTGCGTATCTGTATGTCCTGATCTCCACCGTGGGCATCTGGGAAGAACACGCCTGCAGCGAGCACTATGCATTCCTTGACGACAGcgacga                        <  2:455521/133‑1 (MQ=255)                                                                                                          ccaagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGcgacga                        <  3:269460/133‑1 (MQ=255)                                                                                                          ccaagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGaa                       <  8:156329/134‑1 (MQ=255)                                                                                                            aagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCt                                                                                           >  3:207076/1‑64 (MQ=255)                                                                                                            aagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCt                                                                                           <  4:207076/64‑1 (MQ=255)                                                                                                            aagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCg                    >  2:198825/1‑135 (MQ=255)                                                                                                            aagccaaCGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCg                    >  1:415453/1‑135 (MQ=255)                                                                                                                      cggcTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGTAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCGATCACCAGGc          <  4:399241/135‑1 (MQ=255)                                                                                                                      cggcTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCGATCACCAGGc          >  4:133209/1‑135 (MQ=255)                                                                                                                               caccGATCGCGGAGCGGGCACCTGCCGATCTGGATGTCCTGATCTACAAGGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCGATCACCAGGCTGGTCGAg  <  5:115918/134‑1 (MQ=255)                                                                                                                                  |                                                                                                                                   GGCTCGGCAAGCCAGGGGACACGCTGATCGGCGTATCACGCAAACAGCCGTCAGACTTGGCCCTTGCCCCGGGCATCGACCTGCAGTGGGTCGATGCCGACCTTGCCAAGCCAACGGCGGCTGTAGCACAGATCGCGGATCGGGCACCTGCCGATCTGGATGTCCTGATCTACAACGTTGGCATCTGGGAAGAACACGCCTTCAGCGAGCACTATGCATTCCTTGACGACAGCGACGAATCGATCACCAGGCTGGTCGAG  >  ACB122/2847346‑2847605

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