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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3909311	3910084	774	48 [0]	[0] 47	[pstA]–[pstC]	[pstA], [pstC]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3909310	0 (0.000)	43 (0.900)	37/268	0.2	100%	coding (129/891 nt)	pstA	phosphate ABC transporter permease
	?	NC_000913	3910085 =	0 (0.000)					coding (313/960 nt)	pstC	phosphate ABC transporter permease

ATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909949‑3910227                                                                                                                                         |                                                                                                                                               cagcccgctacccgccagagcgttcgccagaccaccaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaagcgagctcagccagGCGCAAg                                                                                                                                          >  1:872884‑M2/137‑143 (MQ=255)   gccctctgcgcgccagggcgttcgccagaccaccaccttccgtattgagcggcggcttcatttcagtgaggagcgccagcgacgtgtcgtcgacaccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCt                                                                                                                                        <  1:810586‑M2/9‑1 (MQ=255)    cccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccacGCGCAAGCtc                                                                                                                                       >  2:693273‑M2/134‑143 (MQ=255)     ccgctacccgccagagcgttcgccagaccaccaccttccgtcttgggcggcggcgtcacttcagtgaacagcgccagcgacataccgtcgataccgcgagtgaaggtggacattaaaatccagatcagccagGCGCAAGCtca                                                                                                                                      >  1:195743‑M2/133‑143 (MQ=255)     ccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagcccgGCGCAAGctcg                                                                                                                                      >  1:194787‑M2/133‑142 (MQ=255)            ccgccatagcgttcgccggaccacccccttccgtattgggcggcggcgtcagttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGa                                                                                                                               <  2:254219‑M2/18‑1 (MQ=255)              gccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattcacatccagatcagccggGCGCAGGCCCGGTCAGTAAc                                                                                                                             >  2:780978‑M2/124‑143 (MQ=255)                 agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCCGTCAGGAGCAgg                                                                                                                          >  2:577802‑M2/121‑143 (MQ=255)                 agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatacagatcagccagGCGCAAGCTCAGGCAGGAacaag                                                                                                                          >  1:337308‑M2/121‑141 (MQ=255)                          gccagaccaccacctgcggtagtgggcggcggcgtcattgcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                                 <  2:312515‑M2/32‑1 (MQ=255)                               accaccaccttccgtattgggcggcggcgtccttacagtgaacagcgcccgcgacataccgtcgataccgcgagtgagggtggacagtaacatccagatcagcccgCCGCAAGCTCAGTCAGGCAAACGGCAACAccgaacga                                                                                                            >  1:42827‑M2/107‑140 (MQ=255)                               accaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACt                                                                                                            <  1:337156‑M2/37‑1 (MQ=255)                                 caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                          <  2:577333‑M2/39‑1 (MQ=255)                                 caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                          >  2:32088‑M2/105‑143 (MQ=255)                                  accaccttccgtatcgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                         <  1:359371‑M2/40‑1 (MQ=255)                                  accaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                         <  1:656171‑M2/40‑1 (MQ=255)                                              cttgggcggcggcgtcatttgagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                             <  1:150179‑M2/52‑1 (MQ=255)                                                     ggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattacaatccagatcacccagGCGCAAGCTCAGCCAGCACCAGCGCCATACCGAACCTCCCCGGGCCCgcgatcacccgc                                                                                      >  1:204804‑M2/85‑138 (MQ=255)                                                       cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACACACCGGGACGGCGATCAGCAgcgc                                                                                    >  1:514018‑M2/83‑143 (MQ=255)                                                       cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagagcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcgc                                                                                    <  2:416319‑M2/61‑1 (MQ=255)                                                        ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatgggggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                   <  2:860405‑M2/62‑1 (MQ=255)                                                         gcgtcacttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGa                                                                                  <  1:218328‑M2/63‑1 (MQ=255)                                                            tcatttcagtgaacagcgccagccacacacctttgctaccgcgagtgatggtggacattaaagtccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATaa                                                                               <  2:594210‑M2/66‑1 (MQ=255)                                                             gacttcagtgcacagcgccagcgacataccgtcgataccgcgagtcatggtggacatgaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATaaa                                                                              <  1:688166‑M2/67‑1 (MQ=255)                                                                 tcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                          <  2:2266‑M2/71‑1 (MQ=255)                                                                     tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                      >  2:286708‑M2/69‑143 (MQ=255)                                                                     tgaacagcgccagcgacagcccgtcgataccgctagtgatggtgtacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                      <  1:47542‑M2/75‑1 (MQ=255)                                                                          agcgccagcgacatcccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAAc                                                                 <  2:186780‑M2/80‑1 (MQ=255)                                                                           gcgccagcgacatcccggcgataccgggagtgatggtggacattaaaatccagatcagccggGCGCGAGGTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACg                                                                <  2:405481‑M2/81‑1 (MQ=255)                                                                              ccatccacttaccttcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAc                                                             <  1:566677‑M2/84‑1 (MQ=255)                                                                                 gcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGt                                                          >  2:472198‑M2/57‑143 (MQ=255)                                                                                   gacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAg                                                        <  2:414313‑M2/89‑1 (MQ=255)                                                                                    acataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGa                                                       >  2:394493‑M2/54‑143 (MQ=255)                                                                                       taccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTAAGTCAGGAACAGGGCGATACCGACACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCCCCAACGTACCGCAAAtcc                                                    >  2:395304‑M2/51‑142 (MQ=255)                                                                                           gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATAGGCAc                                                >  1:727138‑M2/47‑143 (MQ=255)                                                                                                    gcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATACACGAAGTCACCCACGTACCGTAGTTCGTCACCCGCGcccc                                       >  2:344007‑M2/38‑143 (MQ=255)                                                                                                        gtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGatcggcaccagatcggaagag                                   <  1:833301‑M2/110‑11 (MQ=255)                                                                                                                 gacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCCCCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTGGATCGGCACGAGCGCCCCGTAGGTGTCGGTc                          >  2:178679‑M2/25‑143 (MQ=255)                                                                                                                 gacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGGCACCAACGTACCGTAGATCGGCACCAGCGCCCCGGAGATATCGTTc                          >  2:178678‑M2/25‑143 (MQ=255)                                                                                                                    attaaagtccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCCGCAACGCGATAAACGAAATAACACACGTACCGCAGATCAGCACcagcgcccagtaggactcgtagcga                       >  1:292383‑M2/22‑131 (MQ=255)                                                                                                                    attaaaatccagatccgccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                       <  1:245648‑M2/122‑1 (MQ=255)                                                                                                                    attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                       >  1:450416‑M2/22‑143 (MQ=255)                                                                                                                     tcaaaatccagatcacccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGACACTCCCCGTGACGGCGATCAGCAGCGCGATCAACGTAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTg                      <  2:109‑M2/123‑1 (MQ=255)                                                                                                                                     ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGt      <  2:727138‑M2/139‑1 (MQ=255)                                                                                                                                     ccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAACAGCGCCCCGTAGATATCGTTCGGTGCATCCCActctttagc      >  1:475177‑M2/5‑140 (MQ=255)                                                                                                                                      cagGCACAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTc     <  2:425111‑M2/140‑1 (MQ=255)                                                                                                                                         gcgcAAGCTCAGTCAGGACCAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCGCCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAt  >  2:32404/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                               ATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCCAT  >  NC_000913/3909949‑3910227

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