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

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

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