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	33 [0]	[0] 32	[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)	32 (0.680)	29/270	0.7	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

AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAG  >  NC_000913/3909169‑3909435                                                                                                                                               |                                                                                                                             aaaTTGCCAGCCCGCTACCCGCCAGAGCGTTCGCCGTACCGCCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAG                                                                                                                              <  2:523432‑M1/143‑1 (MQ=255)            tcgCTACCCGCCGGGGCGTTCGCCAGACCTCCACCCTCCGTATGGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctca                                                                                                                   <  2:18522‑M1/142‑12 (MQ=255)                tACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcg                                                                                                               >  1:709390‑M1/1‑128 (MQ=255)                     gCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaac                                                                                                          >  1:136225‑M1/1‑123 (MQ=255)                      ccAGGGCGTTCGCCAGCCCACCACCTTCCGTATCGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaaca                                                                                                         <  2:571108‑M1/143‑22 (MQ=255)                          agCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTAGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggc                                                                                                     >  1:177989‑M1/1‑118 (MQ=255)                            cGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcga                                                                                                   <  1:749181‑M1/143‑28 (MQ=255)                                   cAGCCCACCATCCTCCGTATTGGGCGGCGGCGTCACGGCAGTGAACAGCGCCAGCGACCTACCGTCGATACCGCGAGTGATCGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggacgataccgaa                                                                                            <  2:448651‑M1/143‑35 (MQ=255)                                   cAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggg                                                                                                      >  2:62032‑M1/1‑109 (MQ=255)                                   cAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggg                                                                                                      <  1:62032‑M1/133‑25 (MQ=255)                                     cagcaccgccttcCGTATGGGGCGGCGCCGCCACTTCAGTGAACAGCGCCAGCGACATCCCGTGGATACCGCGAGTGATGGTGGACATTAAAATCCATATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgacac                                                                                          <  1:366892‑M1/140‑37 (MQ=255)                                              cTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccggga                                                                                 >  1:536161‑M1/1‑98 (MQ=255)                                                    tATTGCGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATTGTGTACATTATAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                           <  1:447853‑M1/143‑52 (MQ=255)                                                    gagtgggCGGCGACGGAAGTTCAGTGAACGGCGCCAGCGACATACCGTCGACACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                           <  1:878495‑M1/140‑52 (MQ=255)                                                           cggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcag                                                                    >  1:297526‑M1/1‑85 (MQ=255)                                                             gcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcg                                                                  <  1:721482‑M1/143‑61 (MQ=255)                                                                        tCAGTGGACTGCGCCAGAGAGATAGGGTCGATACCGAGAGTGATTGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                       <  2:318403‑M1/143‑72 (MQ=255)                                                                           gTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatgagcagcgcgacaaacgaagtc                                                    >  2:786944‑M1/1‑69 (MQ=255)                                                                              aaCAGCGCCAGCGACATACGGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcacc                                                 <  1:495887‑M1/143‑78 (MQ=255)                                                                                gcacgccagCGCCTTAGCGTCGCTACCGGGGGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaaggtcagtcaggaacagggctataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaa                                               <  2:762033‑M1/140‑80 (MQ=255)                                                                                  gcgcCAGCGACAGACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacg                                             <  2:34100‑M1/143‑82 (MQ=255)                                                                                      cAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtacc                                         >  1:831780‑M1/1‑58 (MQ=255)                                                                                      aagCGACATACCGTCGATACAGCGAGGGATGGTGGACAGTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtacc                                         <  1:810648‑M1/142‑86 (MQ=255)                                                                                         cGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacggggacggcgatcagcggcgagacaaaagaagacaccaacgtacggta                                      >  2:366751‑M1/1‑55 (MQ=255)                                                                                             aTACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaaccggtcggtaccgaaactcaccgggacgacgatcagcagagggctcaacgacgtcaccaacgtaccgtagagc                                  >  1:563231‑M1/1‑51 (MQ=255)                                                                                              gacCGTCGAGACCGCGAGTGATGGTGGACAAGAAAATCGAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcg                                 <  1:294649‑M1/142‑94 (MQ=255)                                                                                                   tCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagagcggcacc                            >  1:542569‑M1/1‑45 (MQ=255)                                                                                                      aTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagc                         >  1:534311‑M1/1‑42 (MQ=255)                                                                                                      aTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcaacagcgcgacaaacgaagtcaccaacgtacgctagctcggcaccagc                         >  2:416693‑M1/1‑42 (MQ=255)                                                                                                                     gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg          >  2:719445‑M1/1‑27 (MQ=255)                                                                                                                     gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgacaaacgcagtcaccaacgtaccgtagctcggcaccagcgccccggagatatcg          >  2:719451‑M1/1‑27 (MQ=255)                                                                                                                           aTTAAACTCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaagctcaccgagcaggctgtgttcagcgcgagaaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt    <  1:903679‑M1/143‑123 (MQ=255)                                                                                                                             tAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgacaccagcgccccgtagctatcgtgcggcgc  >  1:587969‑M1/1‑19 (MQ=255)                                                                                                                                               |                                                                                                                             AAATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAG  >  NC_000913/3909169‑3909435

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