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] 29	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 (1.350)	26/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

GTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCG  >  NC_000913/3682019‑3682500                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 gtgtTAGACCGCAAGAACATGAAAATGACACGGAAGAAGTGAAATAAACCTGAAATCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacg                                                                                                                                                                                                                                           >  2:169398‑M1/1‑242 (MQ=255)     tAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttc                                                                                                                                                                                                                                     >  1:81790‑M1/1‑238 (MQ=255)                  cATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcgg                                                                                                                                                                                                                        <  2:62835‑M1/251‑27 (MQ=255)                  cATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGCGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGCGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCACGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcgg                                                                                                                                                                                                                        >  2:120873‑M1/1‑225 (MQ=255)                             cACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccg                                                                                                                                                                                                               >  1:66544‑M1/1‑214 (MQ=255)                                                  cTGAACTCACCCAGGTGGTAAGTGCAGGCACAGGGGTCAATTATGCGCACACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGCTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataacgcggtcgataccgaggatcagcgccagacccgctac                                                                                                                                                                                        >  2:90980‑M1/1‑193 (MQ=255)                                                            ccAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGAACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaantgg                                                                                                                                                                              >  2:201717‑M1/1‑183 (MQ=255)                                                                     aaGTTCAGGCACAGGGGTCAATTATGCTCAAACCCCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTTAGTGGATACTTCGTGCGTTTTGCCATCCGTCGCACGATTATTCAGCACATCGTCCATTTTGTTGTGGTCCAGTGCTTTCACTCACTTAGCAACGACAgcgcacgagcttctgtcgtaaagcgttcggtaccgaggatcagcgccagacccgctaccggcaaatggcctaccgca                                                                                                                                                                     <  1:26595‑M1/251‑78 (MQ=255)                                                                          cAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCTCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagag                                                                                                                                                                <  1:168622‑M1/251‑83 (MQ=255)                                                                              cacaGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtc                                                                                                                                                            >  1:151156‑M1/1‑165 (MQ=255)                                                                                     gTCAATTATGCGCAAACGCCCGCACTCGGGGAAGGGAGTGCGTGCATAAGTGATGAGATTAAGAGGATTATTCGTTCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagatcgtcgccgccc                                                                                                                                                     <  2:160428‑M1/251‑94 (MQ=255)                                                                                      tcacttatgcGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGTACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTAACCCACTTAGCACCGACAgcgcacgagcttctgacataaagcggccgataccgaggatcagcgccagacccgctaccgtcaaatggcccaccgcagagagcgtcgccgcc                                                                                                                                                      >  2:65750‑M1/5‑157 (MQ=255)                                                                                                            aCTCGGGGAAGGGCGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCGAGTTTTTTGTGGTACAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccacagcagagagcgtcgccgccagcacgataaagcctctacccgtt                                                                                                                              >  2:96204‑M1/1‑135 (MQ=255)                                                                                                                                  aTAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttaga                                                                                                          >  1:103701‑M1/1‑113 (MQ=255)                                                                                                                                                         aaTTCGTTCTTTTGGCCATCCGGCGCACGAGTATTCAGCACATCGTCCAGTTTTTTTTGGTCCAGTTCTTTCACCCCCTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacag                                                                                  <  2:177468‑M1/250‑161 (MQ=255)                                                                                                                                                             cGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtg                                                                              >  1:49904‑M1/1‑86 (MQ=255)                                                                                                                                                              gTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGATCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctcccgtcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgat                                                                            <  2:39982‑M1/251‑167 (MQ=255)                                                                                                                                                                       gCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgaacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggac                                                                   <  1:36283‑M1/251‑176 (MQ=255)                                                                                                                                                                             cGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatc                                                             >  2:7309‑M1/1‑70 (MQ=255)                                                                                                                                                                              ggCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACGTAGCAACGCCAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatcc                                                            <  1:90980‑M1/251‑183 (MQ=255)                                                                                                                                                                               gcgcACGATTATTCAGCCCATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggtccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatcca                                                           <  1:120873‑M1/251‑184 (MQ=255)                                                                                                                                                                                         attCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaccaacgtgctttgcgggacgatatccatcttcctgtt                                                 >  2:46809‑M1/1‑58 (MQ=255)                                                                                                                                                                                         attCAGCACATCGTCCAGTTTTTTGTGGCCCAGTTCTTTCACCCAATTAACAACGACAgcgtacgatcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagcgcgtcgccgccagcacaataaagccactacccgttacccctgccgcccctttagacgaaagcagcaacacacttatcagcgttatttggtggccgatatccatatagatgtt                                                 >  1:27130‑M1/1‑58 (MQ=255)                                                                                                                                                                                             aGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtg                                             <  1:201609‑M1/251‑198 (MQ=255)                                                                                                                                                                                                        cAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagatcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaccacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgat                                  <  1:427‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                       cAGTTCTTTCACCCACTTAGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgcttaacagcgtgatttggtggacgatatccatctgactgttagtgtcctggtcgataaacaccgccgccat                   >  2:42587‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                    cccttAGCAACGACGgcccacgagcttctgacataaagcggtcgatcccgaggatcagcgcgagacccgcttccggcaaatggcccaccgcagagagcgtcgccgccatcacgataaagcctctacccgttccccctgccgcccctttaaaagaaagcagcaacacgattaacagcgtgatttggtggacgatgtccatcggactgttagtggcctgggcgataaacaccgccgccattgtcaggtatatc      <  1:96204‑M1/249‑237 (MQ=255)                                                                                                                                                                                                                                      ataagaaacgacagctcacgagcttctgacagaaagcggtcgataccggggatcagcgccaaccccgctcccgccaaatgccccaccgcagagggcgtcgccgccagcacgataaagccactacccgttacccctgccgccccttgagaagaaagcagcaacacgattaacagcgtgattgggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtcaggtatatcga    <  2:103701‑M1/245‑239 (MQ=255)                                                                                                                                                                                                                                         aGCAACGACAgcgcacgagcttctgacataaagcggtcgataccgaggatcagcgccagacccgctaccggcaaatggcccaccgcagagagcgtcgccgccagcacgataaagccactacccgttacccctgccgcccctttagaagaaagcagcaacacgattaacagcgtgatttggtggacgatatccatctgactgttagtggcctgggcgataaacaccgccgccattgtcaggtatatcgatg  <  1:193966‑M1/250‑241 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 GTGTTAGACCGCAAGAACATGAAAATGTCACGGAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCG  >  NC_000913/3682019‑3682500

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