Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,148,040	G→A	9.7%	A21V (GCG→GTG)	cgeD ←	putative glycosyltransferase involved in maturation of the outermost layer of the spore

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,148,040	0	G	A	9.7%	122.6 / 3.7	41	A21V (GCG→GTG)	cgeD	putative glycosyltransferase involved in maturation of the outermost layer of the spore
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (2/2);  total (20/21)
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.90e-01

TGCACGAAACTGTTATGATATTGAATGCGCCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCAGAGCCTCTGTCATCATTTAAAAAG  >  BS168_reference_genome/2147902‑2148173                                                                                                                                           |                                                                                                                                      tGCACGAAACTGTTATGATATTGAATGCGCCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCg                                                                                                                                       >  3:293760/1‑139 (MQ=255)              tATGATATTGAATGCGCCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAaat                                                                                                                           >  8:278775/1‑138 (MQ=255)                      tGAATGCGCCGGTCATTGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATTATCGGGt                                                                                                                   <  4:293760/138‑1 (MQ=255)                            cgcCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTttgttg                                                                                                             >  1:324179/1‑138 (MQ=255)                            cgcCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTttgttg                                                                                                             >  7:288915/1‑138 (MQ=255)                                     tGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTa                                                                                                   <  5:190852/139‑1 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGc                                                                                                                               <  8:183133/104‑1 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGc                                                                                                                               >  7:183133/1‑104 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGt                                                                                                                   >  6:122251/1‑116 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGt                                                                                                                   <  5:122251/116‑1 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCACTTTTTGc                                                                                                                               <  6:75135/104‑1 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCACTTTTTGc                                                                                                                               >  5:75135/1‑104 (MQ=255)                                            aCTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCACTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAaataa                                                                                             >  2:72046/1‑138 (MQ=255)                                                gtgAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGa                                                                                         <  1:358527/138‑1 (MQ=255)                                                 tgAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGac                                                                                        <  3:89294/138‑1 (MQ=255)                                                 tgAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGac                                                                                        <  8:288915/138‑1 (MQ=255)                                                       aCCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtc                                                                                  >  2:349054/1‑138 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCTAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  1:349054/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGTGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  6:379585/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  3:23542/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  2:59067/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  7:278775/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  4:250752/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  2:324179/139‑1 (MQ=255)                                                         cGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCACTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  1:72046/139‑1 (MQ=255)                                                                        aTTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGcc                                                                >  3:20949/1‑139 (MQ=255)                                                                                                     cccAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGATCCg                                   >  8:265155/1‑139 (MQ=255)                                                                                                       cAAATATCTGGCTTTGGCTTTATCACGCTTTCCATGGATTTTTGCAAATAATCGGGTTTGTTGTAACTTTTTAAAATAATTCAAACTTTCTCACACATATATACCCCCGCCCGCGGTCAAAAATAAAAACGGCTCCgca                                 <  2:57723/139‑1 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCTGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCaga                        >  2:230218/1‑139 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCTGGTTTGTTGTAACTTGTGAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGCCTCCTCagcacgca                          >  6:120430/1‑137 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  1:224010/84‑1 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  2:144239/84‑1 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               >  2:224010/1‑84 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               >  1:144239/1‑84 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCagcaagca                          >  3:107228/1‑137 (MQ=255)                                                                                                                tGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTCTTGTAACTGGTTAAAATATTAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCaga                        >  6:405385/1‑139 (MQ=255)                                                                                                                      tGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCAGAGCctct                  >  1:319810/1‑139 (MQ=255)                                                                                                                           tATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               >  4:129623/1‑73 (MQ=255)                                                                                                                           tATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTtctct                                                                               <  3:129623/73‑1 (MQ=255)                                                                                                                                     tCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCAGAGCCTCTGTCATCATTTaaaaa   >  4:163795/1‑139 (MQ=255)                                                                                                                                     tCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCAGAGCCTCTGTCATCATTTaaa     >  2:175464/1‑137 (MQ=255)                                                                                                                                       caTCGCTTTTTGAAAATAATCGGGTTTGTTGTAACTTGTTAAAAAAATAGACACTTTCTCCCCCATCTAACCCGCCCCCCGCCGCAAAACATAAACACGGCTCCGCAGCAAGCAGAGCCTCTGTCATCATTTAAAAAg  <  6:85503/137‑1 (MQ=255)                                                                                                                                           |                                                                                                                                      TGCACGAAACTGTTATGATATTGAATGCGCCGGTCATGGAGATACTTGTGAATCACCGCAGTTGTTTCTGCATTCGAATGATCATCCATGATAAAAAGCTCCCAAAGATCATGCGTTTGCTGTATCACGCTTTCGATCGCTTTTTGCAAATAATCGGGTTTGTTGTAACTTGTTAAAATAATAGACACTTTCTCTCCCATCTATACCGCCTCCCGCCGTCAAACATAAAAACGGCTCCGCAGCAAGCAGAGCCTCTGTCATCATTTAAAAAG  >  BS168_reference_genome/2147902‑2148173

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