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	50 [0]	[0] 50	[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)	48 (0.950)	42/268	0.1	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

AATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909169‑3909450                                                                                                                                              |                                                                                                                                             aaTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGACGGGCATCAAACTCCCGATCAGCCAGg                                                                                                                                             >  1:713767‑M1/1‑142 (MQ=255)    ttacaGCCCGCTACCCGCCATAGCGTTCGCCAGACCAGCGCCTTCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgc                                                                                                                                          <  1:609451‑M1/141‑5 (MQ=255)           ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGATATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAAgcgaaagctcg                                                                                                                                   >  1:466183‑M1/1‑132 (MQ=255)                        gagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggg                                                                                                                      <  1:753444‑M1/143‑25 (MQ=255)                                gCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                              >  1:628390‑M1/1‑111 (MQ=255)                                        acccccTCCCGTATCGGGCGGCGGCGTCATTGCAGTGACCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  2:658207‑M1/143‑41 (MQ=255)                                           accTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      >  2:456732‑M1/1‑100 (MQ=255)                                           accTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:456732‑M1/140‑41 (MQ=255)                                             cTTCCGTATTGGGCGGCGGCGTCAGTTCAGGGAGCGGCTCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAAGCCGGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccggaactcgccggga                                                                                                 <  2:621938‑M1/143‑46 (MQ=255)                                              ttCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggac                                                                                                >  1:93202‑M1/1‑97 (MQ=255)                                              ttCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaccaggacgataccgcacctcaccgcgac                                                                                                >  2:649191‑M1/1‑97 (MQ=255)                                                tcctattGGGCGGTGGCGTTTTTTTAGTGAGCTGCGCCAGCGACCTACCGTCGATACCGCGAGTGATGTTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtaaagaacaggcctatacctaaactcaccgggacgg                                                                                              <  2:332218‑M1/140‑49 (MQ=255)                                                  gTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcg                                                                                            >  1:909641‑M1/1‑93 (MQ=255)                                                   tagcgggcggcGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                           <  1:407049‑M1/139‑52 (MQ=255)                                                   tATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                           >  1:101646‑M1/1‑92 (MQ=255)                                                    atggggcggCGGCGTCATTCCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCTAGgcgcaagctcagtcaggaacagggcgataccgaagctcaccgggacggcgat                                                                                          <  1:484574‑M1/140‑53 (MQ=255)                                                    aTTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  1:623220‑M1/143‑53 (MQ=255)                                                       gggCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcag                                                                                       <  1:561916‑M1/143‑56 (MQ=255)                                                            gcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcg                                                                                  <  1:812032‑M1/143‑61 (MQ=255)                                                              ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                <  1:473383‑M1/143‑63 (MQ=255)                                                                     tggcagtGAGCAGCGCCAGCGACATACCGGCGATGCCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCAAGgcgcaagctcagtcaggaacagggcgataccgaaactaaccgggacggcgatcagcagcgcgataaacg                                                                         <  2:253131‑M1/140‑70 (MQ=255)                                                                       tCAGTGAACAGCGCCAGCGACAGACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgagactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  1:818682‑M1/143‑72 (MQ=255)                                                                          gTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaaactcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatcgacgaagtc                                                                    >  2:756123‑M1/1‑69 (MQ=255)                                                                           ttacccgcgccagcgaCAGACCGGCGATACCGCGATCGAGTGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggtacagggcgatatcgaaactcaccgtgacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:103759‑M1/137‑76 (MQ=255)                                                                           tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:248148‑M1/143‑76 (MQ=255)                                                                             aaCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacccggaaggcgttcagcctcgcgataaacgaagtcacc                                                                 >  1:722260‑M1/1‑66 (MQ=255)                                                                               cAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaa                                                               >  2:527728‑M1/1‑64 (MQ=255)                                                                                aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaaggcgataattaaactcaccggcacggagcacagcagcgcaagcaaagaagtaaccaac                                                              >  2:107004‑M1/1‑63 (MQ=255)                                                                                    ccAGCGACATACCGTCGAGACCGCGAGTGAGGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacgacgcaccgccgcgcgataaacgaagtcaccaacgtac                                                          <  2:566187‑M1/143‑85 (MQ=255)                                                                                       gCGACATACCGCCGATACCGCGAGTGACGGTGGACATTAAAATCCAGTTCAGCCAGccgcaagctcagtcaggaacagggcgataccgacactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgt                                                       <  1:110940‑M1/143‑88 (MQ=255)                                                                                             tACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaacagtgcgataccgaaactcaccgggacggcgatcagcatcgcgataaacgaagtcaccgacgtaccgtagatag                                                 >  1:444306‑M1/1‑50 (MQ=255)                                                                                              aCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                <  2:101646‑M1/143‑95 (MQ=255)                                                                                               ccGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggctataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggc                                               >  1:567397‑M1/1‑48 (MQ=255)                                                                                                     aTACCGCGAGTGATGGCGGACATTAAAATCCAGAACAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagc                                         >  1:708962‑M1/1‑42 (MQ=255)                                                                                                         cgcgGGGGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccc                                     <  2:470108‑M1/143‑106 (MQ=255)                                                                                                            gAGTCATGAGGGACAGTACACGCCAGACCAGCCACgcgcaagctcagtcaggaacagcgcggtaccgaacgtcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgt                                  <  2:505363‑M1/143‑109 (MQ=255)                                                                                                             aGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                 <  1:102536‑M1/143‑110 (MQ=255)                                                                                                               tGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtacggtacatcggcaccagcgcccggtgga                               >  2:675023‑M1/1‑32 (MQ=255)                                                                                                                   tgtgtcccttacaatccaGCTCCGCCACgcgcaagctccgtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatc                           <  1:878354‑M1/135‑116 (MQ=255)                                                                                                                     tggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtcgatcggcaccagcgccccgtagctatcgt                         >  1:181914‑M1/1‑26 (MQ=255)                                                                                                                     tggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcaataccaaaacacaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgt                         >  1:264276‑M1/1‑26 (MQ=255)                                                                                                                        aCATTAGAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcg                      <  2:901907‑M1/143‑121 (MQ=255)                                                                                                                         ccttaaAATCCAGATCGCCCAGgcgcaggcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcgg                     <  2:131351‑M1/141‑122 (MQ=255)                                                                                                                          aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       >  2:335040‑M1/1‑21 (MQ=255)                                                                                                                          aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  1:335040‑M1/92‑72 (MQ=255)                                                                                                                          aTTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    <  1:536190‑M1/143‑123 (MQ=255)                                                                                                                                  ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatccca            >  1:365145‑M1/1‑13 (MQ=255)                                                                                                                                    aGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgcagtcaccaacgtcgcgtagatgggcaccagcgccccgtagaaattatttggtgcatcacgcg          >  1:465988‑M1/1‑11 (MQ=255)                                                                                                                                           ccAGgcgcacgcccggccaggaacagagcgataccgtaactcaccgggacggcgatcagcagcgcgataaacgacgtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggt   <  2:592993‑M1/143‑140 (MQ=255)                                                                                                                                            cAGgcgcaagctcagtcaggaacagggcgataccgaacctcacccgggcggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtc  <  1:253601‑M1/143‑141 (MQ=255)                                                                                                                                              |                                                                                                                                             AATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909169‑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: ‑