Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,935,823 Δ1 bp 100% intergenic (+35/+1) ywbE → / ← rlmI conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,935,8230T.100.0% 123.4 / NA 28intergenic (+35/+1)ywbE/rlmIconserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/13);  total (15/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935705‑3935953
                                                                                                                       |                                                                                                                                  
aaaTCAAACTTCAAACCGCATGGCATTAAAGTCCGTCTTGAAGTTGGCCGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                               <  5:57048/139‑1 (MQ=255)
           tCATCCGCATGGCATTAAAGTCCGTCTTGAAGAAGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTATTTCTGCAGAAGGACCACTTTCAgg                                                                                                     >  5:107725/1‑138 (MQ=255)
                                 cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:123781/1‑138 (MQ=255)
                                 cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCGTTTTCATTATTTTTTTTGCAGAAGGACCACTTTGAGGTAGTTGCCTTGCGGAAATGGCg                                                                               >  3:26829/1‑138 (MQ=255)
                                 cGTCTTGAAGATGGCCGCATCGGCCGCGACCACCAGCTCGATACTACATAAAAAAAAGCAGCCCCGGAACAGGGCAGCGTTTTCATTAATTTTAATGCAGAAGGCACCCTTTCAGGTAGTATCCTTCAGGAAATgacc                                                                               >  6:168717/1‑135 (MQ=255)
                                             ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  1:117259/1‑138 (MQ=255)
                                              gCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGGTTTGaaa                                                                 >  1:143803/1‑139 (MQ=255)
                                                      ccccGCGACCAGAATATCGTTCAAACAAAAAAAAACGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCgg                                                         <  2:143803/137‑1 (MQ=255)
                                                        ccGCGACCCGAATTTCGTTTCTAAATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCa                                                       <  2:117259/139‑1 (MQ=255)
                                                            gACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCg                                                    >  4:204653/1‑138 (MQ=255)
                                                                agcagATCGTTTCTACATAAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTaaa                                                <  3:131555/139‑1 (MQ=255)
                                                                  cagATCTTTTCTAAATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             <  7:251845/139‑1 (MQ=255)
                                                                  cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                               <  2:123781/137‑1 (MQ=255)
                                                                        gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  2:127342/138‑1 (MQ=255)
                                                                          ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGa                                     >  1:7931/1‑139 (MQ=255)
                                                                               cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                <  5:204772/139‑1 (MQ=255)
                                                                                   aaaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAg                             >  5:76608/1‑139 (MQ=255)
                                                                                               ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTACGGCAGCGTAAATTCCTCAATGATCGTGTAGCTTTAAttttt                  <  3:204653/138‑3 (MQ=255)
                                                                                               ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCAGCAGAAGGACCACTTTCAGGTAGTTTCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTAAATGATCTTGTAGCGTTGATTTGtt                 >  8:31019/1‑139 (MQ=255)
                                                                                                         gggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   <  7:289995/79‑1 (MQ=255)
                                                                                                         gggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  8:289995/1‑79 (MQ=255)
                                                                                                          ggCATCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCCTAAATTACTCAATGATCGTGTAGCGTTCATTTGTTTATTTAAAg        >  2:4804/1‑137 (MQ=255)
                                                                                                          ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAACTCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  2:243551/1‑139 (MQ=255)
                                                                                                          ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAg        >  6:262381/1‑137 (MQ=255)
                                                                                                          ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGGAGCGTTCATTTGTTTCTTTAAAggc      >  5:292762/1‑139 (MQ=255)
                                                                                                             aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  7:127931/139‑1 (MQ=255)
                                                                                                             aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  7:31019/139‑1 (MQ=255)
                                                                                                              gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  5:262381/139‑1 (MQ=255)
                                                                                                                       |                                                                                                                                  
AAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935705‑3935953

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: