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

GCCACATCAGTGATATGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTCGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCTTTTTCCG  >  1/1568192‑1568689                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          atcagtgatatgtttacgtaacggcagcatttcaatttgagtataaataaaaaaataataaaactcggcacttatctgccgacagcgaattaattcgttaagttctatatcactaagcaaaattaaacggaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataatat                                                                                                                                                                                                                                                         <  1:789096‑M2/3‑1 (MQ=255)                                                                       cttatctgccgacagcgaattaattcgttaagttctatatcactaagcaaaattaaacggaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAg                                                                                                                                                                                   >  1:814587‑M2/179‑251 (MQ=255)                                                                         tatctgccgacagcgaattaattcgttaagttctatatcactaagcaaaattaaacggaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGa                                                                                                                                                                                                    <  1:88137‑M2/56‑1 (MQ=255)                                                                                                    aagttctatatcactaagcaaaattaaacggaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTg                                                                                                                                                       >  1:253250‑M2/150‑250 (MQ=255)                                                                                                         ctatatcactaagcaaaattaaacggaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAAtc                                                                                                                                                 <  1:705590‑M2/107‑1 (MQ=255)                                                                                                                                  gaattgactctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAAc                                                                                                                                <  1:149387‑M2/124‑1 (MQ=255)                                                                                                                                           ctcgtcatgatgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCgg                                                                                                                                                                                                                 >  1:277990‑M2/111‑153 (MQ=255)                                                                                                                                                     tgtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGc                                                                                                                                                                                     >  1:740536‑M2/101‑171 (MQ=255)                                                                                                                                                      gtaaaatctgaattaaatattgttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTgtgt                                                                                                                                                                                                                                      >  1:193183‑M2/100‑121 (MQ=255)                                                                                                                                                                            ttgttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTaa                                                                              >  1:437030‑M2/78‑251 (MQ=255)                                                                                                                                                                              gttcttcaacaggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataatattaata                                                                                                                                                                                                                                                    >  1:587179‑M2/76‑83 (MQ=255)                                                                                                                                                                                         ggaagataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACActgc                                                                  <  1:152808‑M2/186‑1 (MQ=255)                                                                                                                                                                                            agataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataatat                                                                                                                                                                                                                                                         >  1:781965‑M2/62‑64 (MQ=255)                                                                                                                                                                                             gataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGt                                                                                                                                                                             >  1:299606‑M2/61‑139 (MQ=255)                                                                                                                                                                                             gataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAAtc                                                                                                                                                 <  1:543309‑M2/107‑1 (MQ=255)                                                                                                                                                                                              ataatgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCaa                                                            >  1:816438‑M2/60‑251 (MQ=255)                                                                                                                                                                                                  tgcacattcttaataattaaagtccctgattttgtttctttaattagctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGtt                                                                                             >  1:1020‑M2/56‑214 (MQ=255)                                                                                                                                                                                                                                      tctttaattagctcaataaTATTAATATTATCATTTTGTGTACAtt                                                                                                                                                                                                                                 <  1:536106‑M2/27‑1 (MQ=255)                                                                                                                                                                                                                                                gctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTg                                                                                               >  1:501848‑M2/10‑166 (MQ=255)                                                                                                                                                                                                                                                gctcaataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCAc           >  1:804040‑M2/10‑250 (MQ=255)                                                                                                                                                                                                                                                   caataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGAt                                                                                                                        >  1:644087‑M2/7‑138 (MQ=255)                                                                                                                                                                                                                                                   caataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATaaaaa                                                      >  1:321912‑M2/7‑204 (MQ=255)                                                                                                                                                                                                                                                   caataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCttt       >  1:265934‑M2/7‑251 (MQ=255)                                                                                                                                                                                                                                                     ataaTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGgcgc                      <  1:549702‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                                        aTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTaa                                                                                                  >  1:736996‑M2/2‑155 (MQ=255)                                                                                                                                                                                                                                                        aTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCTTTTTCCg  >  1:4846‑M2/2‑251 (MQ=255)                                                                                                                                                                                                                                                         tATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTAGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGt                                                   >  1:674002/1‑201 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          GCCACATCAGTGATATGTTTACGTAACGGCAGCATTTCAATTTGAGTATAAATAAAAAAATAATAAAACTCGGCACTTATCTGCCGACAGCGAATTAATTCGTTAAGTTCTATATCACTAAGCAAAATTAAACGGAATTGACTCTCGTCATGATGTAAAATCTGAATTAAATATTGTTGTTCTTCAACAGGAAGATAATGCACATTCTTAATAATTAAAGTCCCTGATTTTGTTTCTTTAATTAGCTCTATTAATATTATCATTTTGTGTACATTCATGAAAAACCATCGGTAGTGTTTTTCGATCGCTTAATTTATGCAGATAAGTCGCCGCTAAATGACGACCGGTGCCAATCTCTCCCCAAATAAAAACCGGTAGATTCGAATTGGCCAGCTTCTGTAATTGTTTTTTTAACTCGTTAATCCAAACACTGCTGCCAATAAAAACGTCATTAATTTTTTCCATTTGCCATTGGCGCTTTTCAATCACCTTTTTCCG  >  1/1568192‑1568689

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