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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3909312	3910085	774	36 [0]	[0] 35	[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.660)	29/268	0.6	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

GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909944‑3910226                                                                                                                                               |                                                                                                                                             aattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagg                                                                                                                                              >  1:426694‑M2/143‑143 (MQ=255)  attaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggc                                                                                                                                             <  2:371994‑M2/2‑1 (MQ=255)         gcccgctacccgccagtgctttcgccagtcccccaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtctataccgcgagtgatggtgtacattaaaatccagatcagccagGCGCAAGCt                                                                                                                                      <  2:155434‑M2/9‑1 (MQ=255)            cgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgatacggcgagcgatggtggacagtaaaatccagatcagccagGCGCAAGCgcag                                                                                                                                   >  2:218705‑M2/132‑143 (MQ=255)                  ccgccagagcgttcgccagaccaccaccttccgtcttgggcggcggcgccatttcagtgaacagcgccagcgacataccgtcgataccgcgcgtgatggtggacattaaaatccagatcagccagGCGCACGCGCAGTCACGa                                                                                                                             >  2:192678‑M2/126‑143 (MQ=255)                   cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatcaagatcagccagACGCAAGCTCAGTCAGGaa                                                                                                                            >  2:544857‑M2/125‑143 (MQ=255)                             ttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgagcagcgccagcgacataccgtcgataccgcgagtgatggtggaccttaaaatccagatcagccagGCTCAAGCTCAGTCAGGAACAGGGCGATa                                                                                                                  >  1:460909‑M2/115‑143 (MQ=255)                                  cagaccaccaccttccgtattgggcggcggcgtcatttcagtgcacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaagatccagatcagccagGCGAGAGATCAGTCCGGAGCAGGGACATGCCGaa                                                                                                             >  2:246196‑M2/110‑143 (MQ=255)                                         ccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAcc                                                                                                      <  1:351291‑M2/41‑1 (MQ=255)                                                             cggcttcatttcactgaaacgggccagcgacataccgtctataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcgc                                                                                  <  2:59259‑M2/61‑1 (MQ=255)                                                                 gtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATCACGAAACTCACCGGGACGGCGATCAGCAGCGCGATa                                                                              >  2:29419‑M2/79‑143 (MQ=255)                                                                           tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATACACGAAGGCa                                                                    >  1:485751‑M2/69‑143 (MQ=255)                                                                             aacagcgccagcgacataccgtcgataccgcgagtgatgctggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGCTAAACGAAGTCAcc                                                                  >  2:703554‑M2/67‑143 (MQ=255)                                                                                agcgccagcgacataccgtcgataccgcgagtgatggtggacattcaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGCGACGGCGATCAGCAGCGCGATACACGAAGTCACCAAc                                                               >  2:67313‑M2/64‑143 (MQ=255)                                                                                 gcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                                              >  2:451694‑M2/63‑143 (MQ=255)                                                                                         gacacaccgtcgataccgcgagtgagggtggacagtaaaatccagatcagccagGCGCAAGCGCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAg                                                      <  1:53806‑M2/89‑1 (MQ=255)                                                                                            ataccgtcgatcccgcgagtcatcgtgcacatcaaaatccagttttgccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATc                                                   <  1:30115‑M2/92‑1 (MQ=255)                                                                                                 gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAc                                              >  1:2916‑M2/47‑143 (MQ=255)                                                                                                     ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGACCAGCAGCGCGATAAACGCAGACCCCAACGTACCGCACAACGGCACCAgc                                          >  2:246177‑M2/43‑143 (MQ=255)                                                                                                         cgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCGACGTACCGTAGATCGGCACCAGCGccc                                      >  2:221364‑M2/39‑143 (MQ=255)                                                                                                           cgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGGCCCg                                    >  1:264819‑M2/37‑143 (MQ=255)                                                                                                            cagcgatggtgcacatcaaaacccggaccagcccgGCGCAAGCTCAGCCAGGAACAGGGCGATACCGATCCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGt                                   <  2:193072‑M2/108‑1 (MQ=255)                                                                                                               tgatggtggacgtgaaaatccagagcggccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAGCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                <  2:798515‑M2/111‑1 (MQ=255)                                                                                                                    gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           <  1:719922‑M2/116‑1 (MQ=255)                                                                                                                      ggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGTGCGATACCGCAACTCACCAAGACAGCAAACAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACAAGCGACCCGTAGATATCGtt                         >  1:53882‑M2/26‑143 (MQ=255)                                                                                                                        acattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                       >  1:261887‑M2/24‑143 (MQ=255)                                                                                                                          attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                     >  1:635918‑M2/22‑143 (MQ=255)                                                                                                                                  ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCACCGTACAATAGAACGGCACCAGAGCCCCGTAGATAACGCCCGATGCAGCCCa             >  2:374556‑M2/14‑143 (MQ=255)                                                                                                                                    agatcagccggGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAct           <  2:347780‑M2/132‑1 (MQ=255)                                                                                                                                    agatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAct           <  1:680798‑M2/132‑1 (MQ=255)                                                                                                                                         gggcagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  1:726569‑M2/137‑1 (MQ=255)                                                                                                                                         agccagGCGCGGCCTCAGCCGCGAGCAGGGCGATACCGAAACTCCCCGGGACGGCGATCAGCAGCGCGATAAACTAAGTCACCAACGCACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  2:489884‑M2/137‑1 (MQ=255)                                                                                                                                           ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGTAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGGAGATATCGTTCGGTGCATCCCACGCTTTGGt    >  1:722244‑M2/5‑143 (MQ=255)                                                                                                                                            cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACGCTTTGGTc   >  1:510154‑M2/4‑143 (MQ=255)                                                                                                                                             agGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTcc  >  1:359375‑M2/3‑143 (MQ=255)                                                                                                                                               |                                                                                                                                             GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909944‑3910226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑