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	36 [0]	[0] 35	[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)	30 (0.660)	29/268	0.6	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

AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909169‑3909452                                                                                                                                               |                                                                                                                                              aaaGTAGCAGCCCGCTGCCCGCCAGAGCGGTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAG                                                                                                                                               <  2:286023‑M1/143‑1 (MQ=255)  aaTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGg                                                                                                                                              >  1:426694‑M1/1‑142 (MQ=255)   aTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgc                                                                                                                                             <  2:371994‑M1/143‑3 (MQ=255)          gCCCGCTACCCGCCAGTGCTTTCGCCAGTCCCCCACCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCTATACCGCGAGTGATGGTGTACATTAAAATCCAGATCAGCCAGgcgcaagct                                                                                                                                      <  2:155434‑M1/143‑10 (MQ=255)             cGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACGGCGAGCGATGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagcgcag                                                                                                                                   >  2:218705‑M1/1‑131 (MQ=255)                   ccGCCAGAGCGTTCGCCAGACCACCACCTTCCGTCTTGGGCGGCGGCGCCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcacgcgcagtcacga                                                                                                                             >  2:192678‑M1/1‑125 (MQ=255)                    cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCAAGATCAGCCAGacgcaagctcagtcaggaa                                                                                                                            >  2:544857‑M1/1‑124 (MQ=255)                              ttCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAGCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACCTTAAAATCCAGATCAGCCAGgctcaagctcagtcaggaacagggcgata                                                                                                                  >  1:460909‑M1/1‑114 (MQ=255)                                   cAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGCACAGCGCCAGCGACATACGGTCGATACCGCGAGTGATGGTGGACATTAAGATCCAGATCAGCCAGgcgagagatcagtccggagcagggacatgccgaa                                                                                                             >  2:246196‑M1/1‑109 (MQ=255)                                          ccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacc                                                                                                      <  1:351291‑M1/143‑42 (MQ=255)                                                              cggctTCATTTCACTGAAACGGGCCAGCGACATACCGTCTATACCGCGAGTGATGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgc                                                                                  <  2:59259‑M1/143‑62 (MQ=255)                                                                  gTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgatcacgaaactcaccgggacggcgatcagcagcgcgata                                                                              >  2:29419‑M1/1‑78 (MQ=255)                                                                            tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatacacgaaggca                                                                    >  1:485751‑M1/1‑68 (MQ=255)                                                                              aaCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGCTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgctaaacgaagtcacc                                                                  >  2:703554‑M1/1‑66 (MQ=255)                                                                                 aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTCAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgcgacggcgatcagcagcgcgatacacgaagtcaccaac                                                               >  2:67313‑M1/1‑63 (MQ=255)                                                                                  gcgcCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacg                                                              >  2:451694‑M1/1‑62 (MQ=255)                                                                                          gACACACCGTCGATACCGCGAGTGAGGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtag                                                      <  1:53806‑M1/143‑90 (MQ=255)                                                                                             aTACCGTCGATCCCGCGAGTCATCGTGCACATCAAAATCCAGTTTTGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatc                                                   <  1:30115‑M1/143‑93 (MQ=255)                                                                                                  gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcac                                              >  1:2916‑M1/1‑46 (MQ=255)                                                                                                      aTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgaccagcagcgcgataaacgcagaccccaacgtaccgcacaacggcaccagc                                          >  2:246177‑M1/1‑42 (MQ=255)                                                                                                          cgcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccgacgtaccgtagatcggcaccagcgccc                                      >  2:221364‑M1/1‑38 (MQ=255)                                                                                                            cgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcggcccg                                    >  1:264819‑M1/1‑36 (MQ=255)                                                                                                             cagcgatggTGCACATCAAAACCCGGACCAGCCCGgcgcaagctcagccaggaacagggcgataccgatcctcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgt                                   <  2:193072‑M1/139‑109 (MQ=255)                                                                                                                tGATGGTGGACGTGAAAATCCAGAGCGGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgagcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                                <  2:798515‑M1/143‑112 (MQ=255)                                                                                                                     gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                           <  1:719922‑M1/143‑117 (MQ=255)                                                                                                                       ggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagtgcgataccgcaactcaccaagacagcaaacagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcacaagcgacccgtagatatcgtt                         >  1:53882‑M1/1‑25 (MQ=255)                                                                                                                         aCATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcg                       >  1:261887‑M1/1‑23 (MQ=255)                                                                                                                           aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                     >  1:635918‑M1/1‑21 (MQ=255)                                                                                                                                   ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaccgtacaatagaacggcaccagagccccgtagataacgcccgatgcagccca             >  2:374556‑M1/1‑13 (MQ=255)                                                                                                                                     aGATCAGCCGGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccact           <  2:347780‑M1/143‑133 (MQ=255)                                                                                                                                     aGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccact           <  1:680798‑M1/143‑133 (MQ=255)                                                                                                                                          gggcaggcgcaagctcagtcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  1:726569‑M1/140‑138 (MQ=255)                                                                                                                                          aGCCAGgcgcggcctcagccgcgagcagggcgataccgaaactccccgggacggcgatcagcagcgcgataaactaagtcaccaacgcaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  2:489884‑M1/143‑138 (MQ=255)                                                                                                                                            ccAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagtagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccggagatatcgttcggtgcatcccacgctttggt    >  1:722244‑M1/1‑4 (MQ=255)                                                                                                                                             cAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccacgctttggtc   >  1:510154‑M1/1‑3 (MQ=255)                                                                                                                                              aGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtcc  >  1:359375‑M1/1‑2 (MQ=255)                                                                                                                                               |                                                                                                                                              AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909169‑3909452

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