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	48 [0]	[0] 48	[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)	46 (0.950)	40/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

GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909944‑3910225                                                                                                                                               |                                                                                                                                            aattaacagcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatgacgggcatcaaactcccgatcagccagg                                                                                                                                             >  1:713767‑M2/143‑143 (MQ=255)    ttacagcccgctacccgccatagcgttcgccagaccagcgccttccggattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaggcgc                                                                                                                                          <  1:609451‑M2/4‑1 (MQ=255)           ccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgatataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccaaGCGAAAGctcg                                                                                                                                   >  1:466183‑M2/133‑142 (MQ=255)                        gagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAggg                                                                                                                      <  1:753444‑M2/24‑1 (MQ=255)                                gccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                              >  1:628390‑M2/112‑143 (MQ=255)                                        accccctcccgtatcgggcggcggcgtcattgcagtgaccagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      <  2:658207‑M2/40‑1 (MQ=255)                                           accttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      >  2:456732‑M2/101‑140 (MQ=255)                                           accttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                      <  1:456732‑M2/40‑1 (MQ=255)                                             cttccgtattgggcggcggcgtcagttcagggagcggctccagcgacataccggcgataccgcgagtgatggtggacattaaaagccggatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGGAACTCGCCGGGa                                                                                                 <  2:621938‑M2/45‑1 (MQ=255)                                              ttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGAc                                                                                                >  1:93202‑M2/98‑143 (MQ=255)                                              ttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCGCAGTCAGGACCAGGACGATACCGCACCTCACCGCGAc                                                                                                >  2:649191‑M2/98‑143 (MQ=255)                                                tcctattgggcggtggcgtttttttagtgagctgcgccagcgacctaccgtcgataccgcgagtgatgttggacattaaaatccagatcagccagGCGCAAGCTCAGTAAAGAACAGGCCTATACCTAAACTCACCGGGACgg                                                                                              <  2:332218‑M2/48‑1 (MQ=255)                                                  gtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCg                                                                                            >  1:909641‑M2/94‑143 (MQ=255)                                                   tattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                           >  1:101646‑M2/93‑143 (MQ=255)                                                   tagcgggcggcggcgtcatttcagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                           <  1:407049‑M2/51‑1 (MQ=255)                                                    attgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                          <  1:623220‑M2/52‑1 (MQ=255)                                                    atggggcggcggcgtcattccagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagctagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGGGACGGCGAt                                                                                          <  1:484574‑M2/52‑1 (MQ=255)                                                       gggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATcag                                                                                       <  1:561916‑M2/55‑1 (MQ=255)                                                            gcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAgcg                                                                                  <  1:812032‑M2/60‑1 (MQ=255)                                                              ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                <  1:473383‑M2/62‑1 (MQ=255)                                                                     tggcagtgagcagcgccagcgacataccggcgatgccgcgagtgatggtggacattaaaatccagatcagcaagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGGCGATCAGCAGCGCGATAAACg                                                                         <  2:253131‑M2/69‑1 (MQ=255)                                                                       tcagtgaacagcgccagcgacagaccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAGACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       <  1:818682‑M2/71‑1 (MQ=255)                                                                          gtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAACTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATCGACGAAGTc                                                                    >  2:756123‑M2/70‑143 (MQ=255)                                                                           ttacccgcgccagcgacagaccggcgataccgcgatcgagtgtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGTACAGGGCGATATCGAAACTCACCGTGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                   <  1:103759‑M2/75‑1 (MQ=255)                                                                           tgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                   <  1:248148‑M2/75‑1 (MQ=255)                                                                             aacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCCGGAAGGCGTTCAGCCTCGCGATAAACGAAGTCAcc                                                                 >  1:722260‑M2/67‑143 (MQ=255)                                                                               cagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCaa                                                               >  2:527728‑M2/65‑143 (MQ=255)                                                                                agcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAAGGCGATAATTAAACTCACCGGCACGGAGCACAGCAGCGCAAGCAAAGAAGTAACCAAc                                                              >  2:107004‑M2/64‑143 (MQ=255)                                                                                    ccagcgacataccgtcgagaccgcgagtgagggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGACGCACCGCCGCGCGATAAACGAAGTCACCAACGTAc                                                          <  2:566187‑M2/84‑1 (MQ=255)                                                                                       gcgacataccgccgataccgcgagtgacggtggacattaaaatccagttcagccagCCGCAAGCTCAGTCAGGAACAGGGCGATACCGACACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGt                                                       <  1:110940‑M2/87‑1 (MQ=255)                                                                                             taccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCGCAGTCAGGAACAGTGCGATACCGAAACTCACCGGGACGGCGATCAGCATCGCGATAAACGAAGTCACCGACGTACCGTAgatag                                                 >  1:444306‑M2/51‑141 (MQ=255)                                                                                              accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                <  2:101646‑M2/94‑1 (MQ=255)                                                                                               ccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCTATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGc                                               >  1:567397‑M2/49‑143 (MQ=255)                                                                                                     ataccgcgagtgatggcggacattaaaatccagaacagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAgc                                         >  1:708962‑M2/43‑143 (MQ=255)                                                                                                         cgcgggggatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGccc                                     <  2:470108‑M2/105‑1 (MQ=255)                                                                                                             agtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                 <  1:102536‑M2/109‑1 (MQ=255)                                                                                                               tgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACGGTACATCGGCACCAGCGCCCGGTGGa                               >  2:675023‑M2/33‑143 (MQ=255)                                                                                                                     tggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTCGATCGGCACCAGCGCCCCGTAGCTATCGt                         >  1:181914‑M2/27‑143 (MQ=255)                                                                                                                     tggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCAATACCAAAACACACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGt                         >  1:264276‑M2/27‑143 (MQ=255)                                                                                                                        acattagaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCg                      <  2:901907‑M2/120‑1 (MQ=255)                                                                                                                         ccttaaaatccagatcgcccagGCGCAGGCGCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCgg                     <  2:131351‑M2/121‑1 (MQ=255)                                                                                                                          attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       >  2:335040‑M2/22‑92 (MQ=255)                                                                                                                          attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                       <  1:335040‑M2/71‑1 (MQ=255)                                                                                                                          attaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                    <  1:536190‑M2/122‑1 (MQ=255)                                                                                                                                  ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCa            >  1:365145‑M2/14‑143 (MQ=255)                                                                                                                                    agatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGCAGTCACCAACGTCGCGTAGATGGGCACCAGCGCCCCGTAGAAATTATTTGGtgcatcacgcg          >  1:465988‑M2/12‑138 (MQ=255)                                                                                                                                           ccagGCGCACGCCCGGCCAGGAACAGAGCGATACCGTAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGACGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGt   <  2:592993‑M2/139‑1 (MQ=255)                                                                                                                                            cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAACCTCACCCGGGCGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTc  <  1:253601‑M2/140‑1 (MQ=255)                                                                                                                                               |                                                                                                                                            GACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909944‑3910225

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