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

AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909168‑3909450                                                                                                                                               |                                                                                                                                             acattAACACCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAG                                                                                                                                              <  1:248133‑M1/141‑1 (MQ=255)       ccctcccGCTTCCTTCCAGAGCGTCCGCCCGACCACCTCCTCCCGTAATAGGCGGCTGCGTCATTTCACTGCACAGTGCCAGCGACATACCGTCGATACCCCGAGTGATGGTGCACATTAAAATCCAGTTCAGCCAggcgcaa                                                                                                                                        <  1:75107/139‑7 (MQ=255)        cAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAAGCGAGCTCAGCCAGgcgcaag                                                                                                                                       >  1:872884‑M1/1‑136 (MQ=255)          gCCCTCTGCGCGCCAGGGCGTTCGCCAGACCACCACCTTCCGTATTGAGCGGCGGCTTCATTTCAGTGAGGAGCGCCAGCGACGTGTCGTCGACACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagct                                                                                                                                     <  1:810586‑M1/143‑10 (MQ=255)           cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCACgcgcaagctc                                                                                                                                    >  2:693273‑M1/1‑133 (MQ=255)            ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTCTTGGGCGGCGGCGTCACTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGAAGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctca                                                                                                                                   >  1:195743‑M1/1‑132 (MQ=255)            ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCCGgcgcaagctcg                                                                                                                                   >  1:194787‑M1/1‑132 (MQ=255)                   ccGCCATAGCGTTCGCCGGACCACCCCCTTCCGTATTGGGCGGCGGCGTCAGTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcagga                                                                                                                            <  2:254219‑M1/143‑19 (MQ=255)                     gCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTCACATCCAGATCAGCCGGgcgcaggcccggtcagtaac                                                                                                                          >  2:780978‑M1/1‑123 (MQ=255)                        agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctccgtcaggagcagg                                                                                                                       >  2:577802‑M1/1‑120 (MQ=255)                        agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATACAGATCAGCCAGgcgcaagctcaggcaggaacaag                                                                                                                       >  1:337308‑M1/1‑120 (MQ=255)                                 gCCAGACCACCACCTGCGGTAGTGGGCGGCGGCGTCATTGCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                              <  2:312515‑M1/143‑33 (MQ=255)                                      accaccaccTTCCGTATTGGGCGGCGGCGTCCTTACAGTGAACAGCGCCCGCGACATACCGTCGATACCGCGAGTGAGGGTGGACAGTAACATCCAGATCAGCCCGccgcaagctcagtcaggcaaacggcaacaccgaacga                                                                                                         >  1:42827‑M1/1‑106 (MQ=255)                                      accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaact                                                                                                         <  1:337156‑M1/143‑38 (MQ=255)                                        caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactca                                                                                                       <  2:577333‑M1/143‑40 (MQ=255)                                        caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactca                                                                                                       >  2:32088‑M1/1‑104 (MQ=255)                                         accaccTTCCGTATCGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:359371‑M1/143‑41 (MQ=255)                                         accaccGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:656171‑M1/143‑41 (MQ=255)                                                     cttgGGCGGCGGCGTCATTTGAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  1:150179‑M1/142‑53 (MQ=255)                                                            ggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTACAATCCAGATCACCCAGgcgcaagctcagccagcaccagcgccataccgaacctccccgggcccgcgatcacccgc                                                                                   >  1:204804‑M1/1‑84 (MQ=255)                                                              cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaacacaccgggacggcgatcagcagcgc                                                                                 >  1:514018‑M1/1‑82 (MQ=255)                                                              cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGAGCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgc                                                                                 <  2:416319‑M1/143‑62 (MQ=255)                                                               ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGGGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                <  2:860405‑M1/143‑63 (MQ=255)                                                                gcgTCACTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcga                                                                               <  1:218328‑M1/143‑64 (MQ=255)                                                                   tCATTTCAGTGAACAGCGCCAGCCACACACCTTTGCTACCGCGAGTGATGGTGGACATTAAAGTCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataa                                                                            <  2:594210‑M1/143‑67 (MQ=255)                                                                    gacttcagTGCACAGCGCCAGCGACATACCGTCGATACCGCGAGTCATGGTGGACATGAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaa                                                                           <  1:688166‑M1/140‑68 (MQ=255)                                                                        tCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  2:2266‑M1/143‑72 (MQ=255)                                                                            tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   >  2:286708‑M1/1‑68 (MQ=255)                                                                            tGAACAGCGCCAGCGACAGCCCGTCGATACCGCTAGTGATGGTGTACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:47542‑M1/143‑76 (MQ=255)                                                                                 aGCGCCAGCGACATCCCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                              <  2:186780‑M1/143‑81 (MQ=255)                                                                                  gcgcCAGCGACATCCCGGCGATACCGGGAGTGATGGTGGACATTAAAATCCAGATCAGCCGGgcgcgaggtcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacg                                                             <  2:405481‑M1/143‑82 (MQ=255)                                                                                     ccatccacttaccttcgatACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtac                                                          <  1:566677‑M1/134‑85 (MQ=255)                                                                                        gCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgt                                                       >  2:472198‑M1/1‑56 (MQ=255)                                                                                          gACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtag                                                     <  2:414313‑M1/143‑90 (MQ=255)                                                                                           aCATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtaga                                                    >  2:394493‑M1/1‑53 (MQ=255)                                                                                              tACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctaagtcaggaacagggcgataccgacactcaccgggacggcgatcagcagcgcgataaacgaagtccccaacgtaccgcaaatcc                                                 >  2:395304‑M1/1‑50 (MQ=255)                                                                                                  gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagataggcac                                             >  1:727138‑M1/1‑46 (MQ=255)                                                                                                           gcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatacacgaagtcacccacgtaccgtagttcgtcacccgcgcccc                                    >  2:344007‑M1/1‑37 (MQ=255)                                                                                                               gTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagatcggaagag                                <  1:833301‑M1/143‑111 (MQ=255)                                                                                                                        gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactccccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtggatcggcacgagcgccccgtaggtgtcggtc                       >  2:178679‑M1/1‑24 (MQ=255)                                                                                                                        gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaaggcaccaacgtaccgtagatcggcaccagcgccccggagatatcgttc                       >  2:178678‑M1/1‑24 (MQ=255)                                                                                                                           aTTAAAGTCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatccgcaacgcgataaacgaaataacacacgtaccgcagatcagcaccagcgcccagtaggactcgtagcga                    >  1:292383‑M1/1‑21 (MQ=255)                                                                                                                           aTTAAAATCCAGATCCGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    <  1:245648‑M1/143‑123 (MQ=255)                                                                                                                           aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    >  1:450416‑M1/1‑21 (MQ=255)                                                                                                                            tcaaaatcCAGATCACCCAGgcgcaagctcagtcaggaacagggcgataccgacactccccgtgacggcgatcagcagcgcgatcaacgtagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg                   <  2:109‑M1/141‑124 (MQ=255)                                                                                                                                            ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggt   <  2:727138‑M1/143‑140 (MQ=255)                                                                                                                                            ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaacagcgccccgtagatatcgttcggtgcatcccactctttagc   >  1:475177‑M1/1‑4 (MQ=255)                                                                                                                                             cAGgcacaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtc  <  2:425111‑M1/143‑141 (MQ=255)                                                                                                                                               |                                                                                                                                             AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909168‑3909450

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