Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,054,149	Δ14 bp	coding (679‑692/1326 nt)	pepP ←	proline aminopeptidase P II

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3054149	3054162	14	25 [0]	[0] 24	pepP	proline aminopeptidase P II

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3054148	0 (0.000)	24 (0.520)	18/276	1.4	100%	coding (693/1326 nt)	pepP	proline aminopeptidase P II
	?	NC_000913	3054163 =	0 (0.000)					coding (678/1326 nt)	pepP	proline aminopeptidase P II

AGCGTAACCTTTGTATTCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTA  >  NC_000913/3054031‑3054290                                                                                                                                     |                                                                                                                                cgggtaatatcgccagcgtaacctttgtattcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCAc                                                                                                                          >  1:177832‑M2/133‑140 (MQ=255)             ccagcgtaacctttgtattcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTaa                                                                                                              >  2:462821‑M2/121‑140 (MQ=255)                 cgtaacctttgtattcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAattt                                                                                                          >  2:273971‑M2/117‑139 (MQ=255)                    aacctttgtattcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCgtg                                                                                                       <  2:43214‑M2/27‑1 (MQ=255)                    aacctttgtattcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCgtg                                                                                                       <  2:305793‑M2/27‑1 (MQ=255)                              ttcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCg                                                                                             <  2:684841‑M2/37‑1 (MQ=255)                              ttcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCg                                                                                             <  1:625464‑M2/37‑1 (MQ=255)                               tcacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGc                                                                                            <  1:462821‑M2/38‑1 (MQ=255)                                cacaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGcc                                                                                           <  2:421350‑M2/39‑1 (MQ=255)                                 acaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCgtg                                                                                                       <  1:573495‑M2/27‑1 (MQ=255)                                 acaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCgtg                                                                                                       >  2:573495‑M2/101‑127 (MQ=255)                                 acaacccgcgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCt                                                                                          >  2:580881‑M2/101‑140 (MQ=255)                                         cgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCCGATg                                                                                  >  2:45204‑M2/93‑140 (MQ=255)                                         cgtcaatcaacaccaggtcgccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCCCCTTCCAGATg                                                                                  >  2:459350‑M2/93‑140 (MQ=255)                                                             ccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGAc                                                              >  2:189535‑M2/73‑140 (MQ=255)                                                             ccgtcgcgcatttcacactcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGAc                                                              <  1:203574‑M2/68‑1 (MQ=255)                                                                              ctcgttttcggtgtagtgcagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGcc                                             <  1:189535‑M2/85‑1 (MQ=255)                                                                                                cagaatgcagccgttttccccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCg                           <  1:121244‑M2/103‑1 (MQ=255)                                                                                                cagaatgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCg                           <  1:296866‑M2/103‑1 (MQ=255)                                                                                                     tgcagccgttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAAATCGTGGTGAATTTCGCCTTCCAGATGGTCGCCCAACATTACCGGACGGCATTTTTCCATCACCCGTGTATTTGGCAtggcggaggt                      >  2:296866‑M2/33‑136 (MQ=255)                                                                                                             ttttcaccgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCcg              <  1:580881‑M2/116‑1 (MQ=255)                                                                                                                    cgctgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCgcggcg       <  1:92630‑M2/123‑1 (MQ=255)                                                                                                                       tgccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAg    <  1:763749‑M2/126‑1 (MQ=255)                                                                                                                         ccgacaatggtgGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTa  <  1:181600‑M2/128‑1 (MQ=255)                                                                                                                                     |                                                                                                                                AGCGTAACCTTTGTATTCACAACCCGCGTCAATCAACACCAGGTCGCCGTCGCGCATTTCACACTCGTTTTCGGTGTAGTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTA  >  NC_000913/3054031‑3054290

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑