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

ACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909944‑3910224                                                                                                                                              |                                                                                                                                            attaacagcccgctacccgccagagcgttcgccagaccaccaccttccgcattgggcggcggcgtcatatcagtggacagcgccagcgacataccgccgatacggagagtgaggaggcacattaaaatcacgaccagccaggc                                                                                                                                            >  1:190370‑M2/142‑143 (MQ=255)          ccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgaacgccacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtca                                                                                                                                   <  1:252987‑M2/11‑1 (MQ=255)                      agagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAgg                                                                                                                       <  2:89671‑M2/23‑1 (MQ=255)                        agcgttcgccagcccaccaccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGc                                                                                                                     <  1:76091‑M2/25‑1 (MQ=255)                               gccagacaaccacctgccggattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcaataccgcgagtgatggtggacattaaaatcaagatcagccgaGCGCaagctcagtcaggaaagtggacctgccc                                                                                                              >  1:337106‑M2/112‑130 (MQ=255)                                 cagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGACCAGGGCGATACCGaa                                                                                                            >  1:340746‑M2/110‑143 (MQ=255)                                        ccaccttccgtattggtcggcggcgtcatttcagttaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaacatcccgatcagccagGCGCACGCTCAGTCAGGAACAGGGCGATACCGAAACTCAcc                                                                                                     <  1:346294‑M2/41‑1 (MQ=255)                                         caccttccgtattgggcggcggcgtcatttcagtaaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCCAAACTCACCg                                                                                                    >  1:133078‑M2/102‑143 (MQ=255)                                                   gttgggcggcggcgtcattgcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                          <  2:211083‑M2/52‑1 (MQ=255)                                                   attgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGCAACTCACCGGGACGGCGAt                                                                                          >  1:181401‑M2/92‑143 (MQ=255)                                                   attgggcggcggcgtcatttcagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGAt                                                                                          <  2:64167‑M2/52‑1 (MQ=255)                                                     tgggcggtggcgtcattttagtgaacagcgccagcgacgtaccgtcgataccgcgagtgatggtggacattaaaatccaggtcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATca                                                                                        <  1:350383‑M2/54‑1 (MQ=255)                                                         cggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGACAAGGACCATACCGAAACTCACCAGGACGGCGATcagcag                                                                                    >  1:422764‑M2/86‑143 (MQ=255)                                                                    tttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGACACTCACCGGGACGGCCATCAGCAGCGCGCTAGACg                                                                         >  1:462892‑M2/75‑143 (MQ=255)                                                                                  gcaagcgtcttaccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTa                                                           <  2:293925‑M2/83‑1 (MQ=255)                                                                                   ccagcgtcataccgtcgatagcgcgagggagggtggacattaaaagccatatgggccagGCGCAAGCTCCGTCAGGAACGGGGTGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTAc                                                          <  1:523937‑M2/84‑1 (MQ=255)                                                                                           ataccgtcgataccgcgagtgatggtggacattcaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGAGACGGCGACCAGCAGCGCGATAAAcaaagtcaccaaccagccaccacgc                                                  >  1:49845‑M2/52‑131 (MQ=255)                                                                                               cgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGACCAGGGCGATACCGAAACTCACCGAGAAGTAGTTCAGCAGCGCGATAAACGAAGTCACCATCGTACCGTAGATAGGCa                                              >  2:219244‑M2/48‑143 (MQ=255)                                                                                                 tcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATAACGAAATTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCGGCAACGTACCGTAGATCGGCAcc                                            >  2:263248‑M2/46‑143 (MQ=255)                                                                                                                 tggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGAAGacat                            >  1:428072‑M2/30‑143 (MQ=255)                                                                                                                         attaaaatcaagatcagcgtgTCGCAAGCTTAGTCAGCAACAGGGCGATACCGAAACTCACCGGGACGGCGCTCAGCAGCGCGATAAACGAAGCCACCACCGTACCGTAGATCGGCACCAGCGCCCCGTAGAactcgtgccgc                    >  1:347357‑M2/22‑138 (MQ=255)                                                                                                                           taacgtccggatcagccacGCGCAAGCTCCGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGc                  <  2:188913‑M2/124‑1 (MQ=255)                                                                                                                             aaatccagatcagccagTCGCAAGCTCAGTCAGGAACAGGGCGATCACCAATCTCACCGAAACGATGGACAGCAGTACTATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCAt                <  2:266795‑M2/126‑1 (MQ=255)                                                                                                                                  cagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTCGATATCGATCGGTGCATCCCAc           >  1:411717‑M2/13‑143 (MQ=255)                                                                                                                                         gccagGCGCAAGCTCAGTCAGGAGCAGGGCGATACCGATCCGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTgg    <  2:166096‑M2/138‑1 (MQ=255)                                                                                                                                         gccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTgg    <  2:232001‑M2/138‑1 (MQ=255)                                                                                                                                           cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGAGCGAGGAACGGGGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTc  >  2:471689‑M2/4‑143 (MQ=255)                                                                                                                                              |                                                                                                                                            ACATGATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTC  >  NC_000913/3909944‑3910224

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