Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,236,265	T→A	16.1%	intergenic (+46/+157)	yuxK → / ← yufK	putative thiol‑disulfide oxidoreductase/putative integral membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,236,265	0	T	A	16.1%	101.0 / 5.1	41	intergenic (+46/+157)	yuxK/yufK	putative thiol‑disulfide oxidoreductase/putative integral membrane protein
Reads supporting (aligned to +/- strand):  ref base T (5/29);  new base A (5/2);  total (10/31)
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

TCATTGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTCATTGCGTACAGAGCACG  >  BS168_reference_genome/3236131‑3236383                                                                                                                                       |                                                                                                                       tCATTGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCt                                                                                                                     >  7:198371/1‑138 (MQ=255) tCATTGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATACt                                                                                                                     >  8:259790/1‑138 (MQ=255) tCATTGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAAACCt                                                                                                                     >  8:264794/1‑138 (MQ=255)    ttGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTaaaa                                                                                                                 >  6:283591/1‑139 (MQ=255)      gCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCACCCCTTCAATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGATAGGAAGATAGCCATCACAGATTGTTGACAAAAAAATAAAACg                                                                                                               >  6:273550/1‑139 (MQ=255)         aaaaaCCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCt                                                                                                            <  8:37436/139‑1 (MQ=255)         aaaaaCCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCt                                                                                                            <  4:65747/139‑1 (MQ=255)          aaaaCCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCtt                                                                                                           <  5:41542/139‑1 (MQ=255)                 cTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTa                                                                                                    >  8:351697/1‑139 (MQ=255)                        tGGTTTGGTAAAAAAAATGATTGTATGATACCCTCCCCTTCCAATAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCTTCAAAGATTGTTGACAAAATCCTAAATCGGCTTCGTTTTAGGAtttt                                                                                             <  2:228350/139‑1 (MQ=255)                        tGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  5:289909/139‑1 (MQ=255)                        tGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  8:198371/139‑1 (MQ=255)                        tGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  7:159729/139‑1 (MQ=255)                        tGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCATTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGATAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  2:285852/139‑1 (MQ=255)                        tGGTTTGGCAAAAAAAATCAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  3:348198/139‑1 (MQ=255)                         ggTTTGTCAAAAAAAATGAATGTATGTACCCCTCCCATACAATTAAAAAAATATTTTTACCTTAAAAAAAAAATGGCTAGGAAGTTAGCCTTAACTGATTTTTGAAAAAATCCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  2:155919/138‑1 (MQ=255)                          gtatGGCAAAAAAAATGAAAGTATGCTCCCCACCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAAACCTAAAACGGCTTCGTTTTAGGAtttt                                                                                             <  1:172074/134‑1 (MQ=255)                              ggCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCAt                                                                                        <  1:216085/138‑1 (MQ=255)                              ggCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCAt                                                                                        <  3:323931/138‑1 (MQ=255)                              ggCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCAt                                                                                        <  3:157434/138‑1 (MQ=255)                              ggCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCAt                                                                                        <  5:90227/138‑1 (MQ=255)                                  aaaaaaaTGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTAtt                                                                                   <  3:200075/139‑1 (MQ=255)                                      aaaTGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGc                                                                               <  3:66817/139‑1 (MQ=255)                                               tATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGaa                                                                      <  3:163826/139‑1 (MQ=255)                                                 tGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAAAACTAAAAAGGCGGGGTTTTATGATTTTTTGATAGATTAAGAGTGATTTaaa                                                                     >  1:222375/1‑138 (MQ=255)                                                        ctcccctTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAAAACTAAAAAGGCGTCGTGTTAGGTTTTTTTTATATAATAAGAGGGATTGAAAAAttttt                                                              >  7:170321/1‑137 (MQ=255)                                                                ccATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCtt                                                                                                           >  8:367608/1‑85 (MQ=255)                                                                ccATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCtt                                                                                                           <  7:367608/85‑1 (MQ=255)                                                                  aTTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGaaa                                                   <  7:229166/139‑1 (MQ=255)                                                                                       taaataaaAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAAAACAAAAAACGCGTCGTTTTTGGATTTTGTTATGTATTAAGGGGGATTGTTAAAATTTTTAGTTTTGAAAAGGCGGAGACCGGAGaa                                  >  5:249894/1‑135 (MQ=255)                                                                                         aataaaAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAgc                            <  5:283591/139‑1 (MQ=255)                                                                                         aataaaAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAgc                            <  7:351697/139‑1 (MQ=255)                                                                                           taaaAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAgccg                          <  7:264794/139‑1 (MQ=255)                                                                                           taaaAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAgccg                          <  7:259790/139‑1 (MQ=255)                                                                                                aaTGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGtt                     <  2:295833/139‑1 (MQ=255)                                                                                                aaTGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGtt                     <  4:264494/139‑1 (MQ=255)                                                                                                aaTGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGtt                     <  6:346998/139‑1 (MQ=255)                                                                                                  tGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTc                   <  2:222375/139‑1 (MQ=255)                                                                                                     cTACGAAGCAAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTCAtt                <  5:32409/139‑1 (MQ=255)                                                                                                                gCCATCAAATATTGTTGACAAAATCATAAAACGGCTTCGTTTTAGTATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTCATTGCGTACagag      <  4:22695/138‑1 (MQ=255)                                                                                                                   aTCACAGATTGTTGACAAAAACCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGGGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTCATTGCGTACAGAGCACg  <  3:13267/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                       TCATTGCCAAAAACCGCTATAAGTGGTTTGGCAAAAAAAATGAATGTATGCTCCCCTCCCCTTCCATTAAAAAAAGATTTTTACCGTAAATAAAAAATGGCTAGGAAGCTAGCCATCACAGATTGTTGACAAAATCCTAAAACGGCTTCGTTTTAGGATTTTGTCATCTATTCAGCGTGATTGAAAACTTTTGAAGTTGTGAAACGGCGGAAACCTGAGAATCTAGCCGCCGTTTCATTGCGTACAGAGCACG  >  BS168_reference_genome/3236131‑3236383

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