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

ATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCA  >  CP009273/3320760‑3321193                                                                                                                                                                                                         |                                                                                                                                                                                                                                          aTCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACAAATTTATCGCTAAAGTGTTCCTGaa                                                                                                                                                                                                                    >  2:515530/1‑224 (MQ=255)     aTAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTTAACATAttta                                                                                                                                                                                                                                       >  1:267250‑M1/1‑197 (MQ=255)     aTAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCa                                                                                                                                                                                                                <  1:515540/224‑1 (MQ=255)                    cAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCATGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcttgaagccagcggatggagctggcagaacgcttcttacctgtcctttaactttcncatgggtataac                                                                                                                                                >  2:490894‑M1/1‑182 (MQ=255)                                   ccTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAATCCCTTTTTCTGAGCCTGTTCAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGTGGCTGGAGCTGGCAGAACGcttct                                                                                                                                                                                      >  1:616900/1‑220 (MQ=255)                                        gccgccAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAg                                                                                                                                                                                                                   >  1:12799/1‑186 (MQ=255)                                            ccAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGTTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAg                                                                                                                                                                                                               <  2:12799/186‑1 (MQ=255)                                                          cacaTATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGACAGAACGCTGATTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTAATTGGTGat                                                                                                         >  1:335791/1‑274 (MQ=255)                                                           acaTATTTTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAANTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttcttacctgtcatttaactttcccatggggataactcatcgtaaccaattgcgtaaatttatacgcgcctatttgg                                                                                                       <  1:481545‑M1/275‑133 (MQ=255)                                                             atatTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTNTAAACCATGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttcttacctgtcatttaactttcccatggggataactcatcgtaaccaattgcgtaaatttttacgcgcctatttggtg                                                                                                     <  1:490902‑M1/275‑135 (MQ=255)                                                                   tgACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAttta                                                                                                                                                                                                                                       >  1:32396‑M1/1‑135 (MQ=255)                                                                   tgACCAACCACCCGGAGCAGCACCAAGGCCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcct                                                                                                                                                                                                                 >  1:516386‑M1/1‑135 (MQ=255)                                                                       caaccaCCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaag                                                                                                                                                                                                             <  2:516376‑M1/161‑31 (MQ=255)                                                                       caaccaCCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctt                                                                                                                                                                                                                                   <  2:32396‑M1/139‑9 (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‑M1/1‑89 (MQ=255)                                                                                                                     tAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttcttacctgtcatttaactttcccat                                                                                                                                                         <  1:69212‑M1/167‑83 (MQ=255)                                                                                                                            tttGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttct                                                                                                                                                                                >  2:258223‑M1/1‑78 (MQ=255)                                                                                                                                tCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttcttacc                                                                                                                                                                            <  1:258227‑M1/137‑64 (MQ=255)                                                                                                                                          gTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAAAATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCt                                                                                                                                                                                                          >  1:118537/1‑97 (MQ=255)                                                                                                                                              ttCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAg                                                                                                                                                                                                      <  2:118536/97‑1 (MQ=255)                                                                                                                                                catcaAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctgg                                                                                                                                                                                            >  2:335845‑M1/1‑58 (MQ=255)                                                                                                                                                  tcaAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgct                                                                                                                                                                                                                            >  2:399350‑M1/1‑56 (MQ=255)                                                                                                                                                    aaGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcaga                                                                                                                                                                                        <  1:335850‑M1/105‑52 (MQ=255)                                                                                                                                                      gTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaag                                                                                                                                                                                                                        <  1:399355‑M1/71‑20 (MQ=255)                                                                                                                                                           aaCCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagcnggcagaacgcttcttacctgtcatttaactttcccatggggataactcatcgtaaccaattgcgtaaatttttacgcgcctatttggtgatatatgggagatggcggtagaatgacccgttttcaatcccaacgtaagaaaaaaatacgatgaagctgagtactaatctaaaacagcacctga       <  1:491944‑M1/275‑229 (MQ=255)                                                                                                                                                                ggCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagctggcagaacgcttcttacctgtcatttaactttcccatggggataactcatcgtaaccaattgcgtaaatttttacgcgcctatttggtgatatatgggagatggcggtagaatgacccgttttcaatcccaacgtaagacaaaaatacgaagaatgtgagtacgaaacaaaaanaacaccgtaacggc  >  2:491936‑M1/1‑42 (MQ=255)                                                                                                                                                                 gCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGc                                                                                                                                                                                                 >  1:411822/1‑83 (MQ=255)                                                                                                                                                                    cGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATAttt                                                                                                                                                                                                                                        <  1:302490‑M1/41‑4 (MQ=255)                                                                                                                                                                     gggAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTCTCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGaa                                                                                                                                                                                             <  2:411817/83‑1 (MQ=255)                                                                                                                                                                          cGTAACCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctg                                                                                                                                                                                                   >  2:31986‑M1/1‑32 (MQ=255)                                                                                                                                                                            tAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTgg                                                                                                             <  2:325613/156‑1 (MQ=255)                                                                                                                                                                              aCCCCTTTTTCTGTGCCTGTTGAACATAtttatctttatcgctaaagtgttcctgaagccagcggctggagc                                                                                                                                                                                               <  1:31986‑M1/72‑45 (MQ=255)                                                                                                                                                                                                         |                                                                                                                                                                                                                                          ATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACCCGGAGCAGCACCAAGGTCGACAACCGTCATTCCCGGTTTAAAGAGTTTGTCACTTTGCTGTATTTCATCAAGTTTAAACCAGGCACGGGAACGTAACCCCTTTTTCTGTGCCTGTTGAACATATTTATCGCTAAAGTGTTCCTGAAGCCAGCGGCTGGAGCTGGCAGAACGCTTCTTACCTGTCATTTAACTTTCCCATGGGGATAACTCATCGTAACCAATTGCGTAAATTTTTACGCGCCTATTTGGTGATATATGGGAGATGGCGGTAGAATGACCCGTTTTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCA  >  CP009273/3320760‑3321193

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.