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%	211.6 / NA	51	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 . (29/22);  total (29/22)

CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC  >  BS168_reference_genome/3935685‑3935952                                                                                                                                           |                                                                                                                                  cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGAGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCat                                                                                                                                   >  1:76089/1‑139 (MQ=255)  aGTGACGCATATTTCCCCCAAATCAAACTTTCACCCCCAGGCCATAAAAGCCCGCCTTGACGAGGGCCCCCCCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatt                                                                                                                                  <  3:63164/139‑1 (MQ=255)      aaGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatta                                                                                                                                 >  7:116891/1‑135 (MQ=255)              tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTg                                                                                                                       >  1:160487/1‑137 (MQ=255)                  accaaaTAAAACTTACATCCTCTTGCCATTAAAGTGCGAATTCATTTTGGCCGATTCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCa                                                                                                                     <  8:69232/132‑1 (MQ=255)                     aaaTCAAACTTGAATCCGCATGGCATTAACGGCCGTCTTGAAGATGGCCGCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAg                                                                                                                 <  8:169479/136‑1 (MQ=255)                       accAAACTTACACCCGCATGCCATAAACGCCCGTCTTCAAGATGGCCGCACCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGa                                                                                                               <  8:192819/134‑1 (MQ=255)                       aTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGATGGACCa                                                                                                            >  2:342609/1‑139 (MQ=255)                       aTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCa                                                                                                            >  3:350261/1‑139 (MQ=255)                       aTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCa                                                                                                            >  2:114549/1‑139 (MQ=255)                                     cGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAg                                                                                                 >  8:119419/1‑136 (MQ=255)                                                      cGTCTTGAAGATGGCCGCATCGGCCGCGACCCGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGACCAGGGCAGCCTTTTCA‑TTATTTTTTATGCAGAAGGACCACTTACATGTAGTTGCCTTCCGGAAAttcag                                                                              >  5:184956/1‑134 (MQ=255)                                                      cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGCACAGGGCAGCTTTTTCA‑TTATTTTTTATGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  4:354099/1‑138 (MQ=255)                                                      cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  3:251463/1‑138 (MQ=255)                                                      cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTATTCAGAATGACCACTTTAAGGGAGTTGCCTTGCGGCAATGACg                                                                              >  1:346791/1‑138 (MQ=255)                                                        tCTTGAAGATGGCCGCATCGGCAGCGACCAGCAGATCGTTTCTACAAAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTAGCCTTCCGGAAATGCCGaaa                                                                           <  2:218628/139‑1 (MQ=255)                                                         cTTGAAGATGGCCGCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                           <  5:166581/138‑1 (MQ=255)                                                         cTTGAAGATGGACGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAATCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                              >  7:197689/1‑135 (MQ=255)                                                                  ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                  >  3:254255/1‑138 (MQ=255)                                                                       cATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCGGACCAGGGCAGCTTTTTCA‑TTATTTTTTCTGCACAAGGGCCCCGTTCAGGTAGTTGCCTTCCGGAAATGCCGAAAACGTTTTGAAAtccg                                                            >  1:347182/1‑138 (MQ=255)                                                                            gCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcctcc                                                          <  8:116891/136‑1 (MQ=255)                                                                            gCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                       <  4:359701/139‑1 (MQ=255)                                                                             ccGCGACCAGCAGATCGTTTCTCCATAAAAAAACGCAGCCCCTGAACAGCGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcctcct                                                         <  8:197689/136‑2 (MQ=255)                                                                                    cagcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTaaa                                               <  2:76089/139‑1 (MQ=255)                                                                                     agcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAt                                              <  5:190695/139‑1 (MQ=255)                                                                                       cagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                            >  5:230164/1‑139 (MQ=255)                                                                                           tCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTATCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTc                                         >  3:353573/1‑138 (MQ=255)                                                                                               ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCt                                          <  8:33800/133‑1 (MQ=255)                                                                                               ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCt                                          >  7:33800/1‑133 (MQ=255)                                                                                                 cTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                       >  6:361102/1‑118 (MQ=255)                                                                                                 cTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                       <  5:361102/118‑1 (MQ=255)                                                                                                  tACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGAAGAAGGCCCACTAAAAGGTAGTTGCCTTCCGTAAATGCCGAAATCGTTTTGAAATCATCCGGCCGCGTAAATACCTCAATGATCg                                 >  6:96517/1‑139 (MQ=255)                                                                                                    caaaaaaaacaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                               <  1:183262/136‑1 (MQ=255)                                                                                                    ataaaaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                               <  5:101656/136‑1 (MQ=255)                                                                                                          aaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGt                         <  5:325171/139‑1 (MQ=255)                                                                                                                   ccccTGAACAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                 <  4:254255/138‑1 (MQ=255)                                                                                                                         aaCAGGGCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAACTCCTCCGGCAGCGTTAATTCCTCTATGCTCGGGTCTCGTTGATTTGTTTCttt           >  2:101762/1‑138 (MQ=255)                                                                                                                              ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  2:237308/1‑139 (MQ=255)                                                                                                                              ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  5:93532/1‑139 (MQ=255)                                                                                                                              ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  4:176628/1‑139 (MQ=255)                                                                                                                              ggCAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAAACGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  6:247531/1‑139 (MQ=255)                                                                                                                                cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                  >  2:224385/1‑76 (MQ=255)                                                                                                                                cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                  <  1:224385/76‑1 (MQ=255)                                                                                                                                cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTTGCGTTCATTTTTTTGTTTTAAggcgg   >  6:19364/1‑139 (MQ=255)                                                                                                                                cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg      >  8:216243/1‑136 (MQ=255)                                                                                                                                cAGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg      >  8:227483/1‑136 (MQ=255)                                                                                                                                cAGCTTATTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAACCTCCGGCAGAGTAAATTACTCAATTATCGTGTAGCGTTCATTTATTAGTTAAAAggcgg   >  4:332464/1‑139 (MQ=255)                                                                                                                                 aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  2:160487/139‑1 (MQ=255)                                                                                                                                 aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  1:326788/139‑1 (MQ=255)                                                                                                                                 aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  1:301959/139‑1 (MQ=255)                                                                                                                                 aGCTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  4:251463/139‑1 (MQ=255)                                                                                                                                 aGCTTTTTCA‑TTATTTTTTATGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc     >  8:215934/1‑136 (MQ=255)                                                                                                                                   cTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  >  4:66260/1‑137 (MQ=255)                                                                                                                                   cTTTTTCA‑TTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc  <  3:66260/137‑1 (MQ=255)                                                                                                                                           |                                                                                                                                  CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGC  >  BS168_reference_genome/3935685‑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: ‑