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_genome	3,935,823	0	T	.	100.0%	217.0 / NA	48	intergenic (+35/+1)	ywbE/rlmI	conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (33/15);  total (33/15)

GTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC  >  BS168_reference_genome/3935687‑3935952                                                                                                                                         |                                                                                                                                  gTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatta                                                                                                                                 >  8:332367/1‑138 (MQ=255) gTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatta                                                                                                                                 >  7:110129/1‑138 (MQ=255)            tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                       >  8:71344/1‑137 (MQ=255)            tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                       >  8:275844/1‑137 (MQ=255)            tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                       >  3:343077/1‑137 (MQ=255)                   aaaTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCATGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGAc                                                                                                              >  2:208665/1‑139 (MQ=255)                                     cATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGcc                                                                                            >  6:167776/1‑139 (MQ=255)                                            ttAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGa                                                                                     <  5:167776/139‑1 (MQ=255)                                               aaGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAAt                                                                                  >  8:122490/1‑139 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACATGGCAGCTTTTTCA‑TTATTTTTTCTGAAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  7:231931/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  5:285080/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  8:44432/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  7:114843/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  5:76093/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTATGCAGAAGGACCACATTAAGGTAGTTGCCGTCCGGAAATGCCg                                                                              >  2:369265/1‑138 (MQ=255)                                                    cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTATCTACATAAAAAAAAGCTGCACCTGAACAGGGCAGAGTATTCA‑TTATTATTTCTGCAGAATGACCACTTTCAGTTAGTTTCCTACCGGTAATGcag                                                                              >  6:141647/1‑136 (MQ=255)                                                                ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACCGGGCAGCTTTTTGA‑TTATTTTTTCTGCAGAAGGCCCACTTTCAGGTAGTTTCCTTCCGGAAATGCCGAAATCGtttt                                                                    >  3:227496/1‑136 (MQ=255)                                                                     cATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGCACAGGGCAGCTTTTTGA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAGAtcct                                                            >  1:679/1‑139 (MQ=255)                                                                     cATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcct                                                            >  4:217605/1‑139 (MQ=255)                                                                     cATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcct                                                            >  6:196655/1‑139 (MQ=255)                                                                        cGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCg                                                         <  8:114843/139‑1 (MQ=255)                                                                         ggCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCgg                                                        <  7:71344/139‑1 (MQ=255)                                                                          gCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                       <  8:110129/139‑1 (MQ=255)                                                                              cgACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGa                                                                                                               <  2:120914/79‑1 (MQ=255)                                                                              cgACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGa                                                                                                               >  1:120914/1‑79 (MQ=255)                                                                                 ccAGCAGATCGTTTCTACATAAAAAAACGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTaa                                                <  1:369265/139‑1 (MQ=255)                                                                                   agcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                              <  6:76093/139‑1 (MQ=255)                                                                                            tttCTACATAAAAAAAAGCCGCCCCGGATCAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGGAGTTGCCTTCCGGCAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATg                                     <  4:338295/139‑1 (MQ=255)                                                                                                  cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                               >  3:299863/1‑139 (MQ=255)                                                                                                    tAAAAACACGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                             <  7:27695/139‑1 (MQ=255)                                                                                                    tAAAAAAAAGCTGCCCCTGAAAAGGGCAGCTTTTTGA‑TTATTTTTTCGGCAGAAGGCCCCCATTCAGGTAGTTGGCTTACGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                             >  2:296532/1‑139 (MQ=255)                                                                                                      aaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                              >  3:284300/1‑120 (MQ=255)                                                                                                      aaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                              <  4:284300/120‑1 (MQ=255)                                                                                                                 cccaTGAACAGGGCAGCTTTTTCA‑TTATTTTTGCTGCAGTTCTACCACTTTCAGGTAGTTGCCTTCCGGAAATGACGAAATCGTTTTGAATGCCTCCGGCAGCGTAAATTCGTCATTGATCGTGAAGCGTTCATTTGt                 <  5:141647/138‑1 (MQ=255)                                                                                                                           gggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTAAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTATTGAAATCCTCCGCCAGCGAAAATTCATCAATGATCGTCTAGCGATCAATTGTTTCTTAAAAgt      >  5:123448/1‑138 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCCGAAGGACCACTTTCAGGTAGTTGCCTTCCGGTAATGCCGAAATCGTTTTGGAATCATCCGGCAGCGTAAATTCCTCAATGATCGTGTAGAGGTACGTTGGTTCTTTAATggg     >  2:379039/1‑138 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                             <  8:23088/115‑1 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                             >  7:23088/1‑115 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGt                         >  3:144343/1‑119 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  3:331214/1‑139 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  5:372271/1‑139 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  2:224351/1‑139 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  7:31984/1‑139 (MQ=255)                                                                                                                            ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTAAATGATCGTGTAGCGt                         <  4:144343/119‑1 (MQ=255)                                                                                                                              cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg   >  5:204537/1‑139 (MQ=255)                                                                                                                              cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg   >  6:308908/1‑139 (MQ=255)                                                                                                                              cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg   >  4:162615/1‑139 (MQ=255)                                                                                                                              cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACAACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg   >  1:43382/1‑139 (MQ=255)                                                                                                                               aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  7:44432/139‑1 (MQ=255)                                                                                                                               aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  5:308908/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC  >  BS168_reference_genome/3935687‑3935952

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