Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,096,417	C→T	8.9%	A50T (GCG→ACG)	ytaP ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,096,417	0	C	T	8.9%	134.1 / 3.8	45	A50T (GCG→ACG)	ytaP	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base T (2/2);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

GAAGGCGTTTTCAAATAATCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTTTCTTCTT  >  BS168_reference_genome/3096281‑3096553                                                                                                                                         |                                                                                                                                         gaagGCGTTTTCAAATAATCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCt                                                                                                                                        >  7:51047/1‑139 (MQ=255)        ttttCAAATAATCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATg                                                                                                                                 <  1:15705/139‑1 (MQ=255)        ttttCAAATAATCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATg                                                                                                                                 <  2:14678/139‑1 (MQ=255)                  atCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTaaa                                                                                                                       <  8:62134/139‑1 (MQ=255)                                    cTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTc                                                                                                      >  7:170927/1‑138 (MQ=255)                                      cccTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGGTTTCTTCGGCTTCACAAATTAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGGGCGGTTCatc                                                                                                   >  8:252401/3‑137 (MQ=255)                                      cACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCATGAGCGTTTCAAc                                                                                                   >  1:82503/1‑139 (MQ=255)                                           tttCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGata                                                                                               <  6:110933/138‑1 (MQ=255)                                               cccTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATAtttt                                                                                           <  5:143689/138‑1 (MQ=255)                                                            gCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAAAGATATTTTAATCtaattctta                                                                              >  6:234629/1‑137 (MQ=255)                                                              cTCCATGCGAATGCTGCAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCAAAAAATAATCGGGCGCCTTTTATTGTGAAATTAAATAGAGAAGGAGCGTTTCAACGATATTTTATTATCTTTGTTCaa                                                                            >  8:133930/1‑138 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGTCTTTTCATg                                                                                                                                 <  1:181060/85‑1 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGTCTTTTCATg                                                                                                                                 >  2:181060/1‑85 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGTCTTTTCATg                                                                                                                                 <  1:185254/85‑1 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGTCTTTTCATg                                                                                                                                 >  2:185254/1‑85 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCaa                                                                                                    >  1:258547/1‑114 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCaa                                                                                                    <  2:258547/114‑1 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCaa                                                                            >  5:54343/1‑138 (MQ=255)                                                              cTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGCGCGCCATTTCCTGTCCATTTAAAACGAGCCGGAGCGTTACAACGATATTTTCTACTAATTATTCaa                                                                            >  8:48442/1‑138 (MQ=255)                                                                 cATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCGTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCGTTGTTCAATACg                                                                        >  3:265096/1‑139 (MQ=255)                                                                    gagagtGCTGAAACAATACGGCTTGGCACGGTCCTTCTGTTTTCTTAGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCGATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACCTaa                                                                     <  5:88704/134‑1 (MQ=255)                                                                     cGAATTCTGAAACAATACTGCTGCGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGTCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTaa                                                                     <  7:133930/138‑1 (MQ=255)                                                                     cGAAAGCTGAATCAATTCGTCTGGGCACGGTCCTTCTGATTTCTTCTGCTTCACAAAATAATCGGGCTCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTaa                                                                     <  8:184886/138‑1 (MQ=255)                                                                        aTGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGATCGTTTCAACGAAATTTTCTTCTCTTTGTTCAATACGTAATGtc                                                                 >  5:267156/1‑139 (MQ=255)                                                                        aTGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGtc                                                                 >  7:188014/1‑139 (MQ=255)                                                                         tGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTGCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTGTATTTCTTCAAAACGTAATGta                                                                 >  7:114301/1‑137 (MQ=255)                                                                                    tACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGAtt                                                     >  7:181188/1‑139 (MQ=255)                                                                                    tACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTAATGTCCATTTAAATCGAGCAGGAGCGTTTAAACGATATTTTATTCTGTTTAATAAATACGTAATGTCTAAACACTGAtt                                                     >  6:14921/1‑139 (MQ=255)                                                                                    tACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCGTTTAATGTCGATTTAAATCGAGCAGGAGCGTTTCAACGATAATTTGTTCTCTTTCTGAAATACGTAATTTATAAacacatgt                                                      >  6:88704/1‑134 (MQ=255)                                                                                            ggCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCC‑TTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTcc                                              <  4:290154/137‑1 (MQ=255)                                                                                              cACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGAt                                           >  7:268891/1‑139 (MQ=255)                                                                                                       tCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCaca                                   >  6:62262/1‑136 (MQ=255)                                                                                                             ttcttcGGCTTCACAAAATAATCGGGCGCCTTTTAATGTCCATTTAAATAGAGCATGAGCCTTTCCACCATAATTTCTTCTCTATCTTAAAAAAGTAATCTCTACACAATGAATGGCATACCATCAGGCAAAacacaca                            >  2:4327/1‑136 (MQ=255)                                                                                                             ttcttcGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCaa                            >  4:219590/1‑139 (MQ=255)                                                                                                             ttcttcGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCaa                            >  8:85804/1‑139 (MQ=255)                                                                                                             ttcttcGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAAAAGCGCAGGAGCGTTTCAACGATATTTTCTTATCTATATTCAATACATAATGTCACAAAACTTATTGTACTCCCATCAGACAAATAACCa                             >  8:97594/1‑138 (MQ=255)                                                                                                               cttcGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAg                           >  3:101932/1‑138 (MQ=255)                                                                                                                     cTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGa                    <  8:170927/139‑1 (MQ=255)                                                                                                                        cacaAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTATTCAATACGTAATGTCTCAACACTGATTTGCCTCCGATCAGGCAAATCCCCAAGcaaaccaaac                 >  4:221116/1‑139 (MQ=255)                                                                                                                        cacaAAATAAGCGGGCGCCGTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTATCTTTCTTCAATACGTAATGTATCAACACTGATTGGCCTCCGATCCGGCAAATCACCAAGCAATCGcaa                  >  4:216899/1‑138 (MQ=255)                                                                                                                          ctataTAAGCGGTCGCCTTTTCATGTCGATTTATATCGAGCAGGAGCGGTTCAACGATATTTTCTTCTCTTTGTTCAATACGTAATGTCTCAACTCTGATTGGCCTCCGATCATGCAAATCACCAAGCAATCGAAACAg               <  1:4327/135‑1 (MQ=255)                                                                                                                          caAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAg               <  3:219590/139‑1 (MQ=255)                                                                                                                                  gCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTttc       <  8:188014/139‑1 (MQ=255)                                                                                                                                  gCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTttc       <  4:101932/139‑1 (MQ=255)                                                                                                                                  gCGGGCGCCTTTTCATCTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTttc       <  6:267156/139‑1 (MQ=255)                                                                                                                                       cgcCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTTTcttctt  >  7:41927/1‑139 (MQ=255)                                                                                                                                       cgcCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTTTcttctt  >  7:67745/1‑139 (MQ=255)                                                                                                                                       cgcCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTTTcttctt  >  1:208154/1‑139 (MQ=255)                                                                                                                                       cgcCTTTTCATGTCCATTTAAATCGACAAGGAGCGTTTAAACGATATTTTATTCTATTTCTTCAATACCTAATGTAACAAAACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAAAAGTTGGATTcttctt  >  6:158770/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                         GAAGGCGTTTTCAAATAATCCGCTCCCTCAATCAGCTCACTTTTTCCCCTATCATACTGGCCTCCATGCGAATGCTGAAACAATACGGCTGGGCACGGTCCTTCTGTTTTCTTCGGCTTCACAAAATAAGCGGGCGCCTTTTCATGTCCATTTAAATCGAGCAGGAGCGTTTCAACGATATTTTCTTCTCTTTCTTCAATACGTAATGTCTCAACACTGATTGGCCTCCGATCAGGCAAATCACCAAGCAATCGAAACAGTTGCTTTCTTCTT  >  BS168_reference_genome/3096281‑3096553

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