Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	CP009273	2,656,840:1	(GCGCGCA)2→3	81.2%	coding (40/804 nt)	suhB →	inositol monophosphatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 2656826	3 (0.140)	15 (0.670) +GCGCGCA	16/228	NT	81.2%	coding (26/804 nt)	suhB	inositol monophosphatase
	?	CP009273	2656827 =	4 (0.170)					coding (27/804 nt)	suhB	inositol monophosphatase

CATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAT  >  CP009273/2656680‑2656972                                                                                                                                                    |                                                                                                                                                   tttcttcatctaaagattattcaggcatgttatcgtaaaacgaagacagatgcagatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcagc                                                                                                                                                   >  2:239432‑M2/148‑149 (MQ=255)                   attcacgcatcttatcataaaacgaagacagatgccgatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgctgatgacacaccgggtgcgcgcaGCGCGCAGCGCGCAAGGCgg                                                                                                                                 <  2:212039‑M2/20‑1 (MQ=255)                      gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCGGATCTCGCTGCTATACTCTGCGCCGTTTTCCTGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGATCATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                    >  1:102379/1‑143 (MQ=255)                      gATTATTCACGCATCTTATCAGAAAACGAAGACAGATGCCGATCTCGCTGCGATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGta                                                                                                                              <  2:120146/149‑1 (MQ=255)                                   ataaaacgaagacagatgccgatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGCCGGGAAATTTCAATGCCaa                                                                                                                 >  2:110710‑M2/114‑149 (MQ=255)                                   ataaaacgaagacagatgccgatctcgctgctatactctgcgccgttttcccggtctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGAACTTTACTTGACaa                                                                                                                 >  2:301344‑M2/114‑149 (MQ=255)                                    taaaacgaagacagatgccgatctcgctgctatactttgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                      <  1:309513‑M2/31‑1 (MQ=255)                                     aaaacgaagacagatgccgatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                                                                     >  1:46742‑M2/112‑143 (MQ=255)                                          gaagacagatgccgatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGAGCAGCGCCCCAGGCGGGTCATTTCATTGCCcaa                                                                                                                >  1:180550‑M2/107‑143 (MQ=255)                                                   tgccgatctcgctgctatactctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGaa                                                                                                       <  1:123134‑M2/46‑1 (MQ=255)                                                    gccgatctcgctgctatactctgcgccgtttacccgttctttaacatccagtgagcgagaccgatgcatccgatgctgaacatcgccgtgcgcggaGCGCGCAGCGCGCAAGGCGGGTAAGTtaattgccaaaaaacttgtaa                                                                                                      >  1:55445‑M2/97‑135 (MQ=255)                                                         tctcgctgctataccctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCt                                                                                           <  2:29685‑M2/58‑1 (MQ=255)                                                                     actctgcgccgttttcccgttctttaacatccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                     <  1:84221‑M2/64‑1 (MQ=255)                                                                     actctgcgccgttttcccgttctttaacagccagtgagagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                     <  1:125935‑M2/64‑1 (MQ=255)                                                                                             cGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGCAACCCCGGACGCTGTAGAAGCGAGCCAGCAAGGCAGTAAC‑GATTTCGTg                                                      >  2:285115/1‑149 (MQ=255)                                                                                                         agagggaccgatgcatccgaggctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATGGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGAt                                          <  2:60246‑M2/106‑1 (MQ=255)                                                                                                         agagagaccgatgcatccgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGAt                                          >  2:32746‑M2/44‑149 (MQ=255)                                                                                                                       atccgaggctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGc                            <  2:53174‑M2/120‑1 (MQ=255)                                                                                                                          cgatgctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                         >  2:202348‑M2/27‑149 (MQ=255)                                                                                                                              gctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGAtt                     >  2:197622‑M2/23‑149 (MQ=255)                                                                                                                               ctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCCACGTCGATACAGCTGCCGAAGCGGTGATTa                    >  2:30291‑M2/22‑149 (MQ=255)                                                                                                                               ctgaacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GAATTCGTGACCAACGTAGATAAAGCtgccgagtggg                          >  1:241873‑M2/22‑138 (MQ=255)                                                                                                                                  aacatcgccgtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGTCAGAAAGGCAGTAAC‑GATTTCGTGGCCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCg                 >  2:196584‑M2/19‑149 (MQ=255)                                                                                                                                       tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGCAGCGGTGATTATc                  >  1:50099/1‑143 (MQ=255)                                                                                                                                           gtgcgcgcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGaca              >  1:268160‑M2/10‑143 (MQ=255)                                                                                                                                               gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTaa    >  2:375613/1‑149 (MQ=255)                                                                                                                                                 gcaGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTc        >  1:257537‑M2/4‑143 (MQ=255)                                                                                                                                                 cgggcgcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCTGTGATTATCGACACGATTCGTAAAt  <  2:258664/146‑1 (MQ=255)                                                                                                                                                 cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACACCGGCCGACGTAGAAGCGAGCCAGAAAGGAAGCAACGGGCGGCCTG‑CCGACGTAAACAAATCGGCCGAAGCGGTGATTATCGACACGAtcc        >  1:374330/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAAC‑GATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAAT  >  CP009273/2656680‑2656972

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.