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	57 [0]	[0] 57	[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)	54 (0.860)	48/268	0.3	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

GATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909949‑3910226                                                                                                                                          |                                                                                                                                             ccaggccgcgacccgccggcgcgttcgccagaccgccgccgtccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCaa                                                                                                                                         <  2:535385‑M2/6‑1 (MQ=255) acagcccgctacccgccagagccttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacattccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCaa                                                                                                                                         <  1:619210‑M2/6‑1 (MQ=255)   agcccgctacccgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacctaccgtcgataccgcgagtgatggtggacattacaatccagatcagccagGCGCAAGc                                                                                                                                       >  1:72506‑M2/136‑143 (MQ=255)     cccgctacccgccagagcgttcgccagaccacccccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtc                                                                                                                                     <  2:965955‑M2/10‑1 (MQ=255)     cccgctacccgccagagcgttcgccagaccaccaccttccgtactgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtc                                                                                                                                     >  2:1042648‑M2/134‑143 (MQ=255)          tacccgccagagcgttcgccataccaccacctgccgtattgggcggctgcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCtcagtca                                                                                                                                <  2:131317‑M2/15‑1 (MQ=255)              cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGaa                                                                                                                            <  1:1045702‑M2/19‑1 (MQ=255)              cgccagagcgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGaa                                                                                                                            >  2:913357‑M2/125‑143 (MQ=255)               gcaagggcgttcgccagaccaccaccttccggattgggcggcggcgtcatttcagtgaacagcgccagcgagatgccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAAc                                                                                                                           <  2:196400‑M2/20‑1 (MQ=255)                      cgttcgccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacatacggtcgataccgcgagtgatggtggacattaagatccagatcagccagGCGCGAGCTCAGTCAGGAGCAGGGCGa                                                                                                                    >  1:474154‑M2/117‑143 (MQ=255)                           gccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCg                                                                                                               <  2:898663‑M2/32‑1 (MQ=255)                            ccagaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagctacatcccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGa                                                                                                              <  2:1027373‑M2/33‑1 (MQ=255)                              agaccgccgccgcccggagtgggcggcgccgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGaaa                                                                                                            <  2:401286‑M2/35‑1 (MQ=255)                               gaccaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAAc                                                                                                           >  2:412771‑M2/108‑143 (MQ=255)                                accaccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGCTACcgaaaat                                                                                                          >  1:59425‑M2/107‑141 (MQ=255)                                  caccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccatatcagccagGAGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCa                                                                                                        >  1:1052711‑M2/105‑143 (MQ=255)                                   accaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacatacagtcgattccgcgagtgagggtggacattaaaatccagatcagccagACGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAc                                                                                                       >  1:1063936‑M2/104‑143 (MQ=255)                                    ccaccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCAcc                                                                                                      >  1:953137‑M2/103‑143 (MQ=255)                                     caccttccgtattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCg                                                                                                     >  1:365760‑M2/102‑143 (MQ=255)                                          tccgtattgggcggcggcgtcatttcagtgaacagcgacagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCATGAACCGCGCGATACCGACACTCACCGGGGCg                                                                                                >  2:254328‑M2/97‑143 (MQ=255)                                              tattgggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGa                                                                                            <  1:788295‑M2/51‑1 (MQ=255)                                                   ggcggcggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgcgtgatggtggacattaaaatccagatcagccagGCCCAAGCCCAGGCAGGAACAGGGCGATACCGAACCTCACCCGGACGGCGagcacg                                                                                       >  2:458317‑M2/88‑141 (MQ=255)                                                        cggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACAGCGGCCAGCAgcgc                                                                                  >  1:1113453‑M2/83‑143 (MQ=255)                                                         ggcgtcatttcagtgatcagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                 <  1:1148792‑M2/62‑1 (MQ=255)                                                         ggcgtcatttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGcgcg                                                                                 >  1:1014228‑M2/82‑143 (MQ=255)                                                               accacagtgaacagcgccagcgacataccggcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAAc                                                                           <  1:728025‑M2/68‑1 (MQ=255)                                                                tttcagtgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTACGTCACGATCAGGGCGATACCGAAACCCACCGGGCCGGCGATCAGCAGCGAGAGCAACg                                                                          >  2:894440‑M2/75‑143 (MQ=255)                                                                  tgagtgaccagcgccagctgcataccgtcgataccgcgagtgatggtggagattaaaatccagatcagccagGCGCAAGGTCAGTCAGGAACAGGGCGATACCGAAGCTCACCGGGACGGCGATCAGCAGCGCGATAAACGaa                                                                        <  2:1096868‑M2/71‑1 (MQ=255)                                                                      cgaacagcgccagcgacataccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGACCAGGGCCATCCCGGATCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCa                                                                    <  1:612623‑M2/75‑1 (MQ=255)                                                                           agcgccagcgacgtaccgtcgataccgcgagtgatggtggacattaaaagccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAAc                                                               <  1:679095‑M2/80‑1 (MQ=255)                                                                                         accgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGCAGTCACGAACGTAccgtagcgcgg                                                 >  2:84436‑M2/50‑138 (MQ=255)                                                                                         accgccgataccgcgagtgatggtggacattaaaatccagctcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                 <  1:789890‑M2/94‑1 (MQ=255)                                                                                         accgccgataccgcgagggatgggggacattaaaatccagagcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCgg                                                 <  1:992089‑M2/94‑1 (MQ=255)                                                                                          ccgtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGAGAACAGCAGCGCGACAAACGAAGTCACCAACGTACCGTAGACCGGc                                                >  2:468181‑M2/49‑143 (MQ=255)                                                                                            gtcgataccgcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCAc                                              <  1:1042648‑M2/97‑1 (MQ=255)                                                                                              cgataccgcgtgttatgggggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCa                                            <  2:1052711‑M2/99‑1 (MQ=255)                                                                                                   gcgcgagtgacggtggcccttaaacaccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGcc                                       <  1:488571‑M2/104‑1 (MQ=255)                                                                                                   ccgcgagtgatggtggacattaaaatccaaatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTAACCGGGACGACGATCAGCAGCGCGAAAAACGCAGTCACCAACGTACGGTAGACCGGCACCAACGcc                                       >  2:212873‑M2/40‑143 (MQ=255)                                                                                                     gcgagtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCgcccg                                     >  1:158711‑M2/38‑142 (MQ=255)                                                                                                        agtgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATTAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                  <  1:510507‑M2/109‑1 (MQ=255)                                                                                                        agtgatggtggacattaaaacccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTa                                  <  1:412771‑M2/109‑1 (MQ=255)                                                                                                          tgatggtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGa                                >  1:956421‑M2/33‑143 (MQ=255)                                                                                                               gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTGCCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           >  1:148184‑M2/28‑143 (MQ=255)                                                                                                               gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCg                           >  2:939356‑M2/28‑143 (MQ=255)                                                                                                               gtggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGGCACCAACGGACCGTAGATCGGCACCAGCGCCCCGGAGATGTCg                           >  2:784617‑M2/28‑143 (MQ=255)                                                                                                                 ggacattaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCGGCAGCGCGATAAACGAAGTCACCAACGTACCGCAGATAGGCACCAGAGCCCCGTAGACAGCgta                         >  2:92620‑M2/26‑142 (MQ=255)                                                                                                                  gacattaaaatccagatcaaccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTc                        >  1:250987‑M2/25‑143 (MQ=255)                                                                                                                     attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGCTCAGCAGCCCGATAAACGAAGTCACCAACGTACCGTAGCGAGGCACCCGCGCCCCGTAGATATCGTTCGGt                     >  1:599835‑M2/22‑143 (MQ=255)                                                                                                                     attaaaatccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGt                     >  1:398369‑M2/22‑143 (MQ=255)                                                                                                                           gtccagatcagcaagGCGCAAGCTCAGGCAGGAACAGGGCGATACCGAAACGCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATcc               <  1:982752‑M2/128‑1 (MQ=255)                                                                                                                             ccagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCCCCGGGCCAGCAAATAGCAGCGCGATAACCGAAGTCACCAACGTACCGTAGATCGGCGCCAGCGCCCCGTAAATAACGCtcggagccagccc             >  2:780566‑M2/14‑137 (MQ=255)                                                                                                                              cagatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCAc            <  1:772486‑M2/131‑1 (MQ=255)                                                                                                                                gatcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCActc          >  1:953673‑M2/11‑143 (MQ=255)                                                                                                                                  tcagccagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCtt        <  1:458454‑M2/135‑1 (MQ=255)                                                                                                                                    agccagGCGCAAGCTCAGTCCGGGACAGGGCGATACCGAAGCTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTg      <  2:1183175‑M2/137‑1 (MQ=255)                                                                                                                                       cagGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCCATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGAGATCGTTCGGTGCATCGCACTCTTTGGTc   >  1:251356‑M2/4‑143 (MQ=255)                                                                                                                                        agGAGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATTAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGAGCTCCGTAGATATCGTTCGGTGCATCCGACTCTTTggtac  >  1:619033‑M2/3‑141 (MQ=255)                                                                                                                                          |                                                                                                                                             GATATTGCCGACCGGCTCCTGAAAGTAAACGGCGAACAGCGGCGCAAAGATAAACAGGCCCCACATGCCGTAAACGATACTTGGAATGGCTGCCAGCAGCTCAATGGCGATACCCAGCGGGCGTTTCAGCCAGCCAGGCGCAAGCTCAGTCAGGAACAGGGCGATACCGAAACTCACCGGGACGGCGATCAGCAGCGCGATAAACGAAGTCACCAACGTACCGTAGATCGGCACCAGCGCCCCGTAGATATCGTTCGGTGCATCCCACTCTTTGGTCC  >  NC_000913/3909949‑3910226

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