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	29 [0]	[0] 30	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)	29 (0.750)	26/470	0.6	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

TGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATC  >  NC_000913/3682058‑3682541                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       agaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTc                                                                                                                                                                                                                                            <  1:203687‑M2/13‑1 (MQ=255)             catgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCt                                                                                                                                                                                                                                              <  2:80153‑M2/11‑1 (MQ=255)             catgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCt                                                                                                                                                                                                                                              >  1:80153‑M2/226‑236 (MQ=255)               tgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggagataagttcaggcacaggggtcaattatgcgcaaacaccagcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatcaggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttcgcaacgacaGCGAACGAGCATCTGACATAAAGCGGTc                                                                                                                                                                                                                             >  2:165068‑M2/224‑251 (MQ=255)                  aaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGAt                                                                                                                                                                                                                          >  2:351182‑M2/221‑251 (MQ=255)                    atgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATAc                                                                                                                                                                                                                        <  1:152653‑M2/33‑1 (MQ=255)                           tgactcctttccctaaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgccggctaattcagcacatcgtccagttttttgtggtccagttctttcacgcccttctaaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGAt                                                                                                                                                                                                                 >  1:9896‑M2/212‑251 (MQ=255)                                             ctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTAc                                                                                                                                                                                               >  1:211478‑M2/194‑251 (MQ=255)                                                                                                        ctcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgagttttgccatccggcgcccgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGtt                                                                                                                                     >  1:251884‑M2/135‑250 (MQ=255)                                                                                                          cggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTAcc                                                                                                                                  >  1:401105‑M2/133‑251 (MQ=255)                                                                                                            gggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTAcccc                                                                                                                                <  2:357732‑M2/121‑1 (MQ=255)                                                                                                              gaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTg                                                                                                                              >  1:81599‑M2/129‑251 (MQ=255)                                                                                                                    agtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGccc                                                                                                                        <  1:39149‑M2/129‑1 (MQ=255)                                                                                                                    agtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGccc                                                                                                                        <  1:163240‑M2/129‑1 (MQ=255)                                                                                                                       gcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCtt                                                                                                                     <  1:296857‑M2/132‑1 (MQ=255)                                                                                                                                 gtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAgaagaa                                                                                                             >  1:58643‑M2/110‑249 (MQ=255)                                                                                                                                   gatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAagca                                                                                                         <  1:125871‑M2/144‑1 (MQ=255)                                                                                                                                   gatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAagc                                                                                                          >  2:114491‑M2/108‑250 (MQ=255)                                                                                                                                                ggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTa                                                                                             >  1:112926‑M2/95‑250 (MQ=255)                                                                                                                                                                              acgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTg                                                              <  1:213525‑M2/187‑1 (MQ=255)                                                                                                                                                                                          cacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGc                                                  >  1:47446‑M2/53‑251 (MQ=255)                                                                                                                                                                                          cacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACAGAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACTATATCCATCTGACGGTTAGTGGc                                                  <  2:21483‑M2/199‑1 (MQ=255)                                                                                                                                                                                                  ccagtttgttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTCCCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCg                                           <  2:97894‑M2/206‑1 (MQ=255)                                                                                                                                                                                                   ccgttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGCGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACGACCCGTTACCCCTGCCGCCCCTGTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGAt                                         <  2:202850‑M2/208‑1 (MQ=255)                                                                                                                                                                                                       tttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGt                                                     >  1:260728‑M2/40‑235 (MQ=255)                                                                                                                                                                                                       tttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCGTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGt                                                     <  2:260728‑M2/196‑1 (MQ=255)                                                                                                                                                                                                       tttttgtggtccagttatttaacccacttagcaacgacaGCGCACGAGCTTCTGACATAACGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTTGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAAAACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGgcctgggcagtaaac                                     >  2:279992‑M2/40‑244 (MQ=255)                                                                                                                                                                                                             ggggcaagttctttcacccacttagcaacgacaGCGCACGAGCTGCTGACATAAAGCGGTCGATACCGAGGATCAGCGCAAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAccgcc                               <  2:81599‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                       cgttcacccatgtagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAAGGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTCCCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTc                      <  2:404352‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                          gccggcgcACGAGCGTCTGACATAAGGCTGTCGATACCGAGGTTCAGCGCCAGACCCGCTACCGGCAAATGGCCCGCCGCAGAGAGCGTCCCCGCCAGCACGATAAGGCCGCTACCCGTTGCCCCGGCCGCCCCTTTAGAAGAAAGCAGCAGCACGATTAACAGCGGGATTTGGTGGACGATATCCATCGGACTGTTATGGGTCTGGGCGATAAGCACCGCCGCTATTGTCATGTATAGTGATGTGCCATc  <  2:251884/247‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       TGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCATC  >  NC_000913/3682058‑3682541

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