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	27 [0]	[0] 28	[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)	27 (0.670)	24/270	0.7	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

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

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