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] 27	[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)	24 (0.840)	21/270	0.3	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

ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909170‑3909450                                                                                                                                             |                                                                                                                                             aTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGCATTGGGCGGCGGCGTCATATCAGTGGACAGCGCCAGCGACATACCGCCGATACGGAGAGTGAGGAGGCACATTAAAATCACGACCAGCCAGgc                                                                                                                                            >  1:190370‑M1/1‑141 (MQ=255)          ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGAACGCCACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctca                                                                                                                                   <  1:252987‑M1/143‑12 (MQ=255)                      agagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagg                                                                                                                       <  2:89671‑M1/143‑24 (MQ=255)                        agCGTTCGCCAGCCCACCACCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggc                                                                                                                     <  1:76091‑M1/143‑26 (MQ=255)                               gCCAGACAACCACCTGCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCAATACCGCGAGTGATGGTGGACATTAAAATCAAGATCAGCCGAgcgcaagctcagtcaggaaagtggacctgccc                                                                                                              >  1:337106‑M1/1‑111 (MQ=255)                                 cAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaccagggcgataccgaa                                                                                                            >  1:340746‑M1/1‑109 (MQ=255)                                        ccaccTTCCGTATTGGTCGGCGGCGTCATTTCAGTTAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAACATCCCGATCAGCCAGgcgcacgctcagtcaggaacagggcgataccgaaactcacc                                                                                                     <  1:346294‑M1/143‑42 (MQ=255)                                         caccTTCCGTATTGGGCGGCGGCGTCATTTCAGTAAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgatacccaaactcaccg                                                                                                    >  1:133078‑M1/1‑101 (MQ=255)                                                   gttgGGCGGCGGCGTCATTGCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  2:211083‑M1/142‑53 (MQ=255)                                                   aTTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgcaactcaccgggacggcgat                                                                                          >  1:181401‑M1/1‑91 (MQ=255)                                                   aTTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  2:64167‑M1/143‑53 (MQ=255)                                                     tGGGCGGTGGCGTCATTTTAGTGAACAGCGCCAGCGACGTACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGGTCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatca                                                                                        <  1:350383‑M1/143‑55 (MQ=255)                                                         cggcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggacaaggaccataccgaaactcaccaggacggcgatcagcag                                                                                    >  1:422764‑M1/1‑85 (MQ=255)                                                                    tttCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgacactcaccgggacggccatcagcagcgcgctagacg                                                                         >  1:462892‑M1/1‑74 (MQ=255)                                                                                  gcaagcgTCTTACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgta                                                           <  2:293925‑M1/140‑84 (MQ=255)                                                                                   ccAGCGTCATACCGTCGATAGCGCGAGGGAGGGTGGACATTAAAAGCCATATGGGCCAGgcgcaagctccgtcaggaacggggtgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtac                                                          <  1:523937‑M1/143‑85 (MQ=255)                                                                                           aTACCGTCGATACCGCGAGTGATGGTGGACATTCAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgagacggcgaccagcagcgcgataaacaaagtcaccaaccagccaccacgc                                                  >  1:49845‑M1/1‑51 (MQ=255)                                                                                               cGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaccagggcgataccgaaactcaccgagaagtagttcagcagcgcgataaacgaagtcaccatcgtaccgtagataggca                                              >  2:219244‑M1/1‑47 (MQ=255)                                                                                                 tCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataacgaaattcaccgggacggcgatcagcagcgcgataaacgaagtcggcaacgtaccgtagatcggcacc                                            >  2:263248‑M1/1‑45 (MQ=255)                                                                                                                 tggtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgaagacat                            >  1:428072‑M1/1‑29 (MQ=255)                                                                                                                         aTTAAAATCAAGATCAGCGTGtcgcaagcttagtcagcaacagggcgataccgaaactcaccgggacggcgctcagcagcgcgataaacgaagccaccaccgtaccgtagatcggcaccagcgccccgtagaactcgtgccgc                    >  1:347357‑M1/1‑21 (MQ=255)                                                                                                                           taacgtccggaTCAGCCACgcgcaagctccgtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgc                  <  2:188913‑M1/138‑125 (MQ=255)                                                                                                                             aaaTCCAGATCAGCCAGtcgcaagctcagtcaggaacagggcgatcaccaatctcaccgaaacgatggacagcagtactataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcat                <  2:266795‑M1/143‑127 (MQ=255)                                                                                                                                  cAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtcgatatcgatcggtgcatcccac           >  1:411717‑M1/1‑12 (MQ=255)                                                                                                                                         gCCAGgcgcaagctcagtcaggagcagggcgataccgatccgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttgg    <  2:166096‑M1/143‑139 (MQ=255)                                                                                                                                         gCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttgg    <  2:232001‑M1/143‑139 (MQ=255)                                                                                                                                           cAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagagcgaggaacggggtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtc  >  2:471689‑M1/1‑3 (MQ=255)                                                                                                                                             |                                                                                                                                             ATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909170‑3909450

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