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	42 [0]	[1] 43	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)	42 (1.040)	40/486	0.2	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

GAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCAT  >  NC_000913/3682051‑3682540                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      gcgtgttagaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCAc                                                                                                                                                                                                                                                  <  2:1216‑M2/6‑1 (MQ=255)   gtgttagaccgcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggccaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtgcagttctttcacccactgagcaacgacaTCGCACGAg                                                                                                                                                                                                                                               <  2:198714‑M2/9‑1 (MQ=255)             gcaagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACAt                                                                                                                                                                                                                                     >  2:71128‑M2/233‑251 (MQ=255)               aagaacatgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtgttccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATaa                                                                                                                                                                                                                                   <  1:411628‑M2/21‑1 (MQ=255)                 gaacatgaaaatgtcacggaagaagtgaaataaacatgaacacacccaggaggtcagatcaggcaaaggggtcaattatgcgcaaacacccccactcggggcagggagtgcgggcattagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcaaatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATaa                                                                                                                                                                                                                                   >  1:376499‑M2/229‑249 (MQ=255)                    catgaaaatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCgg                                                                                                                                                                                                                              >  1:340360‑M2/226‑251 (MQ=255)                          aatgtcacggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATa                                                                                                                                                                                                                        <  1:289331‑M2/32‑1 (MQ=255)                                  ggaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGAt                                                                                                                                                                                                                >  1:51473‑M2/212‑251 (MQ=255)                                   gaagaagtgaaataaacctgaactcacccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctgtcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATc                                                                                                                                                                                                               <  2:70668‑M2/41‑1 (MQ=255)                                                            acccaggtggtaagttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAAt                                                                                                                                                                                      >  1:261001‑M2/186‑251 (MQ=255)                                                                         gttcaggcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGc                                                                                                                                                                            >  1:265741‑M2/173‑248 (MQ=255)                                                                               gcacaggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGt                                                                                                                                                                   >  2:78854‑M2/167‑251 (MQ=255)                                                                                  caggggtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTcgc                                                                                                                                                                >  1:307227‑M2/164‑251 (MQ=255)                                                                                      ggtcaattatgcgcaaacacccgcactcggggaagggagagcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgaatattcagcacatcgtctagttttttgtgggccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGCTCAGCGCCAGACCCGCTACCGGCAAATGGCCAACCGCAGAGAGCGTCgccgcc                                                                                                                                                            >  1:224885‑M2/160‑251 (MQ=255)                                                                                       gtcaattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCgccgcc                                                                                                                                                            >  1:139547‑M2/159‑250 (MQ=255)                                                                                          tattatgcgcaaacacccgcactcggggaagggagtgcgggcataagtgatgagaataagaggataattcgtgcgttttgccatccggcgcacgattattcagcgcatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAgcc                                                                                                                                                        >  2:15200‑M2/156‑250 (MQ=255)                                                                                               tgcgcaaacaaccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtacagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGAt                                                                                                                                                    >  2:391352‑M2/151‑250 (MQ=255)                                                                                                      acacccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCAc                                                                                                                                            >  2:218505‑M2/144‑251 (MQ=255)                                                                                                        ccccgcactcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTa                                                                                                                                          <  1:242947‑M2/110‑1 (MQ=255)                                                                                                              actcggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGtt                                                                                                                                    <  1:268305‑M2/116‑1 (MQ=255)                                                                                                                 cggggaatggagtgcgtgcatcagtgatgacattaagaggatacttcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtcgagttctttcacccccttatcaacgtcaGCGCCCGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCTTCGCATAGAGCGTCGCCGCCAGCACGATAATGCCACTACCCGTTAc                                                                                                                                  <  2:194864‑M2/118‑1 (MQ=255)                                                                                                                  ggggaagggagtgcgggcataagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTAccc                                                                                                                                >  2:80887‑M2/132‑251 (MQ=255)                                                                                                                           aggncgggcataagtgatgagattaagaggataagtcgtgggttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGTTACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGACCGTCCCCGCCAGCACGATAGTGCCACTATCCGTTACCACtaccgcca                                                                                                                       <  2:376499‑M2/129‑2 (MQ=255)                                                                                                                                     taagtgatgaaattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgccaGCGCACGAGCTTCTGACATCAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGACGCCGCCAGCCCGATAAAGCCACTACCAGTTACCCCTGCCGCCCCTTTaga                                                                                                                >  2:245949‑M2/113‑248 (MQ=255)                                                                                                                                      aagtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAgaagaa                                                                                                            <  1:71128‑M2/140‑1 (MQ=255)                                                                                                                                       agtgatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAgaagaa                                                                                                            >  1:156996‑M2/111‑250 (MQ=255)                                                                                                                                          gatgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGTTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTCCCCGTTACCCCTGCCGCCCCTTTAGAAGAAagca                                                                                                        <  1:245949‑M2/144‑1 (MQ=255)                                                                                                                                           atgagattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCTGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAagcag                                                                                                       <  1:406940‑M2/145‑1 (MQ=255)                                                                                                                                              agattaagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAgcagca                                                                                                     >  2:29272‑M2/104‑250 (MQ=255)                                                                                                                                                   aagaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACg                                                                                                <  1:158257‑M2/152‑1 (MQ=255)                                                                                                                                                    agaggataattcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGAt                                                                                              >  2:176229‑M2/98‑251 (MQ=255)                                                                                                                                                             ttcgtgcgttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACAGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGt                                                                                     >  1:354778‑M2/89‑251 (MQ=255)                                                                                                                                                                    gttttgccatccggcgcacgattattcagcacatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTtgg                                                                              <  1:78854‑M2/170‑1 (MQ=255)                                                                                                                                                                                          tattcagcacatagtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGAAGATATCCATCTGACTGt                                                        >  2:405041‑M2/60‑251 (MQ=255)                                                                                                                                                                                           attcagcccatcgtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGtt                                                       >  1:324929‑M2/59‑251 (MQ=255)                                                                                                                                                                                               agcacatcgtccagttttttgtggtccagttatttcaccaacttagcaaagacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTg                                                   >  1:234604‑M2/55‑251 (MQ=255)                                                                                                                                                                                                       gtccagttttttgtggtccagttctttcacccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGc                                           <  2:139547‑M2/205‑1 (MQ=255)                                                                                                                                                                                                                               tttcacccacttagcaacgacaGCGCACGAGCTTATGACATAAAGCGGTCGATACCGAAGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGCGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGCCGATATCCa                                                                 >  1:368375‑M2/23‑205 (MQ=255)                                                                                                                                                                                                                                 tcacccgcttagcaacgacaCCCCACGAGCTTCTGACATAAAGCGGTCGATACCGAGTCTCAGCGCCAGACCCCCTACCGCCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTACCAGCGTGATTTGGTGGACGATATCCTTCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAgg                  <  2:175549‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                   acccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAccgcc                              >  2:396910‑M2/19‑236 (MQ=255)                                                                                                                                                                                                                                   acccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACAccgcc                              <  1:396910‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                                     ccacttagcaacgacaGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTatat             >  2:89107‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                                cgccAGCGCACGAACTTCGGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCCCCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCAACTACCCGTTACCCCGGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTGGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCAt  <  1:218269/248‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      GAAGAAGTGAAATAAACCTGAACTCACCCAGGTGGTAAGTTCAGGCACAGGGGTCAATTATGCGCAAACACCCGCACTCGGGGAAGGGAGTGCGGGCATAAGTGATGAGATTAAGAGGATAATTCGTGCGTTTTGCCATCCGGCGCACGATTATTCAGCACATCGTCCAGTTTTTTGTGGTCCAGTTCTTTCACCCACTTAGCAACGACAATGGTCGCTACGCCGTTACCGACCAGGTTAGTCAGCGCACGAGCTTCTGACATAAAGCGGTCGATACCGAGGATCAGCGCCAGACCCGCTACCGGCAAATGGCCCACCGCAGAGAGCGTCGCCGCCAGCACGATAAAGCCACTACCCGTTACCCCTGCCGCCCCTTTAGAAGAAAGCAGCAACACGATTAACAGCGTGATTTGGTGGACGATATCCATCTGACTGTTAGTGGCCTGGGCGATAAACACCGCCGCCATTGTCAGGTATATCGATGTGCCAT  >  NC_000913/3682051‑3682540

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