Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,859,935	C→T	6.4%	E143K (GAG→AAG)	thyA ←	Thymidylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,859,935	0	C	T	6.4%	202.5 / 2.7	62	E143K (GAG→AAG)	thyA	Thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base C (22/36);  new base T (2/2);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

TCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGATCTCGCCCAGGTCGT  >  ACB122/5859805‑5860063                                                                                                                                   |                                                                                                                                 tCGTCCAGTTGCGCGCAGGTCCAGCCGTGGAACAGGATGCGGCGGCTACCGGGGTCGTGGGCGATGGGGGCCAGGCACTGGCGTACCTGGGCGAGCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTc                                                                                                                               <  4:127674/134‑1 (MQ=255) tCGTCATGTTGCGCGCTGTTCCAGCCGTGGAACTGGATGCGGCGGCTACCCGTGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACTCGACGCCGTCCTCTTc                                                                                                                               <  8:138989/134‑1 (MQ=255)            cgcgcAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGCTGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCGTCGTCCGTGGCAAt                                                                                                                   >  4:349356/1‑135 (MQ=255)            cgcgcAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGGCCTCTTCGTCCTGGGCAAt                                                                                                                   >  6:428146/1‑135 (MQ=255)                             ggAACAGGCTGCGGCGGCGACCCGGGGCGTGGGCGAGGGTGGCCTGGCACGGGCGTACCGGGGCGAGCGCCTTGTCCAGGATTACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCtt                                                                                                  <  2:142586/135‑1 (MQ=255)                                aCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGGGTACCTGGTCGATCGCCTTGTACAGGCTGACGAACGCGCCGCCGGCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCt                                                                                               >  3:79350/1‑135 (MQ=255)                                aCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTTTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCt                                                                                               >  8:259461/1‑135 (MQ=255)                                aCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTAACTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTTTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCt                                                                                               >  8:266171/1‑135 (MQ=255)                                         gggggcTCCCGGGGTCGTGGCCGATGGTGGCCAGGCACGGGCGTAGCTGGTCGATCGCCTTGTACAGGATGACTACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAATCCGGCTTGCTCGGCCAttt                                                                                      <  2:50831/132‑1 (MQ=255)                                          ggggcTACCGGGGGCGTGGGCGAGGGGGGCCAGCCAGTGGCGTACGGGGGCGATCGCCTTTAACCGGAGCACCACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTc                                                                                     <  5:428146/133‑1 (MQ=255)                                             gcTCCCCGGGTCGTGGGCGATGGGGTCCAGGCACGGGCTTAGCTGGGGGATGGCGTTGTACAGTATGTCGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGAt                                                                                  <  6:176998/135‑1 (MQ=255)                                             gcTACCCGGGTCGTGGGCGATGGTGGCCAGGCCCTGGCGTACCTGGGCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGAt                                                                                  <  3:462551/135‑1 (MQ=255)                                               tACCCGGGTCGGTGGCGAGGGTGGCCCGGCACGGGCGTACGGGGTCGATCGCCTTGTCCCGGATGACGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATgg                                                                                <  3:349356/135‑1 (MQ=255)                                                 cccGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTTCCTGGTCGATCGCCATGTAAATGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAAAGCGGCGGAAAACGGCATGCTCGGCCAGTTCGATGgcc                                                                              >  2:130255/1‑135 (MQ=255)                                                 cccGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGCACGCGCCGCCGTCCTCGTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGGCATTTCGAGGgcc                                                                              >  7:448003/1‑135 (MQ=255)                                                    ggggCGGTGGCGAGGGGGCCCAGGCACTGGCGTCCCTGGGCGATCGCCTTGTACCGCAGGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGgccgcc                                                                           <  3:224659/135‑1 (MQ=255)                                                    gggTGGTTGGCGATGGTGTCCAGGCAGTGGCGTAGCTGGTCGAGCGCCTTGTCCAGGATGACGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGgccgcc                                                                           <  3:441303/135‑1 (MQ=255)                                                                 tGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCg                                                              <  6:212906/135‑1 (MQ=255)                                                                   gtgtCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            <  3:390334/135‑1 (MQ=255)                                                                   gtgtCCAGGCACTGGCGTACCGGGTCGATCGCCTTGTACAGTCTGCCGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            <  7:35529/135‑1 (MQ=255)                                                                    tgtCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGTATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            <  7:8762/134‑1 (MQ=255)                                                                    tgtCCAGGCACTGGCGTACCGGGTCGAGCGCCTTGTACAGTATGACCAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            <  7:42811/134‑1 (MQ=255)                                                                          ggCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCt                                                     >  1:311983/1‑135 (MQ=255)                                                                                   gTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCg                                            >  4:257649/1‑135 (MQ=255)                                                                                   gTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCg                                            >  7:253516/1‑135 (MQ=255)                                                                                   gTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCg                                            >  8:430469/1‑135 (MQ=255)                                                                                   gTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCg                                            >  2:286488/1‑135 (MQ=255)                                                                                     aCCTGGGCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGc                                          <  6:211021/135‑1 (MQ=255)                                                                                        tGGTCGAGCGCCTGGAACCGGAGCACCAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCAc                                       <  2:415550/135‑1 (MQ=255)                                                                                          ggcGATCGCCTGGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTg                                     <  4:112465/133‑1 (MQ=255)                                                                                              aTCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCGTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCGTGTAGCCCGGCCACTGgcgc                                 >  1:215359/1‑135 (MQ=255)                                                                                                 gCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAt                              <  8:185808/135‑1 (MQ=255)                                                                                                  cgttGTACTGGATCCCCACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  1:158738/133‑1 (MQ=255)                                                                                                  ccTTGTACTGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  6:207200/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGCTGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  3:164861/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGTCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGTCGCCAtt                             <  8:364515/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAACGCGACGCCGTCCTTTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  7:266171/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAACGCGACGCCGTCCTTTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  7:259461/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  3:98992/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  4:86279/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  4:228564/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGATGACGAAAGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTAGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             >  6:25478/1‑135 (MQ=255)                                                                                                  ccTTGTACAGGATCACGATCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCGATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  5:95497/135‑1 (MQ=255)                                                                                                  ccTTGTACAGGAGGCCCAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  5:20806/135‑1 (MQ=255)                                                                                                        aCAGGATGACGACCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGcaca                        <  1:286488/134‑1 (MQ=255)                                                                                                            gATGCCGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATa                   <  3:389001/135‑1 (MQ=255)                                                                                                               gACGAACGCGACGCCTTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            <  3:462179/91‑1 (MQ=255)                                                                                                               gACGAACGCGACGCCTTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCggg                                                            >  4:462179/1‑91 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCTGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             >  6:16438/1‑122 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCTGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  5:16438/122‑1 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGgccgcc                                                                           >  3:352017/1‑76 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGgccgcc                                                                           <  4:352017/76‑1 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             <  3:223797/122‑1 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCAtt                             >  4:223797/1‑122 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGa                 >  5:217961/1‑134 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGAt                >  1:372166/1‑135 (MQ=255)                                                                                                               gACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGAt                >  5:170086/1‑135 (MQ=255)                                                                                                                    aCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGATCTCGc           >  3:290714/1‑135 (MQ=255)                                                                                                                    aCGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGATCTCGc           >  2:266651/1‑135 (MQ=255)                                                                                                                         aCGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATgg                                                                                <  2:356408/61‑1 (MQ=255)                                                                                                                         aCGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATgg                                                                                >  1:356408/1‑61 (MQ=255)                                                                                                                             cGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGATCTCGCCCAGgtcgt  <  1:307505/135‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 TCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCGCCATTGCACACCATAGATCTCGCCCAGGTCGT  >  ACB122/5859805‑5860063

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