Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,809,343	G→A	6.3%	P53P (CCG→CCA)	surE →	5'‑nucleotidase SurE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,809,343	0	G	A	6.3%	205.1 / 3.6	63	P53P (CCG→CCA)	surE	5'‑nucleotidase SurE
Reads supporting (aligned to +/- strand):  ref base G (20/39);  new base A (2/2);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

GTGTTACCGCACCTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGCCAACCTGG  >  ACB122/1809213‑1809473                                                                                                                                   |                                                                                                                                   gggtTACCGCACCGGGCATCGCCGCGCTGCACGCTGCGCGGGCGGATTATGCCGCGTGCGCGGTGATCGCCCCGCATCAAGACACGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGa                                                                                                                                <  5:34262/133‑1 (MQ=255)        cGCACCTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTgg                                                                                                                          <  3:72161/134‑1 (MQ=255)        cGCACCTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTgg                                                                                                                          <  6:2715/134‑1 (MQ=255)        cGCACCTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTgg                                                                                                                          <  3:282427/134‑1 (MQ=255)          cACCTGGCATCGCCGCGGTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGAGATCCCCCCGGATCAAGACAAGAGCGGCGCCACCAGTTCGCTGCCGCTGGACCGGCCACTGCATCCGCAGACGCTGccc                                                                                                                        >  7:214858/1‑134 (MQ=255)            ccTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAAt                                                                                                                     <  4:350950/135‑1 (MQ=255)             cTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGCCCCTGGCCAATg                                                                                                                    >  1:254618/1‑135 (MQ=255)                 cATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCtt                                                                                                                <  1:274784/135‑1 (MQ=255)                  ttCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTtc                                                                                                               <  4:126685/134‑1 (MQ=255)                  aTCGCCGCGCTGCACGCTGCGCTGGCGGTTTTTGCGGAGTGCGCGGTGATCGCCCCGGAGCAAGCCAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTtc                                                                                                               <  4:243868/135‑1 (MQ=255)                  aTCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTtc                                                                                                               >  2:206032/1‑135 (MQ=255)                               aCGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACg                                                                                                  <  4:40236/135‑1 (MQ=255)                               aCGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACg                                                                                                  <  7:86066/135‑1 (MQ=255)                               aCGCTGCGCTGGCGGATTATGCCGAGGGCGCGGTGAGCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACg                                                                                                  <  5:6038/135‑1 (MQ=255)                                 gctgcgctgGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGc                                                                                                <  5:134726/135‑1 (MQ=255)                                   tgcgctgGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCAc                                                                                              <  8:219433/135‑1 (MQ=255)                                           cGGATTATGCCGTGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGccgaccg                                                                                      <  2:13598/135‑1 (MQ=255)                                            ggATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGccgaccga                                                                                     <  6:169385/135‑1 (MQ=255)                                              aTTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACt                                                                                   <  2:90470/135‑1 (MQ=255)                                              aTTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACt                                                                                   <  6:211482/135‑1 (MQ=255)                                              aTTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACt                                                                                   <  4:131684/135‑1 (MQ=255)                                                tATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTg                                                                                  <  2:254618/134‑1 (MQ=255)                                                tATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGATGACGCTGGACCGGCCACTGCATCCGCAGACCCTCGCCAATGGCTGCATCAGCCTCAACGGCACGCCGACCGACTGc                                                                                 >  2:24325/1‑135 (MQ=255)                                                     cGAGGGCGCGGTGAGCGCCCCGGATCAAGACAAGTGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCGACTGCATCCGCAGACCCTGGCCACTGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGc                                                                             <  5:344442/134‑1 (MQ=255)                                                        gggcgcGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCt                                                                         <  6:352421/133‑1 (MQ=255)                                                        gTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCt                                                                         <  4:70495/135‑1 (MQ=255)                                                         tGCGCGGTGATCGCCCCGGATCAAGACAAGCGCGGCACCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTg                                                                        <  1:229282/135‑1 (MQ=255)                                                          gcgcGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCCCGCCGACCGACTGCGTGCACCTgg                                                                       >  5:161043/1‑135 (MQ=255)                                                           cgcgGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCACAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTggg                                                                      >  1:368600/1‑135 (MQ=255)                                                           cgcgGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCACAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTggg                                                                      >  1:376334/1‑135 (MQ=255)                                                                   tCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATg                                                              <  5:281326/135‑1 (MQ=255)                                                                     gCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATggg                                                            >  4:347799/1‑135 (MQ=255)                                                                      ccccGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCACCGGCACGCCGACCGCCGGCGTGCACCTGGGGCGCCATGGGc                                                           >  4:68178/1‑135 (MQ=255)                                                                            aTCAAGACCAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGc                                                     <  1:249679/135‑1 (MQ=255)                                                                             tCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAg                                                                                                                                 <  8:348403/58‑1 (MQ=255)                                                                             tCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAg                                                                                                                                 >  7:348403/1‑58 (MQ=255)                                                                             tCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCACAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGcc                                                    <  2:376334/135‑1 (MQ=255)                                                                              cacgacaagaGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCACAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGcc                                                    <  2:368600/131‑1 (MQ=255)                                                                                    aagaGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATg                                             >  1:288139/1‑135 (MQ=255)                                                                                        gcggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCgg                                         <  7:222604/135‑1 (MQ=255)                                                                                         cggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCCTCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGGTGcc                                                    >  5:291167/1‑123 (MQ=255)                                                                                         cggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCg                                                                                <  2:341163/95‑1 (MQ=255)                                                                                         cggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCg                                                                                >  1:341163/1‑95 (MQ=255)                                                                                         cggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGcc                                                    <  6:291167/123‑1 (MQ=255)                                                                                         cggcgCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCgg                                         >  1:314173/1‑134 (MQ=255)                                                                                              ccAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCACTGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATgg                                   >  2:221779/1‑135 (MQ=255)                                                                                              ccAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATgg                                   >  1:236541/1‑135 (MQ=255)                                                                                               cagcagTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGt                                  >  1:8686/1‑135 (MQ=255)                                                                                                   agTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCgtg                              <  5:313447/135‑1 (MQ=255)                                                                                                   agTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCgtg                              <  5:210649/135‑1 (MQ=255)                                                                                                   agTTAGCTGTCGCTGGACCGGCCTCTGCATCCGCTGACCCTGGCCAATGGCTTCATGAGCCTCAACGGGACGCCGACCGACTGCGAGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCgtg                              <  5:313442/135‑1 (MQ=255)                                                                                                    gTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCgtgt                             >  8:231787/1‑135 (MQ=255)                                                                                                         cTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCgggg                        <  7:245014/135‑1 (MQ=255)                                                                                                              gcgggACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCACCGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCaa                   <  2:7538/132‑1 (MQ=255)                                                                                                              gCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCaa                   <  3:347799/135‑1 (MQ=255)                                                                                                              gCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCaa                   <  7:75886/135‑1 (MQ=255)                                                                                                              gCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCaa                   <  3:311899/135‑1 (MQ=255)                                                                                                                     cGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGt            <  7:225660/135‑1 (MQ=255)                                                                                                                     cGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGt            >  7:337576/1‑135 (MQ=255)                                                                                                                        ccACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGcc         <  4:44258/135‑1 (MQ=255)                                                                                                                          aCTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGCCaa       >  6:191120/1‑135 (MQ=255)                                                                                                                          aCTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGCCaa       >  5:283841/1‑135 (MQ=255)                                                                                                                               aTCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGCCAACCTgg  >  4:82233/1‑135 (MQ=255)                                                                                                                                   |                                                                                                                                   GTGTTACCGCACCTGGCATCGCCGCGCTGCACGCTGCGCTGGCGGATTATGCCGAGTGCGCGGTGATCGCCCCGGATCAAGACAAGAGCGGCGCCAGCAGTTCGCTGACGCTGGACCGGCCACTGCATCCGCAGACCCTGGCCAATGGCTTCATCAGCCTCAACGGCACGCCGACCGACTGCGTGCACCTGGGGCTCAATGGGCTGTTGCCGCAGATGCCGGACATGGTCGTGTCGGGGATCAACCTCGGTGCCAACCTGG  >  ACB122/1809213‑1809473

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