Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,910,754 C→T 19.9% S196S (TCC→TCT yndF → putative spore germination lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,910,7540CT19.9% 43.6 / 3.1 20S196S (TCC→TCTyndFputative spore germination lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (9/7);  new base T (2/2);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGGCTGAGCTGTCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAGCGGTCAAGTATTCCGGAGCAGCTGTAATCAACGGAAAATCG  >  BS168_reference_genome/1910627‑1910889
                                                                                                                               |                                                                                                                                       
aaaTCAGGCTGAGCTGTCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGATTTCTTcc                                                                                                                              <  1:287016/139‑1 (MQ=255)
aaaTCAGGCTGAGCTGTCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAAttcttc                                                                                                                               <  3:33543/138‑1 (MQ=255)
     aGGCTGAGCTGTCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTa                                                                                                                         <  5:293849/139‑1 (MQ=255)
                tCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTcgc                                                                                                               <  2:21890/138‑1 (MQ=255)
                          tCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTTCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCTAGAATTCTTcc                                                                                                                              <  1:286845/113‑1 (MQ=255)
                          tCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTTCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCTAGAATTCTTcc                                                                                                                              >  2:290441/1‑113 (MQ=255)
                          tCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTTCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCTAGAATTCTTcc                                                                                                                              >  2:286845/1‑113 (MQ=255)
                          tCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTTCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCTAGAATTCTTcc                                                                                                                              <  1:290441/113‑1 (MQ=255)
                          tCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTg                                                                                                    >  8:83674/1‑139 (MQ=255)
                                        aGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCt                                                                                      >  6:113911/1‑139 (MQ=255)
                                              agaCGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAAttctt                                                                                                                                <  3:271017/91‑1 (MQ=255)
                                              agaCGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAAttctt                                                                                                                                >  4:271017/1‑91 (MQ=255)
                                                aCGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGcagcag                                                                              >  6:240593/1‑139 (MQ=255)
                                                   ctctCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCtt                                                                           >  2:120717/1‑139 (MQ=255)
                                                                   gaACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAAtgtg                                                           >  4:18721/1‑139 (MQ=255)
                                                                      cGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGtt                                                        >  4:202813/1‑139 (MQ=255)
                                                                                   gCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAg                                           <  7:83674/139‑1 (MQ=255)
                                                                                       ttCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAGCGGt                                       >  1:26517/1‑139 (MQ=255)
                                                                                                     ttAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAGCGGTCAAGTATTCCGGag                         <  3:202813/139‑1 (MQ=255)
                                                                                                                             tCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAGCGGTCAAGTATTCCGGAGCAGCTGTAATCAACGGAAAATCg  >  5:255872/1‑138 (MQ=255)
                                                                                                                               |                                                                                                                                       
AAATCAGGCTGAGCTGTCTTGTTTTATCTGCTAAAGGGGAAGCAAGAGACGCTCTCCAGCTGAAAGAGAACGGAGAAATTCCGGCATTTCGGCTGATTGGATTAGGTGAAAATGAACATAAAGTTTCCAGAATTCTTCCCCCTATGACACTCGCGAAACTAATTGGAAAACTGCATTCTGGCAGCAGCTTTTTACTGCAAAATGTGGTTGCGGCGAATGGAGCGGTCAAGTATTCCGGAGCAGCTGTAATCAACGGAAAATCG  >  BS168_reference_genome/1910627‑1910889

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: