Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,471,869	T→C	12.5%	K203R (AAA→AGA)	yvbG ←	putative integral membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,471,869	0	T	C	12.5%	74.1 / 5.3	32	K203R (AAA→AGA)	yvbG	putative integral membrane protein
Reads supporting (aligned to +/- strand):  ref base T (12/16);  new base C (2/2);  total (14/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 = 6.46e-01

GGCTGAGCAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGTGGAAAAAGAGAAAGGTCAG  >  BS168_reference_genome/3471738‑3472002                                                                                                                                    |                                                                                                                                      ggCTGAGCAGGTTGATTGTATGTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGcc                                                                                                                                 <  1:126755/138‑1 (MQ=255) ggCTGAGCAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGcc                                                                                                                                 <  5:170665/138‑1 (MQ=255)  gCTGAGCAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCg                                                                                                                               >  1:133699/1‑139 (MQ=255)      aGCAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCAc                                                                                                                           <  8:205958/139‑1 (MQ=255)        cAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCAcc                                                                                                                          >  1:139365/1‑138 (MQ=255)        cAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACAGGAAACATGCCTTACAGCCCGCCAcc                                                                                                                          >  3:162151/1‑138 (MQ=255)           gTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAatc                                                                                                                      >  4:51413/1‑139 (MQ=255)               attgaTTTTTTTTAGAGGGAGGAAGTCTTTCGGTTTGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAataatc                                                                                                                   <  2:200072/138‑1 (MQ=255)               attgaATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATcc                                                                                                                  <  5:75389/139‑1 (MQ=255)                 tgaATTTTATGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGATGCTTTTTGCATCATGAAGTTAAGACCGGAAACAAGCCTTTCAGCCCGGCACCAaacaaaccc                                                                                                                >  8:124680/1‑134 (MQ=255)                          tGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAAc                                                                                                        <  7:250103/138‑1 (MQ=255)                                tGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCt                                                                                                 <  3:51413/139‑1 (MQ=255)                                                  tgtgCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGt                                                                               <  4:92196/139‑1 (MQ=255)                                                  tgtgCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCACATCAGt                                                                               >  1:151051/1‑139 (MQ=255)                                                              aaaaGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCt                                                                                                 <  1:205069/109‑1 (MQ=255)                                                              aaaaGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCt                                                                                                 >  2:205069/1‑109 (MQ=255)                                                                     tGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTCTCAGCCCGGCAc                                                                                                                           >  5:258191/1‑76 (MQ=255)                                                                     tGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTCTCAGCCCGGCAc                                                                                                                           <  6:262370/76‑1 (MQ=255)                                                                     tGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTCTCAGCCCGGCAc                                                                                                                           <  6:258191/76‑1 (MQ=255)                                                                     tGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTCTCAGCCCGGCAc                                                                                                                           >  5:262370/1‑76 (MQ=255)                                                                                               ctTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGAt                                   >  1:148866/1‑138 (MQ=255)                                                                                               ctTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGAt                                   <  5:130141/138‑1 (MQ=255)                                                                                                       atcatGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGa                          <  3:120026/139‑1 (MQ=255)                                                                                                       atcatGAAGTTAAGACCGGAAACATGCCGTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGa                          <  4:165659/139‑1 (MQ=255)                                                                                                          atGAAGTTATGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAAt                        <  7:124680/138‑1 (MQ=255)                                                                                                             aaGTTAATCCCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACAGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGt                     <  4:162151/138‑1 (MQ=255)                                                                                                             aaGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGt                     <  2:139365/138‑1 (MQ=255)                                                                                                                         gAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAAc                                                                                                        <  3:232767/43‑1 (MQ=255)                                                                                                                         gAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAAc                                                                                                        >  4:232767/1‑43 (MQ=255)                                                                                                                         gAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGTGGAAAAAgaga          >  4:119763/1‑137 (MQ=255)                                                                                                                          aaaCATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGTGGAAAAAGAGAAAg       >  5:78182/1‑139 (MQ=255)                                                                                                                               tGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGTGGAAAAAGAGAAAGGTCAg  >  5:156586/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      GGCTGAGCAGGTTGATTGAATTTTTTGAGAGTGAGGAAGTCTTTCGGTATGTGCCGAAAACAAAAGAATGCAGCTCTCTAAAATAAGAGAGCTGCTTTTTGCATCATGAAGTTAAGACCGGAAACATGCCTTTCAGCCCGGCACCAATCATCCCGACCGCAACCACGGCTAAGATTAATCCCATCAGTCTGGTAATCACATTCATCTCGGTTTTTCCGAGTTTGCTGCTGATAAAAGCAGAATAGTGGAAAAAGAGAAAGGTCAG  >  BS168_reference_genome/3471738‑3472002

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