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

ACACTTCGTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTG  >  NC_000913/3682012‑3682496                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             acacTTCGTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgc                                                                                                                                                                                                                                            <  1:376282‑M1/251‑3 (MQ=255)  cacTTCGTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcg                                                                                                                                                                                                                                           >  1:282546‑M1/1‑248 (MQ=255)   acTTCGTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgc                                                                                                                                                                                                                                          <  1:109847‑M1/251‑5 (MQ=255)           ttagacatcaagaacatGAAAATGTCCCGAAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagctt                                                                                                                                                                                                                                  >  1:121302‑M1/9‑239 (MQ=255)               aCCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTATCAACGACAgcgcacgagcttctga                                                                                                                                                                                                                              <  2:37548‑M1/251‑17 (MQ=255)                           tGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtc                                                                                                                                                                                                                  >  1:15164‑M1/1‑223 (MQ=255)                                    cactgaagaagtgaAAGAAACCTGAACTCACCCAGGTGGTAAGGTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCCCGCGTGGAAGGGAGTGCGGGCATAAGTTATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCGGCACATCGTCCAGTTTTTTGTGGGCCAGGTCTTTCACCCACTTAGCAACGACAgggcactaggttctgtcataaattggtcgattccgat                                                                                                                                                                                                         <  2:337475‑M1/247‑38 (MQ=255)                                       ggAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggat                                                                                                                                                                                                      >  2:94547‑M1/1‑211 (MQ=255)                                                gAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgcca                                                                                                                                                                                              <  1:94547‑M1/250‑49 (MQ=255)                                                                  cccAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatg                                                                                                                                                                           >  1:304552‑M1/1‑184 (MQ=255)                                                                            aaGTTCAGGAACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgc                                                                                                                                                                  >  2:303040‑M1/1‑174 (MQ=255)                                                                                                                      cGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacc                                                                                                                       <  1:215752‑M1/251‑120 (MQ=255)                                                                                                                         ggAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCCCCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccct                                                                                                                    <  2:53666‑M1/251‑123 (MQ=255)                                                                                                                                      ggCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgccccttt                                                                                                         >  1:5018‑M1/1‑116 (MQ=255)                                                                                                                                          tAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttaga                                                                                                      >  2:249518‑M1/1‑112 (MQ=255)                                                                                                                                             gTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaag                                                                                                <  2:304552‑M1/251‑143 (MQ=255)                                                                                                                                                          gagGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgatt                                                                                   >  1:65977‑M1/1‑96 (MQ=255)                                                                                                                                                            ggATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaa                                                                                 <  1:55840‑M1/251‑158 (MQ=255)                                                                                                                                                             gATAATACGTGAGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCAACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaac                                                                                >  2:48902‑M1/1‑93 (MQ=255)                                                                                                                                                               tAATTCGTGCGTTTTGCCATCAGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaac                                                                                >  2:12443‑M1/1‑91 (MQ=255)                                                                                                                                                                 aTTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTTGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcg                                                                            <  2:228114‑M1/251‑163 (MQ=255)                                                                                                                                                                         gTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttcgaagaaagcagcaacacgattaacagcgtgatttgg                                                                    >  2:203633‑M1/1‑81 (MQ=255)                                                                                                                                                                         gTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttgg                                                                    >  2:14802‑M1/1‑81 (MQ=255)                                                                                                                                                                                 aTCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagcaactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgat                                                            >  2:219320‑M1/1‑73 (MQ=255)                                                                                                                                                                                            gATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgac                                                 >  2:376760‑M1/1‑62 (MQ=255)                                                                                                                                                                                                        caTCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcct                                     <  2:65977‑M1/251‑202 (MQ=255)                                                                                                                                                                                                           cGTCCAGTTTTTTGTGGACCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctggg                                  >  2:153828‑M1/1‑47 (MQ=255)                                                                                                                                                                                                                      ttGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacacc                       >  1:60709‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                        aCCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtccggtat     >  2:374215‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                           cACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtcagggatatc  >  2:26922‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                                  cAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtg                                                                          >  2:331768‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                                  cAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtg                                                                          <  1:331768‑M1/172‑165 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             ACACTTCGTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTG  >  NC_000913/3682012‑3682496

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: ‑