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] 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)	30 (1.090)	29/486	0.1	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

TAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAC  >  NC_000913/3682023‑3682508                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         tAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttc                                                                                                                                                                                                                                             >  1:168527‑M1/1‑238 (MQ=255)  aGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttc                                                                                                                                                                                                                                             <  1:278703‑M1/250‑14 (MQ=255)           gAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaag                                                                                                                                                                                                                                   <  1:87433‑M1/251‑24 (MQ=255)                            ggAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggat                                                                                                                                                                                                                  >  1:172146‑M1/1‑211 (MQ=255)                                    tGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgcc                                                                                                                                                                                                           >  1:111843‑M1/1‑203 (MQ=255)                                     gAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgcc                                                                                                                                                                                                           >  1:274185‑M1/1‑202 (MQ=255)                                         tAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagaccc                                                                                                                                                                                                     >  1:202632‑M1/1‑198 (MQ=255)                                                               gtAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccg                                                                                                                                                                               <  1:27950‑M1/251‑76 (MQ=255)                                                                   gTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgc                                                                                                                                                                              >  1:275620‑M1/1‑172 (MQ=255)                                                                                     aTTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcac                                                                                                                                                         >  1:226792‑M1/1‑154 (MQ=255)                                                                                      ttATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACa                                                                                                                                                                                                                                                          <  2:278103/153‑1 (MQ=255)                                                                                      ttATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACa                                                                                                                                                                                                                                                          >  1:278103/1‑153 (MQ=255)                                                                                                acacCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccac                                                                                                                                              >  1:70181‑M1/1‑143 (MQ=255)                                                                                                                            gCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggccaaccgcagagagcgtcgccgccagcacgataaagccaatacccgttacccctgccgccccctaagc                                                                                                                  >  2:173998‑M1/1‑115 (MQ=255)                                                                                                                            gCATAAGTGATGAGATTAAGAGGAGAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccggtacccctgccgcccctttaga                                                                                                                  <  2:198069‑M1/251‑137 (MQ=255)                                                                                                                             cATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttaga                                                                                                                  >  2:155918‑M1/1‑114 (MQ=255)                                                                                                                                                    tAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACCgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcgcaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaac                                                                                            >  2:48822‑M1/1‑91 (MQ=255)                                                                                                                                                                    ccATCCGGCGCACGATTTTTCAGCCCATCGTCCCGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagtgtgatttggtggacg                                                                          <  1:92289‑M1/251‑177 (MQ=255)                                                                                                                                                                        cctgcgcacgATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatat                                                                      <  1:59865‑M1/248‑181 (MQ=255)                                                                                                                                                                          ggCGCACGATTATCCAGCACATCGTCCAGTTTTTTGTGTGCCAGTTCTTTCACCCACTTAGCCACGACAgcgcacgagcttctgacatacagcggtcgatacagaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccgctacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttgcttgacgatatcc                                                                    >  2:225527‑M1/1‑69 (MQ=255)                                                                                                                                                                             gcACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatc                                                                 >  1:178650‑M1/1‑66 (MQ=255)                                                                                                                                                                              cACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatct                                                                >  1:216006‑M1/1‑65 (MQ=255)                                                                                                                                                                                cGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctg                                                               >  2:167907‑M1/1‑63 (MQ=255)                                                                                                                                                                                          gCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACTACAgcgcacgatcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactaccctttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtgg                                                    <  2:275620‑M1/251‑199 (MQ=255)                                                                                                                                                                                            acaTCGTACAGTTTTTTGTGGTACAGTTCTTTCACCCACTGAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcacatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtggtttggtggacgatctccatatgacagttagtggcc                                                  >  1:274457‑M1/1‑51 (MQ=255)                                                                                                                                                                                                   ccAGTTTTTTGTGGACCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcg                                            >  1:87574‑M1/1‑44 (MQ=255)                                                                                                                                                                                                    gagTTATTTGGGGTCCAGTTCTTGCACCCACTTAGCAACGACAgcgcaagagcttgtgacataaagcgggcgatcccgaggatcagcgccagaccctctaccggcaaatggcccaccgcagagggcgtcgccgccagcacgataaagccgctaccctttacccctgccgccccttgagaagaaagcagcaactcgatgaacagcgtgatttggtggacgatatccatctgactgttagtgtcctgggcgat                                          <  2:117164‑M1/250‑209 (MQ=255)                                                                                                                                                                                                      gTTTTTTGTGGTCCAGTTCTTTCCCCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggagcagcgccagacccgctcccggcaaatggccccccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataa                                        <  2:202632‑M1/251‑211 (MQ=255)                                                                                                                                                                                                         ttttGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaac                                      <  1:229700‑M1/250‑213 (MQ=255)                                                                                                                                                                                                                                        aaCGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtcaggtatatcgatgtgcc      >  1:218525‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                                        aaCGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgcaagacccgctaccggcaaatggcccaccgcagagagcgtagccgccagcacgataaagccactacccactacacctgccgcccctttagaagaaagcagcaacacgattaaaagagtgatttggtggaagatatccatctgactgttcgtggcctgcgagatcaacaccggcgccattgtcaggtatatcgctgtgcc      >  1:285232‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                                            aCAgcgcacgagtttctgacaaaaaggggtcgatcccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcggcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaataaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcttgggcgataaacaccgccgccattgtcaggtatatcgatgtgccatca  <  2:172146‑M1/251‑249 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         TAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAC  >  NC_000913/3682023‑3682508

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