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	53 [0]	[2] 56	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)	53 (1.050)	49/476	0.1	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

ATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258435‑258933                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           cctcatcatctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                                                                                                                                                                                                                                           <  2:192141‑M2/2‑1 (MQ=255) cctcatcactttctgacggtaaattattatacatcccatttaacccggtgtatcggaatttcgacaaaggttttcaggaaatccccgactggatcgaaaaacgcttgtgttaataatttgctaaattttgccaatgttgtaaaacagttgcatcacaacttgcgatagcaatgacttcaccgagtggacccccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                                                                                                                                                                                                                                           <  1:488721‑M2/2‑1 (MQ=255)            gctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATaa                                                                                                                                                                                                                                               <  1:336473‑M2/13‑1 (MQ=255)                acggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGa                                                                                                                                                                                                                                           <  1:38845‑M2/17‑1 (MQ=255)                   gttaattagtagatatcccagttaacccgttgtatctgaattttgtcaaagttcttcaggaaatcactgactggaacgaaaaacgcttttgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggatatagcaatgacgttacctagtggactcctgaagtgcagattgatcaaaaaatttaccgcactaggcccggatattcgtgaaggtaATGGCAAAG‑ATAATCGATTc                                                                                                                                                                                                                                        <  2:70813‑M2/20‑1 (MQ=255)                             agagatcccatttaacccggtgtatcggacttttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacacgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTttcgctt                                                                                                                                                                                                                              >  2:162801‑M2/222‑248 (MQ=255)                               agatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGt                                                                                                                                                                                                                            <  1:220634‑M2/32‑1 (MQ=255)                                gatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgttcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGCTTCTTTTTCGAATGTc                                                                                                                                                                                                                           >  1:7312‑M2/219‑251 (MQ=255)                                   cccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTgg                                                                                                                                                                                                                        >  1:347521‑M2/216‑251 (MQ=255)                                           accgggggtatcggaattttgacaaaggtcttcaggaaatcaccgcccggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgcgtggacacccgaagagcagattgatcaaaaactttaccgaactaggcccttatattcgtgaaggaaAGTGCAAAGCAT‑ATAGATTCTTTTTCGTATGCCTGGCGGTAtgg                                                                                                                                                                                                                >  2:543719‑M2/208‑250 (MQ=255)                                             ccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCtttt                                                                                                                                                                                                                                    >  1:96951‑M2/206‑229 (MQ=255)                                             ccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCtttt                                                                                                                                                                                                                                    <  2:96951‑M2/24‑1 (MQ=255)                                                         aattttgtcaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                                                                                                                                                                                                  >  2:392230‑M2/194‑251 (MQ=255)                                                                  aaaggtcttcaggaaataaccgactggatcgaaaaccgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaaaggcagattgatcaaaaaatttacggcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTATGGCTGTATGCGTCAACGTGAAACCGGCACCgg                                                                                                                                                                                          >  2:270368‑M2/185‑250 (MQ=255)                                                                    aggtattcaggaagtgaacgactggatcgaaaaacgattgtgttaaaaatttgctaaatattgccaatttggtaaaacagtagcatcacaaaaggagatagcaatgacgttaccgagtggacacccgaagagaagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggaaAGGGAAAAG‑ATAATCGATTCTTTTTGGATTGTCTGGCTGTATGAGTCAACGagaaaccgggaacgga                                                                                                                                                                                         >  2:529505‑M2/183‑242 (MQ=255)                                                                                      cgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCt                                                                                                                                                                     <  2:425105‑M2/87‑1 (MQ=255)                                                                                      cgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcataacaacaggagatagcaatgacgttaccgagtggacacccgaagagaagattgataaaaaaatttaccgcactaggcccgtatattcgggaaggtaAGTGAAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGAAGGCGTCAAGGTGAAACCGCCACCGGAAGTGAGAGAGttctggtgcg                                                                                                                                                                     >  2:464053‑M2/165‑247 (MQ=255)                                                                                             aactaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcccaacaggagatagcaatgacgttaccgagtggacacccgaagagcatattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGAGTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTTAAGTCTGGGGCTGgtggag                                                                                                                                                              <  2:47735‑M2/94‑2 (MQ=255)                                                                                               cgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggatgg                                                                                                                                                            >  1:384705‑M2/156‑251 (MQ=255)                                                                                                  aaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGc                                                                                                                                                         <  1:139791‑M2/99‑1 (MQ=255)                                                                                                    aacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgctcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTg                                                                                                                                                       >  1:351895‑M2/151‑251 (MQ=255)                                                                                                          tgtgttaaaaaattgctaaattttgccaatttggtaaaacagtttcatcacaacaggagatagcaatgacgtgtccgagtggacacccgaagagcagattgatcaaaaaattttccgcattaggcccttatagtcgtgaaggtaAGTGCAAAG‑ATAATCGCTTCTTTTTCGATTGTCTGGCTTTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGTAagg                                                                                                                                                  <  2:343953‑M2/106‑2 (MQ=255)                                                                                                           gtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcgc                                                                                                                                                >  1:84248‑M2/144‑251 (MQ=255)                                                                                                           gtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcg                                                                                                                                                 <  1:290757‑M2/107‑1 (MQ=255)                                                                                                               taaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAgg                                                                                                                                             >  1:470604‑M2/140‑250 (MQ=255)                                                                                                                    atttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATccc                                                                                                                                       >  1:432826‑M2/135‑251 (MQ=255)                                                                                                                              ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCt                                                                                                                              >  1:386038‑M2/125‑250 (MQ=255)                                                                                                                                ttgccaatttggtaaaacagttgcatcacaacaggagctagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                            >  1:197539‑M2/123‑250 (MQ=255)                                                                                                                                  gccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggaaacacgaagagaagattgataaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCCAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGGGAATTCTGAGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGt                                                                                                                         >  2:484431‑M2/121‑251 (MQ=255)                                                                                                                                   ccaatttggtaacacagttgcatcacaacaggagatagcaatgacgttaccgggtggacacccgaagagcagattgatccaaacatttaccgcactaggcccgtatcttcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGAGGTTGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGtt                                                                                                                        >  1:50658‑M2/120‑251 (MQ=255)                                                                                                                                      atttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGAGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGTAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAcc                                                                                                                     >  1:369897‑M2/117‑251 (MQ=255)                                                                                                                                                     tgcatcgcaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgaacaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTc                                                                                                       <  2:200497‑M2/149‑1 (MQ=255)                                                                                                                                                       catcacaacaggagatagcaatgacgttaccgagtggacaaccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaatgtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGAt                                                                                                    >  1:465080‑M2/100‑251 (MQ=255)                                                                                                                                                            caacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAACCCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGGGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                               >  2:70142‑M2/95‑251 (MQ=255)                                                                                                                                                             aacaggagatagcaataacttgaccgaggttacacctgaagagcagattgatcaaaaaagttaccgcactagtcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTAGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGTGTTGGTAGCTCAGGAATCCCGTTTTACCTACAGTTCCCAGTTTGGTCTGGTCGATAAAGCa                                                                                              <  2:298808‑M2/158‑1 (MQ=255)                                                                                                                                                              acaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGGTACCAGTTTGGTCTGTTCGATAAAGc                                                                                               >  2:153081‑M2/93‑249 (MQ=255)                                                                                                                                                                   agatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGAc                                                                                        >  1:478184‑M2/88‑251 (MQ=255)                                                                                                                                                                        gcaatgacgttaccgagtggacacccgaagagccgattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTAAACGTGAAACCGGCACCGGAAGTGCTTGAATTCTGGGGCTGGTGGATGGGGCTGGAAGCGACGGAATCCCGTTTTACCTACAGTTACCAGTTTGGGCTGTTCGATCAAGCAGGCGACtggc                                                                                    >  2:388220‑M2/83‑249 (MQ=255)                                                                                                                                                                          aaagacgtgaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgggaaggtaAGTGCAAAG‑ATAAGCGATTCTTTTTCGATTGTCTGGCGGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGAGGGAGCTGGAAGCGCAGGAATCCCGTTTTTCCTACAGTTACCGGTTTGGGCGGTTCGATAAATCAGGCGACTGGAaga                                                                                 <  2:275330‑M2/171‑1 (MQ=255)                                                                                                                                                                           atgacgttaccgagtggacacacgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTATGTTCGATAAAGCAGGCGACTGGAagag                                                                                <  1:162801‑M2/172‑1 (MQ=255)                                                                                                                                                                             gacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAgtg                                                                              <  1:179445‑M2/174‑1 (MQ=255)                                                                                                                                                                                            acacccgaagagcagattgaacaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTTAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGaca                                                               <  1:543719‑M2/189‑1 (MQ=255)                                                                                                                                                                                             cacccgaagagcagattgatcaaaaaatttaccgcactaggccagtatattcgtgaaggtaAGTGCAAAG‑CTAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGACTTCTGGGGCTGGGGGATGGAGCTGGTAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGATGGGTGTTCGGGTAAAAGgcac                                                              >  2:436567‑M2/62‑251 (MQ=255)                                                                                                                                                                                                  gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTg                                                            <  2:192959‑M2/192‑1 (MQ=255)                                                                                                                                                                                                  gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTg                                                            >  1:192959‑M2/57‑248 (MQ=255)                                                                                                                                                                                                    agagcagattgatcaaaaaatttaccgcactagacccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTg                                                       <  1:368288‑M2/197‑1 (MQ=255)                                                                                                                                                                                                                     aaatttaccgcactaggcccgtatatgcgtgaaggtaAGTGCAAAG‑ATAAGCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGca                                      <  1:407625‑M2/214‑1 (MQ=255)                                                                                                                                                                                                                        tttaccgcactaggcccgtatattcgtgcaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCacac                                   >  1:336178‑M2/35‑251 (MQ=255)                                                                                                                                                                                                                           accgcactaggcccgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCt                                <  1:23566‑M2/220‑1 (MQ=255)                                                                                                                                                                                                                                        cgtatattcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCa                    >  1:51506‑M2/19‑250 (MQ=255)                                                                                                                                                                                                                                          tatattagtgaaggtaAGAGCAAAG‑ATATTCGATTCTTTTTCTATTATCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAGGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGCTAAAGCAGGCGACTGGAAGCGTGTTCCGGTAAAAGACACTGAAGTGGTTGACCGCCCGGAGCACACCCTGAGTACGTTTCAcgg                 >  2:378024‑M2/17‑250 (MQ=255)                                                                                                                                                                                                                                             attcgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTATGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACgaga               <  1:270368‑M2/237‑1 (MQ=255)                                                                                                                                                                                                                                                cgtgaaggtaAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCt           <  1:57016‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                                    aaGGTAAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGt       >  1:451271/1‑251 (MQ=255)                                                                                                                                                                                                                                                      ggTAAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGg     <  1:252852‑M2/251‑2 (MQ=255)                                                                                                                                                                                                                                                          aGTGCAAAG‑ATAATAGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGCAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAACAGACCATGAAGTGGTTGAACGAATGGAGAACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGggtg  <  2:503780‑M2/250‑5 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           ATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAG‑ATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258435‑258933

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