Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,695,698	C→T	6.7%	D112D (GAC→GAT)	ribZB →	5‑amino‑6‑ribitylamino‑2,4(1H, 3H)‑pyrimidinedione 5'‑phosphate phosphatase (promiscuous activity)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,695,698	0	C	T	6.7%	190.9 / 3.4	60	D112D (GAC→GAT)	ribZB	5‑amino‑6‑ribitylamino‑2,4(1H, 3H)‑pyrimidinedione 5'‑phosphate phosphatase (promiscuous activity)
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/2);  total (30/30)
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.18e-01

CACGACACGGGTTACCGTTTGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAACCAATCCCTCATA  >  BS168_reference_genome/3695564‑3695825                                                                                                                                       |                                                                                                                                cACGACACGGGTTACCGTTTGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCg                                                                                                                             <  1:255355/139‑1 (MQ=255) cACGACACGGGTTACCGTTTGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCg                                                                                                                             <  7:88287/139‑1 (MQ=255) cACGACACGGGTTACCGTTTGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCg                                                                                                                             <  3:47500/139‑1 (MQ=255)                  tttGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCAcgcg                                                                                                            <  1:244098/139‑1 (MQ=255)                    tGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCAcg                                                                                                              <  5:132158/135‑1 (MQ=255)                    tGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCAcg                                                                                                              >  6:132158/1‑135 (MQ=255)                         aGCCGGACATTAATGGTTCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  7:254565/139‑1 (MQ=255)                         aGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGTTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  4:180403/139‑1 (MQ=255)                         aGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  1:24582/139‑1 (MQ=255)                         aGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  1:4567/139‑1 (MQ=255)                         aGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  7:265274/139‑1 (MQ=255)                         aGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  7:130642/139‑1 (MQ=255)                          gacGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGa                                                                                                    >  7:137856/3‑139 (MQ=255)                           ccGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGAc                                                                                                   <  6:179085/139‑1 (MQ=255)                                     aTGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTaaa                                                                                           >  3:72844/1‑137 (MQ=255)                                                aGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCg                                                                                                                             <  6:44927/92‑1 (MQ=255)                                                aGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCg                                                                                                                             >  5:44927/1‑92 (MQ=255)                                                            gATTGCTGATTACTTACTTTCAAAAACCATTTACTTTGAGGTATATCCAGATGCCCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAACCGAGCAGACg                                                                  >  4:135713/1‑139 (MQ=255)                                                             aTTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCaaa                                                                                  <  2:40599/122‑1 (MQ=255)                                                             aTTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCaaa                                                                                  >  1:40599/1‑122 (MQ=255)                                                               tGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGAt                                                               >  5:161356/1‑139 (MQ=255)                                                                gCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGAt                                                               <  2:14373/138‑1 (MQ=255)                                                                           acttTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGAc                                                                                                   <  6:118144/91‑1 (MQ=255)                                                                           acttTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGAc                                                                                                   >  5:118144/1‑91 (MQ=255)                                                                           acttTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGc                                                     >  6:112219/1‑137 (MQ=255)                                                                              ttCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAg                                                  >  7:281953/1‑137 (MQ=255)                                                                                      aCATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGAc                                        <  3:181372/139‑1 (MQ=255)                                                                                      aCATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGAc                                        <  7:183807/139‑1 (MQ=255)                                                                                                 gAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCaaa                              <  8:281953/138‑1 (MQ=255)                                                                                                     tatataCAGATGACCATCTGCTTTCTCCTTTTGCCGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGtt                         <  2:254565/139‑1 (MQ=255)                                                                                                       tataCAGATGACCATCTGCTTTCTCCTTTTGATGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGa                                                                <  6:52559/98‑1 (MQ=255)                                                                                                       tataCAGATGACCATCTGCTTTCTCCTTTTGATGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGa                                                                >  5:52559/1‑98 (MQ=255)                                                                                                       tataCAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGttt                        <  5:112219/138‑1 (MQ=255)                                                                                                       tataCAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGttt                        <  8:137856/138‑1 (MQ=255)                                                                                                            aGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCAt                   >  3:251867/1‑138 (MQ=255)                                                                                                            aGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGAGGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCAt                   >  3:216808/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGATGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAgg                                                 >  2:56562/1‑106 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGATGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAgg                                                 <  1:56562/106‑1 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGCGCGTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACGTGTGGCAAGGAGCCAGGACTCAATTCAAACAGTTTGGCATTa                 >  7:149310/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     <  1:84225/54‑1 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTg                                                                                                     >  2:84225/1‑54 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGATATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTa                 >  2:76949/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGtt                         >  6:915/1‑130 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGtt                         <  5:915/130‑1 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTa                 >  7:140286/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTa                 >  2:110585/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTa                 >  4:130672/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGAATTa                 >  5:47672/1‑138 (MQ=255)                                                                                                              aTGACCATCTGCTTTCTCCTTTTGACGGAGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGAAAGGAGCCATGACTAAATTCAAACAGTTTGGCATTa                 >  8:232771/1‑138 (MQ=255)                                                                                                                 aCCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTaaa                                                                                           >  6:238418/1‑61 (MQ=255)                                                                                                                 aCCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTaaa                                                                                           <  5:238418/61‑1 (MQ=255)                                                                                                                 aCCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACTAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGTCATTAAAc              >  6:197770/1‑138 (MQ=255)                                                                                                                 aCCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAAc              >  4:280322/1‑138 (MQ=255)                                                                                                                   cATCTGCTTTCTCCTTTTGACGGCGAGGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAACCGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAACCa            >  4:154356/1‑138 (MQ=255)                                                                                                                    aTCTGCTTTCTCCTTTTGACGGCGAATCTAAGCTGCACTCGGAGCTTGACATTCTAAAGAGCGCAAAGCCTAACGAGCAGACTGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAACCAAt          <  7:232771/139‑1 (MQ=255)                                                                                                                        gCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGttt                        >  3:82019/1‑121 (MQ=255)                                                                                                                        gCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGttt                        <  4:82019/121‑1 (MQ=255)                                                                                                                            tctcCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCACACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAACCAATCCCTCata  >  4:104833/1‑139 (MQ=255)                                                                                                                            tctcCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCAt                   <  3:277343/122‑1 (MQ=255)                                                                                                                            tctcCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCAt                   >  4:277343/1‑122 (MQ=255)                                                                                                                                       |                                                                                                                                CACGACACGGGTTACCGTTTGATCAGCCGGACATTAATGGATCAGGAAGCCGGTAAAGCGATTGCTGATTACTTACTTTCAAAAAACATTTACTTTGAGGTATATACAGATGACCATCTGCTTTCTCCTTTTGACGGCGAAGCTAAGCTGCACGCGGAGCTTGACATTCTAAAGAGCGCAAACCCAAACGAGCAGACGGATGACTTGTGGCAAGGAGCCATGACTCAATTCAAACAGTTTGGCATTAAACCAATCCCTCATA  >  BS168_reference_genome/3695564‑3695825

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