Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,913,926	C→A	16.7%	intergenic (‑413/‑83)	rodA ← / → ywcE	glycosyltransferase involved in extension of the lateral walls of the cell/holin‑like protein required for proper spore morphogenesis and germination

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,913,926	0	C	A	16.7%	100.7 / 5.8	43	intergenic (‑413/‑83)	rodA/ywcE	glycosyltransferase involved in extension of the lateral walls of the cell/holin‑like protein required for proper spore morphogenesis and germination
Reads supporting (aligned to +/- strand):  ref base C (24/11);  new base A (2/6);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 4.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCAGAGCCCCATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCCTCTCT  >  BS168_reference_genome/3913795‑3914060                                                                                                                                     |                                                                                                                                       aTTCAGAGCCCCATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAAc                                                                                                                                 <  8:153741/139‑1 (MQ=255) aTTCAGAGCCCCATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAAc                                                                                                                                 <  4:105159/139‑1 (MQ=255)     agagCCCCATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAACACAAt                                                                                                                              >  2:170445/1‑138 (MQ=255)          cccATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACaaa                                                                                                                        <  2:104722/139‑1 (MQ=255)                   aTACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGacac                                                                                                                >  2:220039/1‑138 (MQ=255)                      cATTCTCCTAGCACATTAAACCAAAAAAAACGACGATCGATCCAGCCCATTTCTGCAAAAACACGTATCAATTGCTGGAAAACATGGT‑ATAATAATTGAAACAGAATGTTCAAAAAACAATCACAAAACCGACACAACt                                                                                                            <  6:51610/139‑1 (MQ=255)                       aTTCTCCTAGCACATTAAACTAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTc                                                                                                           <  3:139844/139‑1 (MQ=255)                       aTTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGCAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACACACACAACTc                                                                                                           <  4:105212/139‑1 (MQ=255)                         tctcCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCGATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAg                                                                                                         >  3:164383/1‑139 (MQ=255)                         tctcCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACAAGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAAACACAAAACAGACACAACTCAg                                                                                                         >  8:39772/1‑139 (MQ=255)                           tcCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGtt                                                                                                       >  5:196257/1‑139 (MQ=255)                           tcCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTAAGtt                                                                                                       >  6:135769/1‑139 (MQ=255)                             cTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTc                                                                                                      >  8:190923/1‑138 (MQ=255)                                     ttAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTAttt                                                                                              >  6:133367/1‑138 (MQ=255)                                      tAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATAGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTAATTCa                                                                                            >  4:311026/1‑139 (MQ=255)                                        acccaaaaaaaaCGCCGAAAGATCCAGTCCATTTCTCCAAAAACACGAATCAATTTCTGTAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCAt                                                                                           <  7:20405/134‑1 (MQ=255)                                              aaaaTAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGt                                                                                     >  4:168105/1‑138 (MQ=255)                                                                  gTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGa                                                                  >  3:112188/1‑137 (MQ=255)                                                                  gTCCATTTCTCAAAAAACACGTATAAATTTCTGTAAAAAATTTT‑ATAATAATTGAAAAAGAATGTTAAAAAAACAATCACAAAACCCCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGaaa                                                                <  3:168105/139‑1 (MQ=255)                                                                  gTCCATTTCGGCAAAAACACGTATCAATTTCTGGAACACATGGT‑AAAATAATTGTAACAGAATGTTAAAAAAACAATCACAAAACAGCCACAACTCAGTTCCTCTATTACATTTTTGTAAACGATATCACAAAATGaaa                                                                >  6:295406/1‑139 (MQ=255)                                                                  gTCCATTCCTCAAAAAACACGAACCAATTTATCCAAAAC‑TGGTAATAAAATTTGAAACAGAATGATCAAAAAACAATCACAAACCCCCCACAACTCCGTTCCGCTATTTCATTTTTGTAAACGTTATCACAAAATGaaa                                                                <  7:273879/139‑1 (MQ=255)                                                                  gTCAATTTCTCCAAAAACACGTATCAATTTCTGCAAAACATTGA‑ATAATAATTTAAAAATAATGTAAAAAAAACAAACACAAAACACCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGaaa                                                                <  3:198435/139‑1 (MQ=255)                                                                   tCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATTAAAg                                                               >  5:187021/1‑139 (MQ=255)                                                                   tCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAg                                                               >  1:232674/1‑139 (MQ=255)                                                                   tCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAg                                                               >  6:55131/1‑139 (MQ=255)                                                                   tCCATTTCTGCAAAAACACGTATCAAGTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGt                                                                               >  8:167391/1‑123 (MQ=255)                                                                   tCCATTTCTGCAAAAACACGTATCAAGTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGt                                                                               <  7:167391/123‑1 (MQ=255)                                                                       tttCTGCAAAAACACGTATCAATTTCTGTAAAACATGTT‑AAAATATTTGTAACAGAATGTTCAAAAAACAATCACAAAACACACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGgaga                                                           <  6:121230/139‑1 (MQ=255)                                                                        ttCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGAc                                                          >  5:121230/1‑139 (MQ=255)                                                                        ttCTGAAAAAAAACGTATCAATTTCTGTAAAACATGTT‑ATAATAATTGAAACAGAATGTAAAAAAAACAAACAAAAAACACCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGgaga                                                           <  2:232674/138‑1 (MQ=255)                                                                         tctcCAAACACACGTATACATTTCTCTAAAACATGCT‑ATAAAAATTGAAACAGAATTTTCAAAAAACACACACAAAACACCCACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGAc                                                          <  8:73505/138‑1 (MQ=255)                                                                             cAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGAAGGGTa                                                     >  1:277657/1‑139 (MQ=255)                                                                                            aTTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAACCAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGtttt                                      >  7:203169/1‑139 (MQ=255)                                                                                                cTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGc                                  >  4:276337/1‑139 (MQ=255)                                                                                                  tcaaaaaaTGTT‑AAAATATTTGTACAATATTGTAAAAAAAAAAACAAAAAACCACACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCAAAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGct                                 <  7:190923/136‑1 (MQ=255)                                                                                                                        gTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATTAAAGGAGACGGGTATGATGGACATGTTTTTTGATTAATTATTGGTTGCCAGCGCaa           >  6:152589/1‑139 (MQ=255)                                                                                                                        gTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGAGATGTTTTTTGCTTACTTATTGGTTGCCAGCGCaa           >  8:77222/1‑139 (MQ=255)                                                                                                                        gTAACAGAATGTTAAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATCAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGATTGCCAGCGCaa           >  6:198189/1‑139 (MQ=255)                                                                                                                            aatattGTTAAAAAAAAAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGcc       <  4:164383/134‑1 (MQ=255)                                                                                                                                ttttacAAAAAACAATCAAAAAACACACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCctct    <  3:311026/133‑1 (MQ=255)                                                                                                                                aTGTTAAAAAAACAAAAAAAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCctct    <  4:112188/138‑1 (MQ=255)                                                                                                                                 tGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCctctct  >  6:124081/1‑139 (MQ=255)                                                                                                                                 tGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCctctct  >  3:95518/1‑139 (MQ=255)                                                                                                                                 tGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCctctct  >  3:158977/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                       ATTCAGAGCCCCATGCCCATACATTCTCCTAGCACATTAAACGAAAAAATAAGACGATAGATCCAGTCCATTTCTGCAAAAACACGTATCAATTTCTGGAAAACATGGT‑ATAATAATTGTAACAGAATGTTCAAAAAACAATCACAAAACAGACACAACTCAGTTCCTCTATTTCATTTTTGTAAACGTTATCACAAAATGAAAGGAGACGGGTATGATGGACATGTTTTTTGCTTACTTATTGGTTGCCAGCGCAACGCCTCTCT  >  BS168_reference_genome/3913795‑3914060

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