Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258676	258683	8–776	37 [0]	[1] 38	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	37 (0.680)	35/472	0.3	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGT  >  NC_000913/258438‑258928                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     catcatctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagagtgattaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                                                                                                                                                                                                   <  1:446976‑M2/4‑1 (MQ=255)         gctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATa                                                                                                                                                                                                                                           <  1:373040‑M2/12‑1 (MQ=255)         gctgacggtaaattattagagaacccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccttatattcgtgaaggtaAGTGCAAAGAt                                                                                                                                                                                                                                            <  1:363612‑M2/11‑1 (MQ=255)           tgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                                                                                                                                                                                                                        >  1:471903‑M2/237‑251 (MQ=255)                gtaaattattagagatcccatttaacccggtgtatcggcattttgacaaaggtattcaggcaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggaaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagcttgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                                                                                                                                                                                                                                   >  1:451513‑M2/232‑251 (MQ=255)                gtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaacaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcaaaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcacaaaatttaccgcactaggcccgtatattcgtgaagggaAGTGCAAAGCTAATCGattt                                                                                                                                                                                                                                   >  2:383951‑M2/232‑250 (MQ=255)                  aaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCt                                                                                                                                                                                                                                  <  2:413759‑M2/21‑1 (MQ=255)                            agatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                                                                                                                                                                                                       <  1:302837‑M2/32‑1 (MQ=255)                             gatcccatttaacacggtgtaacgggattttgacaatgttcttccggaaatcaccgactggatcgaaaaacgcttgtgttcaaaatttgctaaattttgccaatttggtaaaacagttgcatcacactaggagagagcaatgacgttaccgagcggacacccgcagagcagattgataaaaaaatttaccgcactaggcccgtatattcgtgaaggtaCGTGAAAAGATACACGGTTCTTTTTCGATTGTc                                                                                                                                                                                                                      >  2:332745‑M2/219‑251 (MQ=255)                                  cattgaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtggtaaaaagttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgaggggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatatgcttgaaagtaAGTGCAAAGATAATCGTGTCGGTTGCGATTGTCTGGCt                                                                                                                                                                                                                 <  2:451513‑M2/38‑1 (MQ=255)                                            ggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaaagcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGCTAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTc                                                                                                                                                                                                        >  2:48883‑M2/204‑250 (MQ=255)                                                 atcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCTATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTg                                                                                                                                                                                                  <  1:242601‑M2/53‑1 (MQ=255)                                                   cggaattttgaccaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGa                                                                                                                                                                                                 >  1:147552‑M2/197‑250 (MQ=255)                                                     gaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc                                                                                                                                                                                              <  1:342506‑M2/57‑1 (MQ=255)                                                        ttttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCgg                                                                                                                                                                                           >  1:82431‑M2/192‑251 (MQ=255)                                                                                                aaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCt                                                                                                                                                   >  2:126912‑M2/152‑251 (MQ=255)                                                                                                                aatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATcc                                                                                                                                   >  2:298124‑M2/136‑251 (MQ=255)                                                                                                                      ctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                                                                                                                             >  1:321963‑M2/130‑251 (MQ=255)                                                                                                                      ctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                                                                                                                             >  2:323801‑M2/130‑251 (MQ=255)                                                                                                                                  aatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGAGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc                                                                                                                 >  1:348714‑M2/118‑251 (MQ=255)                                                                                                                                      tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGt                                                                                                             >  1:162002‑M2/114‑251 (MQ=255)                                                                                                                                       ggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGtt                                                                                                            <  2:430908‑M2/139‑1 (MQ=255)                                                                                                                                           aaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGGGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGt                                                                                                        <  2:452525‑M2/143‑1 (MQ=255)                                                                                                                                                    catcacaacaggagatagcaatgacgttaccgagtggacacccgaagaacagattgatcaaaaaatttaccgcactagtcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTGTTTCGATTGTCTGGCTGTATGCGTCATCGTGAAACCGGCTCCGGAAGTGCGTGAAGTCTGGGGCTGGTGGATGGAGCTGGAAGCGTAGGAATCCCGTTTTACCCACAGTTTTCGTTTTCGTCTGTTCGAt                                                                                               <  2:296239‑M2/152‑1 (MQ=255)                                                                                                                                                                       aatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGGGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGCCTGGa                                                                               >  2:315034‑M2/81‑248 (MQ=255)                                                                                                                                                                          gacgttaccgagtggactcccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAgt                                                                          <  1:257832‑M2/173‑1 (MQ=255)                                                                                                                                                                                 ccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGTAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCGGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGt                                                                  <  2:471903‑M2/181‑1 (MQ=255)                                                                                                                                                                                              cgaagagcagattgatcaaaaaatttaccgcactaggcccgtacattagtgaaggtaAGTGAAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAACCCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGGGTTCAGGTAAAAGACACTGa                                                      >  1:44074‑M2/58‑250 (MQ=255)                                                                                                                                                                                               gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGa                                                      >  2:363170‑M2/57‑249 (MQ=255)                                                                                                                                                                                                 agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTg                                                  >  1:22526‑M2/55‑251 (MQ=255)                                                                                                                                                                                                        attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAAc                                           <  2:168992‑M2/204‑1 (MQ=255)                                                                                                                                                                                                              caaaaaatttaccgcacgaggcccgtatattcggtaaggtaAGTGCAAAGATAATCGATGCGTTTTCGAGTGTCTGGCGGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGGTGGTGGATGGAGCTGGAAGGGCAGGAATCCCGTTTTACCTACAGTTTCCAGTTTGGTCGTTTCGATAAAGGAGGCTACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTgg                                     <  1:280224‑M2/210‑1 (MQ=255)                                                                                                                                                                                                                aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTgg                                     >  2:193741‑M2/40‑249 (MQ=255)                                                                                                                                                                                                                 aaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGc                                  <  1:107710‑M2/213‑1 (MQ=255)                                                                                                                                                                                                                        accgcactaggcccgggtattcgtgaaggtaAGTGCAAAGATAATCGCTTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGCAGTGCGTGAATTCTGGGGCTGGTGGATGTAGATGGATGAGCAGGAATCCCGTTTTAACTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTACGGTAAAAGACACTGCAGGGgttgaacgactgaggaccaccct                           >  2:268733‑M2/32‑240 (MQ=255)                                                                                                                                                                                                                                  gcccgtatattcgtgaaggtaAGTGCAAAGATAAGCGATTCTTTTGCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGGGAGttt                 <  2:48526‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                     cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGGTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACGCTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCAc              <  2:50546‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                                 aaGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGt  <  2:82431/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGT  >  NC_000913/258438‑258928

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