Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	1	1,568,439	+AATAA	coding (31/555 nt) coding (410/420 nt)	unknown ← unknown ←	Phosphoglycerate transport system transcriptional regulatory protein PgtA Phosphoglycerate transport system transcriptional regulatory protein PgtA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	1	= 1568439	0 (0.000)	23 (0.270) +AATAA	19/420	2.8	100%	coding (31/555 nt) coding (410/420 nt)	unknown unknown	Phosphoglycerate transport system transcriptional regulatory protein PgtA Phosphoglycerate transport system transcriptional regulatory protein PgtA
	?	1	1568440 =	0 (0.000)					coding (30/555 nt) coding (409/420 nt)	unknown unknown	Phosphoglycerate transport system transcriptional regulatory protein PgtA Phosphoglycerate transport system transcriptional regulatory protein PgtA

ATCAGTGATATGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTCGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCTTTTTCCG  >  1/1568197‑1568689                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           aTCAGTGATATGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatat                                                                                                                                                                                                                                                    <  1:789096‑M1/251‑9 (MQ=255)           tGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataa                                                                                                                                                                                                                                                       >  1:763249‑M1/1‑233 (MQ=255)                                                                       cTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcag                                                                                                                                                                              >  1:814587‑M1/1‑173 (MQ=255)                                                                         tATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcga                                                                                                                                                                                               <  1:88137‑M1/232‑62 (MQ=255)                                                                                                    aaGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtg                                                                                                                                                  >  1:253250‑M1/1‑144 (MQ=255)                                                                                                      gTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataa                                                                                                                                                                                                                                                       >  1:731207‑M1/1‑142 (MQ=255)                                                                                                         cTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatc                                                                                                                                            <  1:705590‑M1/251‑113 (MQ=255)                                                                                                                                  gAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaac                                                                                                                           <  1:149387‑M1/243‑130 (MQ=255)                                                                                                                                           ctcGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcgg                                                                                                                                                                                                            >  1:277990‑M1/1‑105 (MQ=255)                                                                                                                                                     tgTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgc                                                                                                                                                                                >  1:740536‑M1/1‑95 (MQ=255)                                                                                                                                                      gTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgt                                                                                                                                                                                                                                 >  1:193183‑M1/1‑94 (MQ=255)                                                                                                                                                                            ttgttCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaa                                                                         >  1:437030‑M1/1‑72 (MQ=255)                                                                                                                                                                              gttCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaata                                                                                                                                                                                                                                               >  1:587179‑M1/1‑70 (MQ=255)                                                                                                                                                                                         ggAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaatccaaacactgc                                                             <  1:152808‑M1/250‑192 (MQ=255)                                                                                                                                                                                            aGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatat                                                                                                                                                                                                                                                    >  1:781965‑M1/1‑56 (MQ=255)                                                                                                                                                                                             gATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatc                                                                                                                                            <  1:543309‑M1/167‑113 (MQ=255)                                                                                                                                                                                             gATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagt                                                                                                                                                                        >  1:299606‑M1/1‑55 (MQ=255)                                                                                                                                                                                              aTAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaatccaaacactgctgccaa                                                       >  1:816438‑M1/1‑54 (MQ=255)                                                                                                                                                                                                  tGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtt                                                                                        >  1:1020‑M1/1‑50 (MQ=255)                                                                                                                                                                                                                                      tCTTTAATTAGCTCaataatattaatattatcattttgtgtacatt                                                                                                                                                                                                                            <  1:536106‑M1/46‑33 (MQ=255)                                                                                                                                                                                                                                                gCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaatccaaacactgctgccaataaaaacgtcattaattttttccatttgccattggcgcttttcaatcac      >  1:804040‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                gCTCaataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattg                                                                                          >  1:501848‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                   caataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaatccaaacactgctgccaataaaaacgtcattaattttttccatttgccattggcgcttttcaatcaccttt  >  1:265934‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                                   caataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagattcgaattggccagcttctgtaattgtttttttaactcgttaatccaaacactgctgccaataaaaa                                                 >  1:321912‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                                   caataatattaatattatcattttgtgtacattcatgaaaaaccatcggtagtgtttttcgatcgcttaatttatgcagataagtagccgctaaatgacgaccggtgccaatctctccccaaataaaaaccggtagat                                                                                                                   >  1:644087‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           ATCAGTGATATGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTCGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCTTTTTCCG  >  1/1568197‑1568689

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