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

AGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTT  >  NC_000913/3682024‑3682487                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    aGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttc                                                                                                                                                                                                                        <  1:203687‑M1/250‑14 (MQ=255)             cATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagct                                                                                                                                                                                                                          <  2:80153‑M1/236‑12 (MQ=255)             cATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagct                                                                                                                                                                                                                          >  1:80153‑M1/1‑225 (MQ=255)               tGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGAGATAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCAGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCAGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTCGCAACGACAgcgaacgagcatctgacataaagcggtc                                                                                                                                                                                                         >  2:165068‑M1/1‑223 (MQ=255)                  aaaTGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgat                                                                                                                                                                                                      >  2:351182‑M1/1‑220 (MQ=255)                    aTGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgatac                                                                                                                                                                                                    <  1:152653‑M1/251‑34 (MQ=255)                           tgactcctttccctaaacctgaactcaCCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCCGGCTAATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACGCCCTTCTAAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggat                                                                                                                                                                                             >  1:9896‑M1/14‑211 (MQ=255)                                             cTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctac                                                                                                                                                                           >  1:211478‑M1/1‑193 (MQ=255)                                                                                                        cTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGAGTTTTGCCATCCGGCGCCCGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgtt                                                                                                                 >  1:251884‑M1/1‑134 (MQ=255)                                                                                                          cGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacc                                                                                                              >  1:401105‑M1/1‑132 (MQ=255)                                                                                                            gggAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccc                                                                                                            <  2:357732‑M1/251‑122 (MQ=255)                                                                                                              gAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctg                                                                                                          >  1:81599‑M1/1‑128 (MQ=255)                                                                                                                    aGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgccc                                                                                                    <  1:39149‑M1/251‑130 (MQ=255)                                                                                                                    aGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgccc                                                                                                    <  1:163240‑M1/251‑130 (MQ=255)                                                                                                                       gCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctt                                                                                                 <  1:296857‑M1/251‑133 (MQ=255)                                                                                                                                 gTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaa                                                                                         >  1:58643‑M1/1‑109 (MQ=255)                                                                                                                                   gatgaGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagca                                                                                     <  1:125871‑M1/251‑145 (MQ=255)                                                                                                                                   gatgaGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagc                                                                                      >  2:114491‑M1/1‑107 (MQ=255)                                                                                                                                                ggATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgatta                                                                         >  1:112926‑M1/1‑94 (MQ=255)                                                                                                                                                                              aCGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctg                                          <  1:213525‑M1/251‑188 (MQ=255)                                                                                                                                                                                          cacaTCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggc                              >  1:47446‑M1/1‑52 (MQ=255)                                                                                                                                                                                          cacaTCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacagaaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggactatatccatctgacggttagtggc                              <  2:21483‑M1/251‑200 (MQ=255)                                                                                                                                                                                                  ccAGTTTGTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactccccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcg                       <  2:97894‑M1/250‑207 (MQ=255)                                                                                                                                                                                                   ccgttTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgcgcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccacgacccgttacccctgccgcccctgtagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgat                     <  2:202850‑M1/249‑209 (MQ=255)                                                                                                                                                                                                       tttttGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagt                                 >  1:260728‑M1/1‑39 (MQ=255)                                                                                                                                                                                                       tttttGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcgtctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagt                                 <  2:260728‑M1/235‑197 (MQ=255)                                                                                                                                                                                                       tttttGTGGTCCAGTTATTTAACCCACTTAGCAACGACAgcgcacgagcttctgacataacgcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgttgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaaaacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcagtaaac                 >  2:279992‑M1/1‑39 (MQ=255)                                                                                                                                                                                                             ggggcaagttcttTCACCCACTTAGCAACGACAgcgcacgagctgctgacataaagcggtcgataccgaggatcagcgcaagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgcc           <  2:81599‑M1/245‑219 (MQ=255)                                                                                                                                                                                                                       cgttcaCCCATGTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaagggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttccccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtc  <  2:404352‑M1/248‑228 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    AGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTT  >  NC_000913/3682024‑3682487

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