Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	779,310	A→G	10.6%	N826S (AAC→AGC)	yetA →	putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	779,310	0	A	G	10.6%	97.8 / 5.4	40	N826S (AAC→AGC)	yetA	putative enzyme
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base G (2/2);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

AAAAGCAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCAAAAAATGAA  >  BS168_reference_genome/779172‑779446                                                                                                                                           |                                                                                                                                         aaaaGCAAGATCCTCATCTCGCCGCCATAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGaa                                                                                                                                          <  6:275227/139‑1 (MQ=255) aaaaGCAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGaa                                                                                                                                          <  2:319612/139‑1 (MQ=255)   aaGCAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAaca                                                                                                                                        <  1:55656/139‑1 (MQ=255)   aaGCAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAaca                                                                                                                                        <  8:118797/139‑1 (MQ=255)      cAAGCTCCTCATCTCGCCGCCAAAGCGTGGATTCTTTTGCTAGAAGATCAGCTCAGCCTTCCGCCGCTTCCGATCACCCCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  2:144403/139‑1 (MQ=255)      cAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  8:328673/139‑1 (MQ=255)      cAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  5:22213/139‑1 (MQ=255)      cAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  6:156314/139‑1 (MQ=255)      cAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  6:303752/139‑1 (MQ=255)      cAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTg                                                                                                                                     <  6:329964/139‑1 (MQ=255)         gATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCTGAGCGCATCGAGACATGTACGAAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGtct                                                                                                                                  >  3:15683/1‑139 (MQ=255)         gATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAaca                                                                                                                                        <  3:263618/133‑1 (MQ=255)         gATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAaca                                                                                                                                        >  4:263618/1‑133 (MQ=255)                    cgccgcCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATg                                                                                                                           >  6:69536/1‑135 (MQ=255)                    cgccgcCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATg                                                                                                                           <  5:69536/135‑1 (MQ=255)                     gccgccAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAACCCGGAGCGCATCGAGACATGGACGCAGCTAGAGGCGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGtt                                                                                                                      >  3:101920/1‑139 (MQ=255)                     gccgccAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGtt                                                                                                                      >  7:199857/1‑139 (MQ=255)                        gccAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTaa                                                                                                                   <  4:49346/139‑1 (MQ=255)                         ccAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTaaa                                                                                                                  <  8:323191/139‑1 (MQ=255)                                             cTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGaa                                                                                               <  8:179853/138‑1 (MQ=255)                                                 aagaTAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGAt                                                                                          <  3:276697/139‑1 (MQ=255)                                                             gCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCt                                                                              >  8:289146/1‑139 (MQ=255)                                                                  aCGCCGCTTCCGATCATACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACTACGAACACTGTCTATCAATGGTGTTTAAATGTCATAGCGGCGCTGGACCTGATCGGAGATTCCCTTCCgg                                                                         >  1:290363/1‑139 (MQ=255)                                                                  aCGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCgg                                                                         >  2:132297/1‑139 (MQ=255)                                                                      cgcTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCaaa                                                                     <  3:95307/139‑1 (MQ=255)                                                                      cgcTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCaa                                                                      >  8:71290/1‑138 (MQ=255)                                                                        cTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGAt                                                                                          >  8:148280/1‑116 (MQ=255)                                                                        cTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGAt                                                                                          <  7:148280/116‑1 (MQ=255)                                                                        cTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCaaaa                                                                    >  3:247292/1‑138 (MQ=255)                                                                                   aCCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCg                                                        >  1:312606/1‑139 (MQ=255)                                                                                                gagaCATGGACGCAGCTAGAGGAGCTTCCATGGGTGACACCGAACACTGTCTCTCAATGGTGTTTAAATGTCAAAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGAt                                           <  2:148521/139‑1 (MQ=255)                                                                                                       ggACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTAATAGCGGAGCTGGAACTGATCGGAGATTCCCTTCCGGAAAAAAAAGAAACAACGGGCAAAAAAGGATAACAGcc                                    >  6:105978/1‑139 (MQ=255)                                                                                                        gACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCg                                   >  4:233542/1‑139 (MQ=255)                                                                                                                  gaggagCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCAt                          <  7:289146/138‑1 (MQ=255)                                                                                                                    ggagCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAgg                       <  3:233542/139‑1 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGATCACTGACTCTCTATGCTGTATAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAACAAACATCGTGGAATAAAGGATATCTGCCGGAATTTCATAGgag                     >  8:272076/1‑138 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGAGCACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGAtt                                                                                  <  5:289146/77‑1 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGAGCACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGAtt                                                                                  >  6:289146/1‑77 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGAGCACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGAtt                                                                                  <  5:284112/77‑1 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGAGCACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGAtt                                                                                  >  6:284112/1‑77 (MQ=255)                                                                                                                       gCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGgag                     >  4:352864/1‑138 (MQ=255)                                                                                                                        cTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGgaga                    >  4:229628/1‑138 (MQ=255)                                                                                                                         ttCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGaaa                  >  5:110549/1‑139 (MQ=255)                                                                                                                                    aCAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCaaaa        >  4:3826/1‑138 (MQ=255)                                                                                                                                    aCAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCaaaa        >  8:26334/1‑138 (MQ=255)                                                                                                                                         gAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCAAAAAATGaa  <  1:132297/139‑1 (MQ=255)                                                                                                                                           acacTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCAAAAAATg    >  5:216997/1‑135 (MQ=255)                                                                                                                                           acacTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCAAAAAATg    >  3:191146/1‑135 (MQ=255)                                                                                                                                           |                                                                                                                                         AAAAGCAAGATCCTCATCTCGCCGCCAAAGCGTGGAATCTTTTGCTAGAAGATAAGCTCAGCCATACGCCGCTTCCGATCAAACCGGAGCGCATCGAGACATGGACGCAGCTAGAGGAGCTTCCATGGGTGACAACGAACACTGTCTCTCAATGGTGTTTAAATGTCATAGCGGCGCTGGAACTGATCGGAGATTCCCTTCCGGCAAAAAAAGAAACATCGGGCAAAAAAGGATAACAGCCGGAATTTCATAGGAGAAATGGTGCAAAAAATGAA  >  BS168_reference_genome/779172‑779446

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