Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 3,935,825 Δ1 bp coding (1190/1191 nt) rlmI ← 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.98.1% 230.9 / ‑2.5 53intergenic (+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 . (32/20);  minor base C (1/0);  total (33/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 = 1.00e+00

TTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCA  >  BS168_reference_genome/3935698‑3935953
                                                                                                                             |                                                                                                                                  
tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTATg                                                                                                                        >  3:301739/1‑137 (MQ=255)
tatGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCATGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                        >  2:248612/3‑137 (MQ=255)
    aCGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCACATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCGGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAg                                                                                                                  >  5:62124/1‑139 (MQ=255)
       aaaTAAAACTTTCATCCGCATGGCATTAAAGCCCGGCTTTAATATGGCCGCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               <  8:2476/139‑1 (MQ=255)
         aTCAAATTTTAATCCGCATGTCATTAACTCCCGTTTTTATGTTGGCCGCATCGGCCGCGACCAGCATATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCa                                                                                                             <  4:140447/139‑1 (MQ=255)
         aTCAAACTTTCTACCGCCTGCATTTAATGTCGGTGTTGAATGAGGCCGCATCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCa                                                                                                             <  8:148463/139‑1 (MQ=255)
                   cATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAgg                                                                                                     >  5:171301/1‑137 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGCCCCGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  4:327992/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAATTGCCg                                                                               >  8:321818/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:192239/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTATGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  5:213595/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTT‑GCTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                               >  1:160441/1‑138 (MQ=255)
                                        cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGAACGTTTCTACATAAAAAAAAGCTGCCACTGAACAGGGCAGCTATTTCA‑GTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTTCCTTACGTAAATGCCg                                                                               >  7:229642/1‑138 (MQ=255)
                                          tCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                            >  2:47936/1‑139 (MQ=255)
                                          tCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTGTTTGA‑TTATTTTTTCTGCAGAAGGCCAACTTTCAGGTAGTTGCCTTCCGGGAATGCCGaa                                                                             >  5:213912/1‑138 (MQ=255)
                                                    ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  3:159990/1‑138 (MQ=255)
                                                    ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  1:242476/1‑138 (MQ=255)
                                                      ccGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTTCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAt                                                                >  2:80883/1‑139 (MQ=255)
                                                              gCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGCCCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                        >  4:212546/1‑139 (MQ=255)
                                                              cccGCGACCCGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTATTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                        <  4:19388/138‑1 (MQ=255)
                                                                   gACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGt                                                   <  1:47936/139‑1 (MQ=255)
                                                                       agcatATAGTTTCAAAAAAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                               <  6:213912/139‑1 (MQ=255)
                                                                                 ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGa                                     >  5:109109/1‑139 (MQ=255)
                                                                                  tCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCTGAAGGACCACTTTCAGGTAGTTGCCTGCCGGAAATGCCGAAATCGTTTTGAAATCCTACGGCAGCGTAAATGCCTCAATGAt                                    >  7:220035/1‑139 (MQ=255)
                                                                                  tCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGa                                     <  6:237689/138‑1 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                <  7:5729/139‑1 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                <  7:10662/139‑1 (MQ=255)
                                                                                      cATAAAAAAAAGCTGCCACTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                <  6:219094/139‑1 (MQ=255)
                                                                                           aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               <  7:78569/55‑1 (MQ=255)
                                                                                           aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                               >  8:78569/1‑55 (MQ=255)
                                                                                           aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                              <  7:137598/136‑1 (MQ=255)
                                                                                           aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                              >  8:137598/1‑136 (MQ=255)
                                                                                                aGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCCCTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGGGTTCAt                      >  2:189269/1‑139 (MQ=255)
                                                                                                     ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTTCAGAAGGACCACATTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATACTCCGTCAGCGTAAAATCCTCAATGATCGTGTAGCGTTCATTTGtt                 >  4:293535/1‑139 (MQ=255)
                                                                                                     ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTGCGGCAGCGTAAATTCCTCAATGATc                                   <  8:273173/121‑1 (MQ=255)
                                                                                                     ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTGCGGCAGCGTAAATTCCTCAATGATc                                   >  7:273173/1‑121 (MQ=255)
                                                                                                     ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGtt                 >  8:31599/1‑139 (MQ=255)
                                                                                                     ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                  <  2:242476/138‑1 (MQ=255)
                                                                                                      cccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGttt                >  7:10111/1‑139 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg       >  7:101592/1‑139 (MQ=255)
                                                                                                               gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTAATTTGTTTCTTAAACgg       >  7:199005/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATTCCGAAATCGTTTTGAAATCCTCAGGAAGCGTAAATTCCTCAATGAACGTGTAGCGTTAAATTGTATCATTAAAtgc      >  7:125886/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        >  6:128906/1‑105 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                        <  5:128906/105‑1 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  4:128941/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTAGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCATCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc      >  3:274936/1‑139 (MQ=255)
                                                                                                                ggCAGCTTTTTCA‑TTATTTTTTATGCAGAAGGCCCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTATCCTTCATTTGTTTCTTTAAAtgc      >  8:137741/1‑139 (MQ=255)
                                                                                                                  cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   <  2:131945/76‑1 (MQ=255)
                                                                                                                  cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                   >  1:131945/1‑76 (MQ=255)
                                                                                                                  cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg    <  7:137741/139‑1 (MQ=255)
                                                                                                                   aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  5:124896/139‑1 (MQ=255)
                                                                                                                   aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc   <  7:31599/139‑1 (MQ=255)
                                                                                                                    gCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCa  <  6:109109/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: ‑