Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,255,382	T→C	7.0%	Y204C (TAC→TGC)	yjcD ←	putative ATP‑dependent DNA helicase (sporulation associated)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,255,382	0	T	C	7.0%	163.8 / 4.6	57	Y204C (TAC→TGC)	yjcD	putative ATP‑dependent DNA helicase (sporulation associated)
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (29/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 = 2.46e-01

TTGTTCCCGCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCTGTCAGTACCCGTGTTTTTCCGCTCC  >  BS168_reference_genome/1255246‑1255520                                                                                                                                         |                                                                                                                                           ttGTTCCCGCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTa                                                                                                                                           <  1:265897/138‑1 (MQ=255)    ttCCCGCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTgg                                                                                                                                       >  4:205100/1‑139 (MQ=255)      cccGCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGt                                                                                                                                      <  3:205100/138‑1 (MQ=255)         gCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGc                                                                                                                                  >  3:136159/1‑139 (MQ=255)                cTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGc                                                                                                                           <  8:175242/139‑1 (MQ=255)                   cATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGctct                                                                                                                        >  6:41472/1‑139 (MQ=255)                     ttttCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCtt                                                                                                                       >  1:278895/1‑138 (MQ=255)                                tGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGa                                                                                                           >  6:194836/1‑139 (MQ=255)                                              tttAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTtgttgt                                                                                              >  2:101596/1‑138 (MQ=255)                                                      aaTTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGATCTTTCATTTCATCGACAGCTTTTGTTGTAAAGGTGa                                                                                      >  1:46455/1‑138 (MQ=255)                                                            agTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGt                                                                                                                                      <  2:148918/84‑1 (MQ=255)                                                            agTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGt                                                                                                                                      >  1:148918/1‑84 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGt                                                                                                                                      >  4:131991/1‑80 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGt                                                                                                                                      <  3:131991/80‑1 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCa                                                                                                                                >  1:237014/1‑86 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCa                                                                                                                                <  2:237014/86‑1 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCata                                                                           >  4:92037/1‑139 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCata                                                                           >  5:222442/1‑139 (MQ=255)                                                                aTACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCata                                                                           >  3:130508/1‑139 (MQ=255)                                                                          ttATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGTCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCg                                                                 <  5:41472/139‑1 (MQ=255)                                                                           tataAAACAGACTGTGGAATGTTCCTGTCACTATTCTTCTGACTTTTGCCGGCTGTAATCCGTCCTGGTTCGCGATACGCTGTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCg                                                                 <  1:10479/138‑1 (MQ=255)                                                                           tataAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    >  6:209706/1‑103 (MQ=255)                                                                           tataAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    <  5:209706/103‑1 (MQ=255)                                                                            ataAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCg                                                                 <  4:77673/137‑1 (MQ=255)                                                                            ataAAACAGAATGTGGAATGTTCATGTCACAATTGTTAGCAGATTTGCCGGCTGTAATCCGTCCTGTTTCGCATTACGAACTTTCATTACAGCCACAGCTTTTGAGGTAATGTTCACAAGAAGCATATTATccgccgt                                                                >  2:10479/1‑132 (MQ=255)                                                                             taAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGTCAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               <  4:130508/138‑1 (MQ=255)                                                                             taAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               <  5:194836/138‑1 (MQ=255)                                                                               aaaCAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATg                                                            <  6:222442/139‑1 (MQ=255)                                                                                  cAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTATAGGTGACAAGCAGCATATTTTCAGGAGGTATGcc                                                          >  6:262148/1‑138 (MQ=255)                                                                                  cAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCt                                                         <  5:217141/139‑1 (MQ=255)                                                                                    gACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCt                                                         >  6:227537/1‑137 (MQ=255)                                                                                    gACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGCACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    >  8:186046/1‑94 (MQ=255)                                                                                    gACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGCACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    <  7:181710/94‑1 (MQ=255)                                                                                    gACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGCACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    >  8:181710/1‑94 (MQ=255)                                                                                    gACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGCACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    <  7:186046/94‑1 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCg                                                                 >  1:56622/1‑128 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCg                                                                 <  2:56622/128‑1 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               >  6:21553/1‑130 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               >  8:215408/1‑130 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               >  6:17438/1‑130 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               <  5:21553/130‑1 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               <  5:17438/130‑1 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGt                                                               <  7:215408/130‑1 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAg                                                      >  8:161990/1‑139 (MQ=255)                                                                                     aCTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTCACAAGCAGCATATTTTCAGTCGGTATGCCTAAg                                                      >  1:212600/1‑139 (MQ=255)                                                                                      cTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAgt                                                     <  2:278895/139‑1 (MQ=255)                                                                                                ttCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    >  4:242970/1‑82 (MQ=255)                                                                                                ttCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    <  3:242970/82‑1 (MQ=255)                                                                                                 tCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    >  7:246199/1‑81 (MQ=255)                                                                                                 tCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCttt                                                                                                    <  8:246199/81‑1 (MQ=255)                                                                                                   cTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTATCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGACTAAGTGTTCGATCAtgtg                                        >  5:248406/1‑139 (MQ=255)                                                                                                      tCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCAt                                           <  4:135303/133‑1 (MQ=255)                                                                                                      tCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCAt                                           >  3:135303/1‑133 (MQ=255)                                                                                                      tCACAATTCTTCTGACTTTTGCCGGCTCTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTgg                                       <  4:288213/137‑1 (MQ=255)                                                                                                                 ctGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACggg                             <  3:124625/136‑1 (MQ=255)                                                                                                                 ctGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCt                           <  1:101596/138‑1 (MQ=255)                                                                                                                    aCTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCTGt                         <  5:227537/137‑1 (MQ=255)                                                                                                                                         tACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACggg                             <  2:153323/112‑1 (MQ=255)                                                                                                                                         tACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACggg                             >  1:153323/1‑112 (MQ=255)                                                                                                                                         tACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCTGTCAGTACCCGTGTTTTTCCGCTcc  >  1:8981/1‑139 (MQ=255)                                                                                                                                         tACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCTGTCAGAACCCGTGTTTTTCCGCTcc  >  3:37384/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                           TTGTTCCCGCTGCCACTCCATTTTCAGCAGATGCTCACCGTTCCATTTAGCAGAATTAGAGTGATACAAGATCTTATAAAACAGACTGTGGAATGTTCCTGTCACAATTCTTCTGACTTTTGCCGGCTGTAATCCGTACTGGTTCGCCATACGCTCTTTCATTTCAGCGACAGCTTTTGTTGTAAAGGTGACAAGCAGCATATTTTCAGGCGGTATGCCTAAGTGTTCGATCATGTGGGCGGCACGGGCTGTCAGTACCCGTGTTTTTCCGCTCC  >  BS168_reference_genome/1255246‑1255520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑