Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,019,150	C→T	8.2%	Q103Q (CAG→CAA)	ytfI ←	conserved sporulation‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,019,150	0	C	T	8.2%	144.7 / 3.2	49	Q103Q (CAG→CAA)	ytfI	conserved sporulation‑related protein
Reads supporting (aligned to +/- strand):  ref base C (26/18);  new base T (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ATACAGCACCGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGTCTTCTTTGTTCACCTTA  >  BS168_reference_genome/3019026‑3019288                                                                                                                             |                                                                                                                                           aTACAGCACCGGTTAACACACCTGTCACTGCCGCGTCCATAATTCCCATTACTGTAACCCAATCTAGGTTTGTGATATGTAAATTCGCGAGAAAAGACCCGCTGATTTTCCGCAGTCTGGTGAAATGATGCAAAAGCAt                                                                                                                              >  7:163292/1‑139 (MQ=255)       cACCGGTTACCCAACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGGGATATGTAAATGCGCGAGAAACGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        <  5:205688/139‑1 (MQ=255)       cACCGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTATATGCGCGAGAAAAGACGCGCAGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        <  5:194374/139‑1 (MQ=255)       cACCGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        >  7:208061/1‑139 (MQ=255)       cACCGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        <  5:186096/139‑1 (MQ=255)       cACCGGTTAACACACCTGTCACTGCCGCGGCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        >  4:19248/1‑139 (MQ=255)       cACCGGCTAACACACCTGTCACTGCCGCGTCCTGAATTCCTATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGCGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        <  7:266664/139‑1 (MQ=255)          cGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAtttt                                                                                                                     >  3:225518/1‑139 (MQ=255)          cGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAtttt                                                                                                                     >  6:299928/1‑139 (MQ=255)            gTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGCGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTAGAtttttt                                                                                                                   >  6:47545/1‑137 (MQ=255)            gTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGCAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAATCATTTCGATTttat                                                                                                                   >  7:181678/1‑139 (MQ=255)            gTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTttat                                                                                                                   >  8:300895/1‑139 (MQ=255)            gTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTttat                                                                                                                   >  7:89663/1‑139 (MQ=255)            gTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTttat                                                                                                                   >  2:17363/1‑139 (MQ=255)                      cTGTCACTGCCGCGTGCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTtg                                                                                                         <  1:255927/139‑1 (MQ=255)                      cTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTGTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCCGTCTGGAGACCTGATGCAAAAGCCTTTCGATTTTAGTTATGCTTtg                                                                                                         >  2:93040/1‑139 (MQ=255)                      cTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGCCCTGATGCAAAAGCATTTCGATTTTATTAATGCTTtg                                                                                                         >  1:111406/1‑139 (MQ=255)                      cTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTtg                                                                                                         >  5:56746/1‑139 (MQ=255)                       tGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATTCTTtgc                                                                                                        >  4:165774/1‑139 (MQ=255)                            cTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTAGCTGca                                                                                                    >  3:51696/1‑138 (MQ=255)                               ccGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGcaca                                                                                                  <  1:191641/137‑1 (MQ=255)                                  cgTCCTGAAATGCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGTAAAGACGCGCTGATTTTCCGCTGTCTGGTGACCTGTTTCAAAAGCATTTGGATTTTATTAATGCTTTGCTGCACATCGtt                                                                                             <  7:122914/139‑1 (MQ=255)                                       tGAATTCCGATTACTGTTACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGTTTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTtgtg                                                                                        <  8:163292/139‑1 (MQ=255)                                              cGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAAtttt                                                                                 <  3:78006/139‑1 (MQ=255)                                                    cTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTtttct                                                                            <  8:47996/138‑1 (MQ=255)                                                         aCCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCgc                                                                      >  5:47457/1‑139 (MQ=255)                                                         aCCCATTCTAGTTTTGAGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTAGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCgc                                                                      >  2:95804/1‑139 (MQ=255)                                                          cccATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATg                                                                                                              <  1:48606/98‑1 (MQ=255)                                                          cccATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATg                                                                                                              >  2:48606/1‑98 (MQ=255)                                                             aTTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCgctgct                                                                  >  3:281388/1‑139 (MQ=255)                                                             aTTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCgctgct                                                                  >  7:39349/1‑139 (MQ=255)                                                              ttCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGCCCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCtt                                                                 >  4:102274/1‑139 (MQ=255)                                                              ttCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCtt                                                                 >  6:201680/1‑139 (MQ=255)                                                               tCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACTTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCCGCttt                                                                >  7:25506/1‑139 (MQ=255)                                                               tCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACTTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCCGCttt                                                                >  7:30106/1‑139 (MQ=255)                                                               tCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCttt                                                                >  8:255737/1‑139 (MQ=255)                                                               tCTAGATTTGTCATATGAAAATGCCCGAGAAAAGACGCGCTGATTTTCCGCACTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGATTTGCTGCACATCGTTTTGTCAAATTATTTTCTTTCCgctgcata                                                                >  2:112577/1‑136 (MQ=255)                                                                         tgATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        >  8:41944/1‑73 (MQ=255)                                                                         tgATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGAt                                                                                                                        <  7:41944/73‑1 (MQ=255)                                                                                         gagaAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTtttgt                                      <  8:39349/139‑1 (MQ=255)                                                                                         gagaAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTtttgt                                      <  4:281388/139‑1 (MQ=255)                                                                                            aaaaGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTtttgtttg                                   <  8:89663/139‑1 (MQ=255)                                                                                              aaGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTACACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGat                                 <  7:300895/139‑1 (MQ=255)                                                                                                   gcgcTGATTTTCCGCAGTCTGGGGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCa                            <  4:80625/139‑1 (MQ=255)                                                                                                              ccGCAGTCTGGTGACTTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCCGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGtc                  <  8:25506/138‑1 (MQ=255)                                                                                                              ccGCAGTCTGGTGACTTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCCGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGtc                  <  8:30106/138‑1 (MQ=255)                                                                                                               cGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTtgct                                                                                                       <  2:101267/52‑1 (MQ=255)                                                                                                               cGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTtgct                                                                                                       >  1:101267/1‑52 (MQ=255)                                                                                                                        gTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATATTTTTCTTTCCGCTGCTTTTTTTCATTTTGGAGTCAGTTTTTTGTTTGATATCCACTGTTCCGTCTTGTTTGTTCa       >  6:154018/1‑139 (MQ=255)                                                                                                                            ccTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGTCTTCTTTGTTCACCt    >  6:82071/1‑138 (MQ=255)                                                                                                                             cTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGTCTTCTTTGTTCACCTTa  <  5:299928/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                           ATACAGCACCGGTTAACACACCTGTCACTGCCGCGTCCTGAATTCCGATTACTGTAACCCATTCTAGTTTTGTGATATGTAAATGCGCGAGAAAAGACGCGCTGATTTTCCGCAGTCTGGTGACCTGATGCAAAAGCATTTCGATTTTATTAATGCTTTGCTGCACATCGTTTTGTGTAATTTTTTTCTTTCCGCTGCTTTTTTTCATTTTGGATTCAGTTTTTTGTTTGATATCCACTGTTCCGTCTTCTTTGTTCACCTTA  >  BS168_reference_genome/3019026‑3019288

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