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	25 [0]	[0] 26	[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)	25 (0.650)	22/270	0.9	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

ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTA  >  NC_000913/3909171‑3909441                                                                                                                                             |                                                                                                                                   cttaACCGCCCGCGACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACCGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgc                                                                                                                                  <  1:349932‑M1/142‑3 (MQ=255)    accagCCCGCCACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgca                                                                                                                               <  1:572870‑M1/141‑6 (MQ=255)     aCAGCCCGCTACCCGCCAGAGCGGTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                              <  2:643416‑M1/143‑7 (MQ=255)         cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTACTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGCTGGGGGACATCGAAATCCAGACCAGCCAGgggcgagcac                                                                                                                          >  2:575422‑M1/1‑133 (MQ=255)           cGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcag                                                                                                                        <  1:687720‑M1/143‑13 (MQ=255)                      agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCGGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagg                                                                                                             >  1:9040‑M1/1‑120 (MQ=255)                        ggggttcGCCAGACACCCGCCTGCGCTATTGGGCGGCGGCGTGATGTCAGTGAACAGGGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggc                                                                                                           <  1:523408‑M1/140‑26 (MQ=255)                             tCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgatac                                                                                                      <  1:559093‑M1/143‑31 (MQ=255)                                ccAGACCCCCACCTTCCGTATGGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccga                                                                                                   <  1:466834‑M1/143‑34 (MQ=255)                                   gACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGGTACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgacaccgaaac                                                                                                >  2:357226‑M1/1‑107 (MQ=255)                                    accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAAGCCAGCTCAGCCAGgcgcaggctcagtcaggaccagggcgataccgcgacg                                                                                               >  1:195502‑M1/1‑106 (MQ=255)                                                          ggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgaccagcagc                                                                         >  1:514962‑M1/1‑84 (MQ=255)                                                                 tCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcggtcagccgcgcgataa                                                                  >  2:73371‑M1/1‑77 (MQ=255)                                                                               agcgccagcgccATACCGTCGATACCGCGAGTGATGGTGGCCATTAAAATCCAGATCAGCCAGgcgccagcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaggccacccac                                                    >  1:246871‑M1/1‑63 (MQ=255)                                                                               aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                    <  1:27196‑M1/143‑81 (MQ=255)                                                                                  gcCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacacggcgataccgaaactcaccgggacggcgctcagcagcgcgataaacgaagccaccaccgcc                                                 >  2:183705‑M1/1‑60 (MQ=255)                                                                                                        cgcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgagcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccc                           >  2:651714‑M1/1‑38 (MQ=255)                                                                                                              tGATGGTGGACATTAACATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                     <  1:193784‑M1/143‑112 (MQ=255)                                                                                                              cgaTGGTGCACAGTAACAGCCAGATCAGCCAGgcgcaagcgcagtcaggagcagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                     <  2:862‑M1/142‑112 (MQ=255)                                                                                                                   gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggcgcagggcgataccggccctcaccgggacggcgatcagcagcgcgataagcgaagtcaccaacggaccgtagatcggcaccagcgccccgtagatatcg                <  1:272435‑M1/143‑117 (MQ=255)                                                                                                                   gtggACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                <  1:656913‑M1/143‑117 (MQ=255)                                                                                                                      gACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttc             >  1:519480‑M1/1‑24 (MQ=255)                                                                                                                        cATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcgg           >  1:362207‑M1/1‑22 (MQ=255)                                                                                                                          ttAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtg         <  1:694599‑M1/143‑124 (MQ=255)                                                                                                                                 ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatccca  <  1:663545‑M1/143‑131 (MQ=255)                                                                                                                                             |                                                                                                                                   ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTA  >  NC_000913/3909171‑3909441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑