Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,176,353	(GCA)4→3	coding (479‑481/828 nt)	cpdA ←	3',5' cAMP phosphodiesterase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3176344	3176346	3	29 [0]	[0] 31	cpdA	3',5' cAMP phosphodiesterase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3176347 =	0 (0.000)	29 (0.860)	24/242	0.6	100%	coding (487/828 nt)	cpdA	3',5' cAMP phosphodiesterase
	?	NC_000913	= 3176343	0 (0.000)					coding (491/828 nt)	cpdA	3',5' cAMP phosphodiesterase

TCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCT  >  NC_000913/3176208‑3176493                                                                                                                                            |                                                                                                                                                   agatccagctcctgatgaatatgaccgcacagcaagtatttgacgtgcggaaacttcgccagcacggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCaa                                                                                                                                           <  4:125289‑M1/10‑1 (MQ=255)                  aatatgaccgcacagcaagtatttgacgtgcggaaacttcgccagcacggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGgc                                                                                                                           >  1:106916‑M1/123‑148 (MQ=255)                                       tttgacgtgcggaaacttcgccagcacggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGtt                                                                                                              >  8:21095‑M1/102‑140 (MQ=255)                                       tttgacgtgcggaaacttcgccagcacggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGtt                                                                                                              <  7:21095‑M1/39‑1 (MQ=255)                                                         cgccagcacggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGaaa                                                                                     >  3:125289‑M1/84‑147 (MQ=255)                                                                acggtatccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCgctcag                                                                            <  2:87222‑M1/73‑1 (MQ=255)                                                                     atccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCTGCCa                                                                                                 >  6:122110‑M1/72‑123 (MQ=255)                                                                     atccagttcgcccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCa                                                                                                 <  5:122110‑M1/52‑1 (MQ=255)                                                                               cccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGgc                                                                                                                           <  6:84485‑M1/26‑1 (MQ=255)                                                                               cccgcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGgc                                                                                                                           >  5:84485‑M1/62‑87 (MQ=255)                                                                                  gcgttacgcagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGc                                                           <  5:35619‑M1/90‑1 (MQ=255)                                                                                          cagactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCt                                                  >  1:87222‑M1/51‑149 (MQ=255)                                                                                           agactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAac                                                      <  7:115621‑M1/95‑1 (MQ=255)                                                                                           agactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAac                                                      >  8:115621‑M1/50‑144 (MQ=255)                                                                                            gactgtgttgatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                                                >  6:35619‑M1/49‑149 (MQ=255)                                                                                                     gatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTg                                                 <  5:163604‑M1/100‑1 (MQ=255)                                                                                                     gatcgagccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTg                                                 >  6:163604‑M1/40‑139 (MQ=255)                                                                                                            ccaactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGa                                <  2:121121‑M1/117‑1 (MQ=255)                                                                                                              aactacagcccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAtt                              <  2:91005‑M1/119‑1 (MQ=255)                                                                                                              aactacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGAtt                              <  8:74818‑M1/119‑1 (MQ=255)                                                                                                                 tacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                                                <  6:42920‑M1/101‑1 (MQ=255)                                                                                                                 tacaacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTgg                                                >  5:42920‑M1/28‑128 (MQ=255)                                                                                                                   caacccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGc                           >  1:121121‑M1/26‑147 (MQ=255)                                                                                                                      cccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGgccaatt                                                                                                              >  4:66069‑M1/23‑58 (MQ=255)                                                                                                                      cccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGgccaatt                                                                                                              <  3:66069‑M1/39‑4 (MQ=255)                                                                                                                      cccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCa                         >  2:25016‑M1/23‑146 (MQ=255)                                                                                                                      cccgcaggtagcggatgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCa                         <  1:25016‑M1/124‑1 (MQ=255)                                                                                                                                    atgatgatGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAaca        >  1:91005‑M1/9‑149 (MQ=255)                                                                                                                                          atGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGCCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCt  <  2:191683‑M1/147‑1 (MQ=255)                                                                                                                                            GCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGc                           <  2:47784‑M1/122‑1 (MQ=255)                                                                                                                                            GCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGc                           >  1:47784‑M1/1‑122 (MQ=255)                                                                                                                                            |                                                                                                                                                   TCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCT  >  NC_000913/3176208‑3176493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑