Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,682,261	Δ34 bp	coding (1154‑1187/1287 nt)	dctA ←	C4 dicarboxylate/orotate:H(+) symporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3682261	3682294	34	32 [0]	[0] 32	dctA	C4 dicarboxylate/orotate:H(+) symporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3682260	0 (0.000)	32 (0.880)	29/488	0.4	100%	coding (1188/1287 nt)	dctA	C4 dicarboxylate/orotate:H(+) symporter
	?	NC_000913	3682295 =	0 (0.000)					coding (1153/1287 nt)	dctA	C4 dicarboxylate/orotate:H(+) symporter

TCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATC  >  NC_000913/3682046‑3682530                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            acacttcgtgttagaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacagc                                                                                                                                                                                                                                            <  1:376282‑M2/2‑1 (MQ=255)  cacttcgtgttagaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacagcg                                                                                                                                                                                                                                           >  1:282546‑M2/249‑251 (MQ=255)   acttcgtgttagaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacagcgc                                                                                                                                                                                                                                          <  1:109847‑M2/4‑1 (MQ=255)           ttagacatcaagaacatgaaaatgtcccgaaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCtt                                                                                                                                                                                                                                  >  1:121302‑M2/240‑251 (MQ=255)               accgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttatcaacgacaGCGCACGAGCTTCTGa                                                                                                                                                                                                                              <  2:37548‑M2/16‑1 (MQ=255)                           tgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTc                                                                                                                                                                                                                  >  1:15164‑M2/224‑251 (MQ=255)                                    cactgaagaagtgaaagaaacctgaactcacccaggtggtaaggtcaggcacaggggtcaattatgcgcaaacacccgcccgcgtggaagggagtgcgggcataagttatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcggcacatcgtccagttttttgtgggccaggtctttcacccacttagcaacgacaGGGCACTAGGTTCTGTCATAAATTGGTCGATTCCgat                                                                                                                                                                                                         <  2:337475‑M2/37‑2 (MQ=255)                                       ggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGAt                                                                                                                                                                                                      >  2:94547‑M2/212‑251 (MQ=255)                                                gaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCa                                                                                                                                                                                              <  1:94547‑M2/48‑1 (MQ=255)                                                                  cccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATg                                                                                                                                                                           >  1:304552‑M2/185‑251 (MQ=255)                                                                            aagttcaggaacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGc                                                                                                                                                                  >  2:303040‑M2/175‑250 (MQ=255)                                                                                                                      cggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTAcc                                                                                                                       <  1:215752‑M2/119‑1 (MQ=255)                                                                                                                         ggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcccccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCt                                                                                                                    <  2:53666‑M2/122‑1 (MQ=255)                                                                                                                                      ggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCttt                                                                                                         >  1:5018‑M2/117‑249 (MQ=255)                                                                                                                                          taagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTaga                                                                                                      >  2:249518‑M2/113‑248 (MQ=255)                                                                                                                                             gtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAag                                                                                                <  2:304552‑M2/142‑1 (MQ=255)                                                                                                                                                          gaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGAtt                                                                                   >  1:65977‑M2/97‑251 (MQ=255)                                                                                                                                                            ggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTaa                                                                                 <  1:55840‑M2/157‑1 (MQ=255)                                                                                                                                                             gataatacgtgagttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcaccaacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAAc                                                                                >  2:48902‑M2/94‑251 (MQ=255)                                                                                                                                                               taattcgtgcgttttgccatcaggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAAc                                                                                >  2:12443‑M2/92‑249 (MQ=255)                                                                                                                                                                 attcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgttgtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCg                                                                            <  2:228114‑M2/162‑1 (MQ=255)                                                                                                                                                                         gttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTCGAAGAAAGCAGCAACACGATTAACAGCGTGATTtgg                                                                    >  2:203633‑M2/82‑251 (MQ=255)                                                                                                                                                                         gttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTtgg                                                                    >  2:14802‑M2/82‑251 (MQ=255)                                                                                                                                                                                 atccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCAACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGat                                                            >  2:219320‑M2/74‑251 (MQ=255)                                                                                                                                                                                            gattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGAc                                                 >  2:376760‑M2/63‑251 (MQ=255)                                                                                                                                                                                                        catcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCt                                     <  2:65977‑M2/201‑1 (MQ=255)                                                                                                                                                                                                           cgtccagttttttgtggaccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTggg                                  >  2:153828‑M2/48‑251 (MQ=255)                                                                                                                                                                                                                      ttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAcc                       >  1:60709‑M2/37‑251 (MQ=255)                                                                                                                                                                                                                                        acccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCCGGtat     >  2:374215‑M2/19‑251 (MQ=255)                                                                                                                                                                                                                                           cacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGGATATc  >  2:26922‑M2/16‑251 (MQ=255)                                                                                                                                                                                                                                                  caacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTg                                                                          >  2:331768‑M2/9‑172 (MQ=255)                                                                                                                                                                                                                                                  caacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTg                                                                          <  1:331768‑M2/164‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            TCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATC  >  NC_000913/3682046‑3682530

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