Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	CP000730	62,818	+CAATTT	intergenic (+83/+228)	USA300HOU_0056 → / ← speG	hypothetical protein/diamine N‑acetyltransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000730	= 62817	0 (0.000)	26 (0.690) +TCAATT	22/266	0.4	100%	intergenic (+82/+229)	USA300HOU_0056/speG	hypothetical protein/diamine N‑acetyltransferase
	?	CP000730	62818 =	0 (0.000)					intergenic (+83/+228)	USA300HOU_0056/speG	hypothetical protein/diamine N‑acetyltransferase

AAAAATAGAATTATTAAAACAACGCAAACAAGGCTTACTACAAAAAATGTTTGTGTAATTCTGATAACGTCCTATTATGTTAATTGAAAATAGACAAAAATAAATTGAGAGCATCCCTTACCGCAAAGTGAAGGATGCTCTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTCGTAGATTTTATCAAGAAT  >  CP000730/62678‑62945                                                                                                                                             |                                                                                                                                agaattattaaaacaacgcaaacaaggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattt                                                                                                                                 >  1:34839‑M2/141‑141 (MQ=255)       attaaaacaacgcaaacaaggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTAttttt                                                                                                                           <  1:4764‑M2/7‑1 (MQ=255)            aacaacgcaaacaaggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagaaaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGa                                                                                                                      >  1:227167‑M2/130‑141 (MQ=255)             acaacgcaaacaaggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGaa                                                                                                                     >  1:199670‑M2/129‑141 (MQ=255)                acgcaaacaaggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAAtt                                                                                                                  <  2:34839‑M2/16‑1 (MQ=255)                         aggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGc                                                                                                         <  2:114015‑M2/25‑1 (MQ=255)                          ggcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCt                                                                                                        <  2:127465‑M2/26‑1 (MQ=255)                           gcttactacaaaaaatgtttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCtt                                                                                                       >  2:187910‑M2/115‑141 (MQ=255)                                            tttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTc                                                                                      >  1:399366‑M2/98‑141 (MQ=255)                                             ttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCa                                                                                     >  1:143246‑M2/97‑141 (MQ=255)                                             ttgtgtaattctgataacgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCa                                                                                     >  1:16814‑M2/97‑141 (MQ=255)                                                             acgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTttgaa                                                                     <  1:379709‑M2/61‑3 (MQ=255)                                                              cgtcctattatgttaattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCaa                                                                    >  1:56996‑M2/80‑141 (MQ=255)                                                                           taattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGt                                                       <  1:72420‑M2/75‑1 (MQ=255)                                                                            aattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTa                                                      >  2:55836‑M2/66‑141 (MQ=255)                                                                            aattgaaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTa                                                      >  2:244049‑M2/66‑141 (MQ=255)                                                                                 aaaatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCgtt                                                 >  1:427498‑M2/61‑141 (MQ=255)                                                                                   aatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCgttgt                                               >  1:47918‑M2/59‑141 (MQ=255)                                                                                   aatagacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAACATAAGCTTCATTCTGATAGATAGTCATTCATAATTGTTTGCAAATCGTTATTAAGTACCgttgt                                               <  1:380183‑M2/83‑1 (MQ=255)                                                                                      agacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAACTTAATATAAGCTTCATTCTGATCGATAGTCATTCCTCCTTGTTTGTCAATCGTTATTAAGTAccgttgttgc                                            >  1:19495‑M2/56‑139 (MQ=255)                                                                                      agacaaaaataaattgagagcatcccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTAc                                            >  1:253424‑M2/56‑141 (MQ=255)                                                                                                              ccttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTc                    >  2:29604‑M2/32‑141 (MQ=255)                                                                                                               cttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTTATCATTCATAAAATCCTCg                   >  2:14232‑M2/31‑141 (MQ=255)                                                                                                               cttaccgcaaagtgaaggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTCg                   <  1:187910‑M2/111‑1 (MQ=255)                                                                                                                              aggatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTCGTAGATTTTATCAAGa    >  2:378561‑M2/16‑141 (MQ=255)                                                                                                                                gatgctctcaattTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTCGTAGATTTTATCAAGAAt  <  2:199670‑M2/128‑1 (MQ=255)                                                                                                                                             |                                                                                                                                AAAAATAGAATTATTAAAACAACGCAAACAAGGCTTACTACAAAAAATGTTTGTGTAATTCTGATAACGTCCTATTATGTTAATTGAAAATAGACAAAAATAAATTGAGAGCATCCCTTACCGCAAAGTGAAGGATGCTCTATTTTTGTTGAAATTAATATAAGCTTCATTCTGATAGATAGTCATTAATAATTGTTTGCAAATCGTTATTAAGTACCGTTGTTACAGCGTCATCATTCATAAAATCCTCGTAGATTTTATCAAGAAT  >  CP000730/62678‑62945

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