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

ACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909175‑3909452                                                                                                                                         |                                                                                                                                              ccaGGCCGCGACCCGCCGGCGCGTTCGCCAGACCGCCGCCGTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                                         <  2:535385‑M1/142‑7 (MQ=255) aCAGCCCGCTACCCGCCAGAGCCTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATTCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaa                                                                                                                                         <  1:619210‑M1/143‑7 (MQ=255)   aGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACCTACCGTCGATACCGCGAGTGATGGTGGACATTACAATCCAGATCAGCCAGgcgcaagc                                                                                                                                       >  1:72506‑M1/1‑135 (MQ=255)     cccGCTACCCGCCAGAGCGTTCGCCAGACCACCCCCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctc                                                                                                                                     <  2:965955‑M1/143‑11 (MQ=255)     cccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTACTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctc                                                                                                                                     >  2:1042648‑M1/1‑133 (MQ=255)          tACCCGCCAGAGCGTTCGCCATACCACCACCTGCCGTATTGGGCGGCTGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtca                                                                                                                                <  2:131317‑M1/143‑16 (MQ=255)              cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaa                                                                                                                            <  1:1045702‑M1/143‑20 (MQ=255)              cGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaa                                                                                                                            >  2:913357‑M1/1‑124 (MQ=255)               gcaagggcgttcgCCAGACCACCACCTTCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGAGATGCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaac                                                                                                                           <  2:196400‑M1/137‑21 (MQ=255)                      cGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACGGTCGATACCGCGAGTGATGGTGGACATTAAGATCCAGATCAGCCAGgcgcgagctcagtcaggagcagggcga                                                                                                                    >  1:474154‑M1/1‑116 (MQ=255)                           gCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                               <  2:898663‑M1/143‑33 (MQ=255)                            ccAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCTACATCCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccga                                                                                                              <  2:1027373‑M1/143‑34 (MQ=255)                              aGACCGCCGCCGCCCGGAGTGGGCGGCGCCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaa                                                                                                            <  2:401286‑M1/143‑36 (MQ=255)                               gACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaac                                                                                                           >  2:412771‑M1/1‑107 (MQ=255)                                accaccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgctaccgaaaat                                                                                                          >  1:59425‑M1/1‑106 (MQ=255)                                  caccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCATATCAGCCAGgagcaagctcagtcaggaacagggcgataccgaaactca                                                                                                        >  1:1052711‑M1/1‑104 (MQ=255)                                   accaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACAGTCGATTCCGCGAGTGAGGGTGGACATTAAAATCCAGATCAGCCAGacgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                       >  1:1063936‑M1/1‑103 (MQ=255)                                    ccaccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacc                                                                                                      >  1:953137‑M1/1‑102 (MQ=255)                                     caccTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccg                                                                                                     >  1:365760‑M1/1‑101 (MQ=255)                                          tCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGACAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcatgaaccgcgcgataccgacactcaccggggcg                                                                                                >  2:254328‑M1/1‑96 (MQ=255)                                              tATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                            <  1:788295‑M1/143‑52 (MQ=255)                                                   ggcggcggcGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGCGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcccaagcccaggcaggaacagggcgataccgaacctcacccggacggcgagcacg                                                                                       >  2:458317‑M1/1‑87 (MQ=255)                                                        cggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacagcggccagcagcgc                                                                                  >  1:1113453‑M1/1‑82 (MQ=255)                                                         ggcgTCATTTCAGTGATCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                 <  1:1148792‑M1/143‑63 (MQ=255)                                                         ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                 >  1:1014228‑M1/1‑81 (MQ=255)                                                               accacagtgAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaac                                                                           <  1:728025‑M1/139‑69 (MQ=255)                                                                tttCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctacgtcacgatcagggcgataccgaaacccaccgggccggcgatcagcagcgagagcaacg                                                                          >  2:894440‑M1/1‑74 (MQ=255)                                                                  tgagtGACCAGCGCCAGCTGCATACCGTCGATACCGCGAGTGATGGTGGAGATTAAAATCCAGATCAGCCAGgcgcaaggtcagtcaggaacagggcgataccgaagctcaccgggacggcgatcagcagcgcgataaacgaa                                                                        <  2:1096868‑M1/141‑72 (MQ=255)                                                                      cgaACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaccagggccatcccggatctcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                    <  1:612623‑M1/142‑76 (MQ=255)                                                                           aGCGCCAGCGACGTACCGTCGATACCGCGAGTGATGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaac                                                               <  1:679095‑M1/143‑81 (MQ=255)                                                                                         aCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgcagtcacgaacgtaccgtagcgcgg                                                 >  2:84436‑M1/1‑49 (MQ=255)                                                                                         aCCGCCGATACCGCGAGTGATGGTGGACATTAAAATCCAGCTCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                 <  1:789890‑M1/143‑95 (MQ=255)                                                                                         aCCGCCGATACCGCGAGGGATGGGGGACATTAAAATCCAGAGCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                 <  1:992089‑M1/143‑95 (MQ=255)                                                                                          ccGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggagaacagcagcgcgacaaacgaagtcaccaacgtaccgtagaccggc                                                >  2:468181‑M1/1‑48 (MQ=255)                                                                                            gTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcac                                              <  1:1042648‑M1/143‑98 (MQ=255)                                                                                              cGATACCGCGTGTTATGGGGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcacca                                            <  2:1052711‑M1/143‑100 (MQ=255)                                                                                                   gcgcgaGTGACGGTGGCCCTTAAACACCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgcc                                       <  1:488571‑M1/142‑105 (MQ=255)                                                                                                   ccGCGAGTGATGGTGGACATTAAAATCCAAATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactaaccgggacgacgatcagcagcgcgaaaaacgcagtcaccaacgtacggtagaccggcaccaacgcc                                       >  2:212873‑M1/1‑39 (MQ=255)                                                                                                     gcgAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgcccg                                     >  1:158711‑M1/1‑37 (MQ=255)                                                                                                        aGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgattagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                  <  1:510507‑M1/143‑110 (MQ=255)                                                                                                        aGTGATGGTGGACATTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                  <  1:412771‑M1/143‑110 (MQ=255)                                                                                                          tGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtaga                                >  1:956421‑M1/1‑32 (MQ=255)                                                                                                               gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtgccgtagatcggcaccagcgccccgtagatatcg                           >  1:148184‑M1/1‑27 (MQ=255)                                                                                                               gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcg                           >  2:939356‑M1/1‑27 (MQ=255)                                                                                                               gtggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaaggcaccaacggaccgtagatcggcaccagcgccccggagatgtcg                           >  2:784617‑M1/1‑27 (MQ=255)                                                                                                                 ggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcggcagcgcgataaacgaagtcaccaacgtaccgcagataggcaccagagccccgtagacagcgta                         >  2:92620‑M1/1‑25 (MQ=255)                                                                                                                  gACATTAAAATCCAGATCAACCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttc                        >  1:250987‑M1/1‑24 (MQ=255)                                                                                                                     aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgctcagcagcccgataaacgaagtcaccaacgtaccgtagcgaggcacccgcgccccgtagatatcgttcggt                     >  1:599835‑M1/1‑21 (MQ=255)                                                                                                                     aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                     >  1:398369‑M1/1‑21 (MQ=255)                                                                                                                           gtcCAGATCAGCAAGgcgcaagctcaggcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcc               <  1:982752‑M1/142‑129 (MQ=255)                                                                                                                             ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactccccgggccagcaaatagcagcgcgataaccgaagtcaccaacgtaccgtagatcggcgccagcgccccgtaaataacgctcggagccagccc             >  2:780566‑M1/1‑13 (MQ=255)                                                                                                                              cAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccac            <  1:772486‑M1/143‑132 (MQ=255)                                                                                                                                gATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactc          >  1:953673‑M1/1‑10 (MQ=255)                                                                                                                                  tCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctt        <  1:458454‑M1/143‑136 (MQ=255)                                                                                                                                    aGCCAGgcgcaagctcagtccgggacagggcgataccgaagctcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttg      <  2:1183175‑M1/143‑138 (MQ=255)                                                                                                                                       cAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggccatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagagatcgttcggtgcatcgcactctttggtc   >  1:251356‑M1/1‑3 (MQ=255)                                                                                                                                        aGgagcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgattagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagagctccgtagatatcgttcggtgcatccgactctttggtac  >  1:619033‑M1/1‑2 (MQ=255)                                                                                                                                         |                                                                                                                                              ACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCCT  >  NC_000913/3909175‑3909452

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