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	44 [0]	[0] 44	[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)	40 (0.880)	36/268	0.2	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

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

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