Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,614,614	Δ3 bp	coding (681‑683/849 nt)	focB →	putative formate transporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2614614	2614616	3	30 [0]	[0] 28	focB	putative formate transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2614613	0 (0.000)	28 (0.550)	23/232	0.9	100%	coding (680/849 nt)	focB	putative formate transporter
	?	NC_000913	2614617 =	0 (0.000)					coding (684/849 nt)	focB	putative formate transporter

TTGCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGT  >  NC_000913/2614479‑2614747                                                                                                                                           |                                                                                                                                   ttattgcgggcgttcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccc                                                                                                                                    >  2:910254‑M2/139‑139 (MQ=255)     tgcgggcgttcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCt                                                                                                                                >  1:712209‑M2/135‑139 (MQ=255)         ggcgttcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTg                                                                                                                            >  2:486122‑M2/131‑139 (MQ=255)                   atgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGgcg                                                                                                                  <  1:742977‑M2/19‑1 (MQ=255)                        acaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAg                                                                                                             <  1:208388‑M2/24‑1 (MQ=255)                          aaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTa                                                                                                           <  1:910254‑M2/26‑1 (MQ=255)                                            ttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCgg                                                                                         <  2:712209‑M2/44‑1 (MQ=255)                                            ttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCgg                                                                                         <  2:298708‑M2/44‑1 (MQ=255)                                                   tcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGAc                                                                                  >  2:101919‑M2/89‑139 (MQ=255)                                                          gtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                           <  1:613505‑M2/58‑1 (MQ=255)                                                                cgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCAtttt                                                                     >  2:561043‑M2/76‑139 (MQ=255)                                                                          tttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCa                                                           <  2:961372‑M2/74‑1 (MQ=255)                                                                                cactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGc                                                     <  1:383561‑M2/80‑1 (MQ=255)                                                                                 actgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCt                                                    <  1:561043‑M2/81‑1 (MQ=255)                                                                                           caatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                          <  1:786440‑M2/91‑1 (MQ=255)                                                                                           caatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                          <  1:610911‑M2/91‑1 (MQ=255)                                                                                             atttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCt                                        >  2:47938‑M2/47‑139 (MQ=255)                                                                                             atttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCt                                        <  1:199440‑M2/93‑1 (MQ=255)                                                                                             atttgtttgtgattccgttcgccattgccattcgccatttcgacccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCt                                        >  2:7842‑M2/47‑139 (MQ=255)                                                                                              tttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTg                                       <  2:589443‑M2/94‑1 (MQ=255)                                                                                                               tcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATcggcgg                      <  2:698396‑M2/111‑1 (MQ=255)                                                                                                                 gccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtg                    <  1:898724‑M2/113‑1 (MQ=255)                                                                                                                 gccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtg                    <  1:770362‑M2/113‑1 (MQ=255)                                                                                                                 gccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtg                    <  1:274822‑M2/113‑1 (MQ=255)                                                                                                                     ttgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtgcggt                <  1:47938‑M2/117‑1 (MQ=255)                                                                                                                       gccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCggtgcggtgc              <  1:101919‑M2/119‑1 (MQ=255)                                                                                                                              gccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAg       <  1:7842‑M2/126‑1 (MQ=255)                                                                                                                                   ttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTAtgt  >  1:412034‑M2/9‑139 (MQ=255)                                                                                                                                           |                                                                                                                                   TTGCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTATGT  >  NC_000913/2614479‑2614747

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