Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,909,312	Δ774 bp	100%		[pstA]–[pstC]	[pstA], [pstC]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3909312	3910085	774	34 [0]	[0] 34	[pstA]–[pstC]	[pstA],[pstC]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3909311	0 (0.000)	30 (0.910)	28/484	NT	100%	coding (129/891 nt)	pstA	phosphate ABC transporter permease
	?	NC_000913	3910086 =	0 (0.000)					coding (313/960 nt)	pstC	phosphate ABC transporter permease

CATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTAC  >  NC_000913/3909799‑3910364                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                       acagaccaaccacaatcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggc                                                                                                                                                                                                                                                                                       >  1:27100‑M2/288‑289 (MQ=255)    gaccaaccacaatcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCa                                                                                                                                                                                                                                                                                    >  2:266206‑M2/285‑289 (MQ=255)      caaacaccaatggccggcgcagagagcggaaggtcgttaatgaagcgaatcactcctgccagccaggatttccgaccatattccgccagataaacccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAg                                                                                                                                                                                                                                                                                  <  2:61463‑M2/7‑1 (MQ=255)               atcgacggcgcagagagcagaatgtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtcagtaag                                                                                                                                                                                                                                                                         >  1:317866‑M2/274‑289 (MQ=255)                        caagaaggcagaaggtcgtaagtgaagcgaatcccttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                                                                                                                                                                                 <  2:317866‑M2/24‑1 (MQ=255)                                      gtcgttaatgaagcgaatcacttctgccagccaggatttacgaccatattccgccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTc                                                                                                                                                                                                                                                   >  2:304820‑M2/251‑288 (MQ=255)                                                                                           ccagataaatccccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATc                                                                                                                                                                                             >  1:165986‑M2/198‑289 (MQ=255)                                                                                                      cccgccataatgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgcgc                                                                                                                                                                                  <  1:266206‑M2/103‑1 (MQ=255)                                                                                                                tgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                                                                                                                                                                               >  2:93433‑M2/177‑234 (MQ=255)                                                                                                                tgcccagcggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                                                                                                                                                                               <  1:93433‑M2/58‑1 (MQ=255)                                                                                                                       cggcgtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                                                                                                                                                                 >  2:43365‑M2/170‑289 (MQ=255)                                                                                                                           gtaccgaataccgtggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                                                                                                                                                              <  1:275053‑M2/123‑1 (MQ=255)                                                                                                                           gcaccaaaaaccgtggcccacaaaattaacagcccccgccccgccaaagcgttcgccagaccaccaccttccgtattgggcgtcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacagtaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                                                                                                                                                              <  2:165986‑M2/123‑1 (MQ=255)                                                                                                                             accgaataccgtcgcccacaaaataaccagccccctacacgccagggcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCa                                                                                                                                                            <  2:27100‑M2/125‑1 (MQ=255)                                                                                                                                        tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                                                                                                                                                               <  1:21935‑M2/122‑1 (MQ=255)                                                                                                                                        tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                                                                                                                                                               >  2:21935‑M2/153‑274 (MQ=255)                                                                                                                                        tggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGACATCGTTCGGTGCATCCCACTCTTTg                                                                                                                                                >  2:318121‑M2/153‑289 (MQ=255)                                                                                                                                         ggcccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggccattaaaatccaaatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCCCCAAAGTACCGTAGCTCGGCAACAGCGCCCAGTAGATATCGTTCGGTGCATCCCACTCTTTgg                                                                                                                                               >  2:213695‑M2/152‑289 (MQ=255)                                                                                                                                           cccacaaaattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaccagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATGAGCAGCGCGATAAACTAAGTCACCAACGCACCGTAGATCGCCACCAGCGCCACGTAGATATCGATCGGTGCATCCACCTCTTTGGt                                                                                                                                              >  1:298530‑M2/150‑288 (MQ=255)                                                                                                                                                          gcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAACCGAAGt                                                                                                                                                                                                                >  1:276534‑M2/135‑207 (MQ=255)                                                                                                                                                          gcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGt                                                                                                                                                                                                                <  2:276534‑M2/73‑1 (MQ=255)                                                                                                                                                             cgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGAACCAGCGCCCCGTAGATATCGTTCGGTTCATCCCACTCTTTGGTCCATAGGAAAGCACGACCa                                                                                                                           >  2:89616‑M2/132‑289 (MQ=255)                                                                                                                                                                          agcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATg                                                                                                               >  2:182459‑M2/119‑288 (MQ=255)                                                                                                                                                                                     taccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCaggag                                                                                                   <  1:304820‑M2/182‑1 (MQ=255)                                                                                                                                                                                       ccaccacctgccgtattgggcggcggcgtcatgtcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAggagga                                                                                                 <  1:89616‑M2/184‑1 (MQ=255)                                                                                                                                                                                                         ggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACaat                                                                               >  2:83483‑M2/88‑289 (MQ=255)                                                                                                                                                                                                              cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCag                                                                                        >  2:71859‑M2/83‑275 (MQ=255)                                                                                                                                                                                                              cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCag                                                                                        <  1:71859‑M2/193‑1 (MQ=255)                                                                                                                                                                                                                                        gcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAgc                                                >  1:6826‑M2/57‑289 (MQ=255)                                                                                                                                                                                                                                               accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCa                                         <  1:43365‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                                                 gatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCt                         <  1:182459‑M2/256‑1 (MQ=255)                                                                                                                                                                                                                                                                                  tccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACACTAATGCCACCCAACATCAATAGCACAATCAGCCCCGCCAGTTTTACCAGCCAGCTGAACATTCTGTCGACCttt      >  2:26743‑M2/15‑289 (MQ=255)                                                                                                                                                                                                                                                                                       atcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCAGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCCGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCCCCGCCAGTTTTACCAGCACGCTGAAAATTAAGTCCCCATTTTTAc  >  2:109736‑M2/10‑288 (MQ=255)                                                                                                                                                                                                                                                                                       atcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCACAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACACAACATCAATAGCCCAATCAGCGCCGCAAGTTTTACCAGAACGCTGAAAATAATGTAGCCCTTTTTAc  >  1:244222‑M2/10‑288 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                       CATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGAACAGCGCGCCAACAATCGGGATATTCGACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCATAGGAAAGCCAGACCAAATTTCTGAATGCTCGGCCAGGAGGAGATGATCAGAGAGACAATAATGCCACCCAACATCAATAGCACAATCAGCGCCGCCAGTTTTACCAGCACGCTGAAAATTATGTCGCCCTTTTTAC  >  NC_000913/3909799‑3910364

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