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	33 [0]	[0] 32	[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)	32 (0.680)	29/270	0.7	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

TGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGC  >  NC_000913/3909954‑3910209                                                                                                                                     |                                                                                                                            tcgctacccgccggggcgttcgccagacctccaccctccgtatggggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtca                                                                                                                   <  2:18522‑M2/11‑1 (MQ=255)     tacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGctcagtcg                                                                                                               >  1:709390‑M2/129‑142 (MQ=255)          gccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAAc                                                                                                          >  1:136225‑M2/124‑143 (MQ=255)           ccagggcgttcgccagcccaccaccttccgtatcgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACa                                                                                                         <  2:571108‑M2/21‑1 (MQ=255)               agcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtagataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGc                                                                                                     >  1:177989‑M2/119‑143 (MQ=255)                 cgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGa                                                                                                   <  1:749181‑M2/27‑1 (MQ=255)                        cagcccaccatcctccgtattgggcggcggcgtcacggcagtgaacagcgccagcgacctaccgtcgataccgcgagtgatcgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGACGATACCGaa                                                                                            <  2:448651‑M2/34‑1 (MQ=255)                        cagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                      >  2:62032‑M2/110‑133 (MQ=255)                        cagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                      <  1:62032‑M2/24‑1 (MQ=255)                          cagcaccgccttccgtatggggcggcgccgccacttcagtgaacagcgccagcgacatcccgtggataccgcgagtgatggtggacattaaaatccatatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGacac                                                                                          <  1:366892‑M2/36‑1 (MQ=255)                                   cttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGa                                                                                 >  1:536161‑M2/99‑143 (MQ=255)                                         tattgcgcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgattgtgtacattataatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                           <  1:447853‑M2/51‑1 (MQ=255)                                         gagtgggcggcgacggaagttcagtgaacggcgccagcgacataccgtcgacaccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                           <  1:878495‑M2/51‑1 (MQ=255)                                                cggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcagcag                                                                    >  1:297526‑M2/86‑143 (MQ=255)                                                  gcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcg                                                                  <  1:721482‑M2/60‑1 (MQ=255)                                                             tcagtggactgcgccagagagatagggtcgataccgagagtgattgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                       <  2:318403‑M2/71‑1 (MQ=255)                                                                gtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATGAGCAGCGCGACAAACGAAGTc                                                    >  2:786944‑M2/70‑143 (MQ=255)                                                                   aacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCAcc                                                 <  1:495887‑M2/77‑1 (MQ=255)                                                                     gcacgccagcgccttagcgtcgctaccgggggtgatggtggacattaaaatccagatcagccagGCGCAAGGTCAGTCAGGAACAGGGCTATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCaa                                               <  2:762033‑M2/79‑1 (MQ=255)                                                                       gcgccagcgacagaccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                             <  2:34100‑M2/81‑1 (MQ=255)                                                                           cagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAcc                                         >  1:831780‑M2/59‑143 (MQ=255)                                                                           aagcgacataccgtcgatacagcgagggatggtggacagtaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAcc                                         <  1:810648‑M2/85‑1 (MQ=255)                                                                              cgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACGGGGACGGCGATCAGCGGCGAGACAAAAGAAGACACCAACGTACGGTa                                      >  2:366751‑M2/56‑143 (MQ=255)                                                                                  ataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACCGGTCGGTACCGAAACTCACCGGGACGACGATCAGCAGAGGGCTCAACGACGTCACCAACGTACCGTagagc                                  >  1:563231‑M2/52‑141 (MQ=255)                                                                                   gaccgtcgagaccgcgagtgatggtggacaagaaaatcgagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCg                                 <  1:294649‑M2/93‑1 (MQ=255)                                                                                        tcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGAGCGGCAcc                            >  1:542569‑M2/46‑143 (MQ=255)                                                                                           ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgc                         >  1:534311‑M2/43‑143 (MQ=255)                                                                                           ataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAACAGCGCGACAAACGAAGTCACCAACGTACGCTAGCTCGGCACCAgc                         >  2:416693‑M2/43‑143 (MQ=255)                                                                                                          gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg          >  2:719445‑M2/28‑143 (MQ=255)                                                                                                          gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGACAAACGCAGTCACCAACGTACCGTAGCTCGGCACCAGCGCCCCGGAGATATCg          >  2:719451‑M2/28‑143 (MQ=255)                                                                                                                attaaactccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGAGCAGGCTGTGTTCAGCGCGAGAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt    <  1:903679‑M2/122‑1 (MQ=255)                                                                                                                  taaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGACACCAGCGCCCCGTAGCTATCGTGCGgcgc  >  1:587969‑M2/20‑143 (MQ=255)                                                                                                                                     |                                                                                                                            TGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGC  >  NC_000913/3909954‑3910209

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: ‑