Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	4,019,488	A→C	17.1%	L41R (CTT→CGT)	yxiK ←	conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	4,019,488	0	A	C	17.1%	81.5 / 4.1	44	L41R (CTT→CGT)	yxiK	conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (28/7);  new base C (2/6);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 5.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATATGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTAT  >  BS168_reference_genome/4019354‑4019625                                                                                                                                       |                                                                                                                                           gTAATATGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTTGGCAAGccc                                                                                                                                        >  1:110543/1‑139 (MQ=255)  tAATATGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGccc                                                                                                                                        <  3:55345/138‑1 (MQ=255)     tatGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGcc                                                                                                                                    >  7:4727/1‑139 (MQ=255)     tatGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTACGCAAGCCCGGcc                                                                                                                                    >  1:116423/1‑139 (MQ=255)     taaggATTGGAAATGGAATTCCCTTTTGGCCAGCTACGTGACCTTGTAAGCATCACCTTCAAAATACTCGGACACGTCTTCTATAAATGAACCAGTCGTAATTAACAGAATATTGAGTCGAACGTAGGCACTCCCAGcc                                                                                                                                    >  4:195945/4‑139 (MQ=255)         gATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCTCCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTAACTGTATAAAGAATAGATCGTAGGCAACCCCTGCCCaa                                                                                                                                 >  8:28391/1‑138 (MQ=255)         gATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCaaa                                                                                                                                <  7:344767/139‑1 (MQ=255)         gATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCaaa                                                                                                                                <  3:61866/139‑1 (MQ=255)          aTTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTAAAAATGAACCAGTCGTTATTAACTGTATCATGATTCGATCGTAGGCAAGCCCAGCCCaaa                                                                                                                                >  5:372116/1‑138 (MQ=255)            tGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTTTTACCGGTATCATGATTCGCTCCTAGGCAAGCCCGGCCCAAAAAc                                                                                                                             >  7:202901/1‑139 (MQ=255)              gAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGa                                                                                                                           >  7:156536/1‑139 (MQ=255)                                      cTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGt                                                                                                  >  1:60032/1‑139 (MQ=255)                                          ggcACCTTTTACCCATCACATTAAAAATACTCGGCCCGGTCTTATAAAAATCACCCAGGCTTTATTACCTTTATCTTTTTTCTATCGTAGCCAAGCCCGCCCAAAAAACCAAAAC‑ACGACTCTGCATCTTCAAGTTCg                                                                                               <  2:60032/136‑1 (MQ=255)                                                 tGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCGGTATGATGATTCGATAGTAGGCAAACCCTGCCCAAAAACGAAAAA‑ACCACAATGCATATTCAAGTTAGCAAGCCa                                                                                        >  4:71232/1‑138 (MQ=255)                                                 tGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTAACTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCa                                                                                        >  4:51639/1‑138 (MQ=255)                                                 tGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTGATTAACTGTATAATGATTCGATCGTAGGCAAGCCAGGCCCAAAAACGATAAA‑AAGACTCTGCATCTTCAAGTTAGCATGCCa                                                                                        >  5:356716/1‑138 (MQ=255)                                                   tAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTTCCGGTTAAATGTTTAGATCGTAGGCACGCCCTGCCCCAAAACGATAAA‑ACGACTCTGGCTCTTACAGGTCGCCTGGCCTaa                                                                                     >  2:45804/1‑139 (MQ=255)                                                         tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGcc                                                                                                                                         <  3:246276/82‑1 (MQ=255)                                                         tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGcc                                                                                                                                         >  4:246276/1‑82 (MQ=255)                                                         tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCCAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCCTGCCATAACGGttt                                                                               >  3:249275/1‑139 (MQ=255)                                                         tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGttt                                                                               >  2:66169/1‑139 (MQ=255)                                                         tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACAAGTATAATGATTCGATCGTAGGCAAGCCCTGCCCCAAAACGATAAA‑ACGACTAAGCATCTTCAAGTTAGCATGACATAAAGGttt                                                                               >  3:12832/1‑139 (MQ=255)                                                            cattaaaaaTACTCGTACAGGTCTTATAAAAATCAACCAGGCGTTATAACCTGTATCATGATTCGATCGTAGGCAAGCCCGGCCAAAAAACCATAAC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGt                                                                            <  7:223321/134‑1 (MQ=255)                                                             cataaaaaTCCTCGGAATGTTATTATATATATTAACCGGTCTTTATTTCTTGTTTGATTATTGTTTCTTAGTCTTGCCCTGCCAAAAAAC‑ATAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGt                                                                            <  8:4727/134‑1 (MQ=255)                                                                               ggTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCAt                                                         >  3:74652/1‑139 (MQ=255)                                                                               ggTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCAt                                                         >  7:373284/1‑139 (MQ=255)                                                                                   ttctATAAATGAACCAGTCGTTATTAACTTTATAATGATTCGATAGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGTTCAGGGTTAAGCTAGATCATGTc                                                      >  4:256701/1‑138 (MQ=255)                                                                                        taAATCAACCTGTCGTTAATACCTGTATAATTATGCTATAGTACCAACGCCCTGCAAAAAACGGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAAtt                                                <  3:193699/139‑1 (MQ=255)                                                                                          aatcacccAGTCGTTATTACATGTATCATCATTCGATCGTCGCCAACCCCTCCCAAAAACCCAAAAC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGc                                              <  3:51639/133‑1 (MQ=255)                                                                                           aTGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTAGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGc                                              >  6:82888/1‑138 (MQ=255)                                                                                           aTGAACCAGTCGTTATTACCTGTATCATGATTCGAAACTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTTCATCTTCAAGTTCGCATGCCATAACGGTTAGGTCAGGCTTAAGCTCGATCATGTCCAAAATgg                                              >  3:302744/1‑137 (MQ=255)                                                                                           aTGAACCAGTCGTTATTACCAGTATCATGCTTAGATCGTAGGCAAGCCCAGCCCAAAAAAGCTAAA‑ACGACTCAGCATCTTCAAGTTCGCAAGCCATTAACGTTTGGTCACGCTTAATCTCGATAATGACCAAATAGcc                                             >  8:412730/1‑139 (MQ=255)                                                                                               ccccgTCGTTATTACCTTTACAATTATTCGAACTTCCGCACCCCCGCCCCAAAAACTAAAAC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATa                                         <  5:384622/135‑1 (MQ=255)                                                                                               cccAGCCTTTATAACCTGTCACATGATTCCATCGGCCCCACCCCCCCCCAAAAACCCAAACC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATa                                         <  4:336617/138‑1 (MQ=255)                                                                                               aCCAGTCTTTATTACATGAATAATGATTCTATCGTCGCCACCCCCTCCCCAAAAACTAAAAC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATa                                         <  8:373284/139‑1 (MQ=255)                                                                                                cccgTCGTTTTTTCCGTTATCATGATTCGTGCGTCGGCCCCCCCTGCCCAAAAACGAAAAC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATAc                                        <  8:202901/136‑1 (MQ=255)                                                                                                ccAGTCGTGATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCCTAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATAc                                        >  4:391125/1‑139 (MQ=255)                                                                                                         tttaCCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTAGATAATGTACAAATTGCCAATACTGCTta                                  >  3:336617/2‑136 (MQ=255)                                                                                                                 tATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTCCTTATAGAGGTTAtt                       >  3:402973/1‑139 (MQ=255)                                                                                                                  aTCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAttt                      >  6:396101/1‑139 (MQ=255)                                                                                                                  aTCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAttt                      >  8:257251/1‑139 (MQ=255)                                                                                                                        aTTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATa                 >  7:71732/1‑138 (MQ=255)                                                                                                                                ttAGCCACGCCCTCCAAAAAAACCATACC‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCtt          <  3:71232/136‑1 (MQ=255)                                                                                                                                  aGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTccgc      >  1:192850/1‑137 (MQ=255)                                                                                                                                       aGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTAt  >  3:382143/1‑138 (MQ=255)                                                                                                                                       |                                                                                                                                           GTAATATGGATTGGAAATGGAATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAA‑ACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTAT  >  BS168_reference_genome/4019354‑4019625

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