Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	CP009273	3,320,960	+TTTATC	55.2%	coding (64/630 nt)	rlmE ←	23S rRNA U2552 2'‑O‑ribose methyltransferase, SAM‑dependent

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 3320960	15 (0.580)	18 (0.700) +TTTATC	17/202	NT	55.2%	coding (64/630 nt)	rlmE	23S rRNA U2552 2'‑O‑ribose methyltransferase, SAM‑dependent
	?	CP009273	3320961 =	15 (0.550)					coding (63/630 nt)	rlmE	23S rRNA U2552 2'‑O‑ribose methyltransferase, SAM‑dependent

ATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGT  >  CP009273/3320760‑3321187                                                                                                                                                                                                          |                                                                                                                                                                                                                                   aTCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACAAATTTATCGCTAAAGTGTTCCTGaa                                                                                                                                                                                                              >  2:515530/1‑224 (MQ=255)     aTAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCa                                                                                                                                                                                                          <  1:515540/224‑1 (MQ=255)              caagcgatgatgcggcctttgccgccaatttgggtgaccacatattgtgaccaaccacccggagcagcaccaaggtcgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaagtttaaaccatgcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCTTGAAGCCAGCGGATGGAGCTGGCAGAACGCTTCTTACCTGTCCTTTAACTTTCNCATGGGTATAAc                                                                                                                                                >  2:490894‑M2/189‑273 (MQ=255)                                   ccTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAATCCCTTTTTCTGAGCCTGTTCAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGTGGCTGGAGCTGGCAGAACGcttct                                                                                                                                                                                >  1:616900/1‑220 (MQ=255)                                        gccgccAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAg                                                                                                                                                                                                             >  1:12799/1‑186 (MQ=255)                                            ccAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGTTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAg                                                                                                                                                                                                         <  2:12799/186‑1 (MQ=255)                                                     acatattttgaccaaccacccggagcagcaccaaggtcgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaantttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTATACGCGCCTATTTgg                                                                                                       <  1:481545‑M2/126‑1 (MQ=255)                                                       atattgtgaccaaccacccggagcagcaccaaggtcgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaagtntaaaccatgcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTg                                                                                                     <  1:490902‑M2/128‑1 (MQ=255)                                                          cacaTATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGACAGAACGCTGATTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTAATTGGTGat                                                                                                   >  1:335791/1‑274 (MQ=255)                                                             tgaccaaccacccggagcagcaccaaggccgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCt                                                                                                                                                                                                                 >  1:516386‑M2/142‑161 (MQ=255)                                                                 caaccacccggagcagcaccaaggtcgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatctt                                                                                                                                                                                                                                   <  2:32396‑M2/2‑1 (MQ=255)                                                                 caaccacccggagcagcaccaaggtcgacaaccgtcattcccggtttaaagagtttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAg                                                                                                                                                                                                             <  2:516376‑M2/24‑1 (MQ=255)                                                                                      gcaCCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCg                                                                                                                                                                                                       >  1:158526/1‑146 (MQ=255)                                                                                       caCCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGaa                                                                                                                                                                                                              >  1:575271/1‑138 (MQ=255)                                                                                          cAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTg                                                                                                                                                                                                   <  2:158523/146‑1 (MQ=255)                                                                                           aaGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCa                                                                                                                                                                                                          <  2:575260/138‑1 (MQ=255)                                                                                                           ggtttaaagagtttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTc                                                                                                                                                             >  2:69212‑M2/96‑167 (MQ=255)                                                                                                               taaagagtttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCAt                                                                                                                                                         <  1:69212‑M2/76‑1 (MQ=255)                                                                                                                      tttgtcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGcttct                                                                                                                                                                                >  2:258223‑M2/85‑137 (MQ=255)                                                                                                                          tcactttgctgtatttcatcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTAcc                                                                                                                                                                            <  1:258227‑M2/57‑1 (MQ=255)                                                                                                                                          gTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAAAATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCt                                                                                                                                                                                                    >  1:118537/1‑97 (MQ=255)                                                                                                                                          catcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTgg                                                                                                                                                                                            >  2:335845‑M2/65‑105 (MQ=255)                                                                                                                                            tcaagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCt                                                                                                                                                                                                                            >  2:399350‑M2/63‑71 (MQ=255)                                                                                                                                              ttCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAg                                                                                                                                                                                                <  2:118536/97‑1 (MQ=255)                                                                                                                                              aagtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGa                                                                                                                                                                                        <  1:335850‑M2/45‑1 (MQ=255)                                                                                                                                                gtttaaaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAg                                                                                                                                                                                                                        <  1:399355‑M2/13‑1 (MQ=255)                                                                                                                                                     aaccaggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCNGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGAAAAAAATACGATGAAGCTGAGTACTAATCTAAAACAGCACCTGa       <  1:491944‑M2/222‑1 (MQ=255)                                                                                                                                                          ggcacgggaacgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGACAAAAATACGAAGAATGTGAGTACGAAACAAaaanaacaccgtaacggc  >  2:491936‑M2/49‑267 (MQ=255)                                                                                                                                                                 gCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGc                                                                                                                                                                                           >  1:411822/1‑83 (MQ=255)                                                                                                                                                                    cgtaacccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTg                                                                                                                                                                                                   >  2:31986‑M2/39‑72 (MQ=255)                                                                                                                                                                     gggAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTCTCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGaa                                                                                                                                                                                       <  2:411817/83‑1 (MQ=255)                                                                                                                                                                        acccctttttctgtgcctgttgaacatatttatcTTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGc                                                                                                                                                                                               <  1:31986‑M2/38‑1 (MQ=255)                                                                                                                                                                            tAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTgg                                                                                                       <  2:325613/156‑1 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                                                   ATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGT  >  CP009273/3320760‑3321187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.