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.97.8% 182.5 / ‑2.1 46intergenic (+35/+1)ywbE/rlmIconserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (27/18);  minor base A (1/0);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935698‑3935953
                                                                                                                             |                                                                                                                                  
tttGACGAAATCACACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTGTTTCA‑TTATTTTTTATg                                                                                                                        >  7:114959/1‑137 (MQ=255)
tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                        >  7:14537/1‑137 (MQ=255)
tttGACGAAATCAAACTTTCATCAGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTATCTACATAAAAAAAAGCTGCCCATGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                        >  4:21364/1‑137 (MQ=255)
       aaaTCAAAATTTCATCCCCATGCAATTAATGTCGGTCTTGATGATGGCCCCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               <  7:284426/139‑1 (MQ=255)
                  tCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAgg                                                                                                     >  3:231308/1‑138 (MQ=255)
                    aTCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAg                                                                                                  >  1:108354/1‑139 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:199707/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  6:88467/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGAAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  5:263424/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTT‑AATTATTTTTTCTGCAGAAGGACAACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  7:39455/1‑138 (MQ=255)
                                          tCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                            >  5:34471/1‑139 (MQ=255)
                                                                    aCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCACCGTTTTCA‑TTAAATTTTCTGCAGAAGGACCACTTTCATGTAGTTGCCTTCCGGAAATGCAGAAATCGTATTGAAATGCTCCGGAAGCGt                                                   >  5:49977/1‑138 (MQ=255)
                                                                         cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             <  3:231228/139‑1 (MQ=255)
                                                                         cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             <  6:151256/139‑1 (MQ=255)
                                                                         cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             <  8:130354/139‑1 (MQ=255)
                                                                         cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                             >  4:34560/1‑139 (MQ=255)
                                                                               gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  6:49977/138‑1 (MQ=255)
                                                                               gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  8:60278/138‑1 (MQ=255)
                                                                               gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  2:44181/138‑1 (MQ=255)
                                                                               gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  4:231308/138‑1 (MQ=255)
                                                                               gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGAAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  6:31595/138‑1 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               >  1:189390/1‑60 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               <  2:189390/60‑1 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                >  6:254078/1‑139 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGaaa                                                                 >  8:217902/1‑81 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATTAACGTTTTGGGTTCATTTTTTTTTTTaaatg       >  2:183513/1‑137 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg       >  7:223694/1‑139 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTATCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGaaa                                                                 <  7:217902/81‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGCCCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGAGTAAATTCCTCAATGAACGAGTAGCGTTCATTTGTTTCTTTAAAggc      >  2:2098/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCATGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  <  4:70074/127‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                            >  1:76954/1‑69 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                            <  2:76954/69‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  8:202111/105‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        >  7:202111/1‑105 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATg                                      <  5:290875/107‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATg                                      >  6:290875/1‑107 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  <  6:56413/127‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  >  5:56413/1‑127 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTTAAggc      >  6:175290/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  4:89053/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATGTGTTTATTTaactgc      >  5:157352/1‑135 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTGCATTTGTTTCTTTAAAggc      >  5:161417/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGGTGATTTGTTTCTTTAAAggc      >  2:204282/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCACAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  >  3:70074/1‑127 (MQ=255)
                                                                                                                   aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  8:223694/139‑1 (MQ=255)
                                                                                                                    gCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  3:89053/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                  
TTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935698‑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: ‑