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	70 [0]	[0] 73	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)	70 (0.840)	57/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

CTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCA  >  NC_000913/257660‑258143                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             cTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                           >  1:71771‑M1/1‑248 (MQ=255)           gctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATATACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                                                                                                                                                                 <  2:517291‑M1/251‑14 (MQ=255)           gctgACGGTAAATTATTAGAGAGCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAATTCACCGACTGGATCTAAATACGCTTGTTTTAAAAATTTGCTAAATTTTGCCAATTTGGTGAAACTGTTGCATCACAACAGGAGTTAGCAAGGACGTTACCGAGTGGGC‑ACACGAAGAGGAGATTGATCAAAAAATTTACCTCTATAGGCACGTATATTCGTGAAGGTAagtgcaaagataa                                                                                                                                                                                                                                 <  1:701767‑M1/251‑14 (MQ=255)               aCGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGG‑ATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAAGTTTGCCAATTTTGTAAAACAGTTGGATCACAACAGGAGATAGCAAGGACGTTACCGAGTGGAC‑ACGCGAAGAGTAGATTGATCAAAAAAGTTACCTCACTAGGGCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                                                                                                                                                            <  1:64920‑M1/251‑19 (MQ=255)                cGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                                                                                                                                                            <  1:564516‑M1/251‑19 (MQ=255)                 ggTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGACTTTTGACAAAGGTCTTCAGGAACTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGCGCAGATTGATCAAACAATTTACAGCACTAGGCCCGTATATTCGTGAAGGCAagtgcaaagataatcgatt                                                                                                                                                                                                                           >  1:721506‑M1/1‑232 (MQ=255)                  gTAAATTATTAGAGATCCCATGTAACCCGGTTTATCTGAATTTTGACTAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATGTGGTAAAACAGTTGCATCACAACAGGAGAAAGCAATGACGTTTCCGAGTGGTC‑GCTCGATGAGCATATTGATCAAAAAATTTACCTCACTAGGCCCGTATATTCGTAAAGGTAagtgcaaagataatcgattc                                                                                                                                                                                                                          <  1:668276‑M1/251‑21 (MQ=255)                            agagATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                                                                                                                                                                <  1:634394‑M1/251‑31 (MQ=255)                               gaTCCAATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACAGACTGGATCGCACAACGCTTGTGTTAAAAATTTGATCAATTTTGACAATTTGGTAAAACAGTTGCATCACACCAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCGAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                                                                                                                                                             >  2:318788‑M1/1‑218 (MQ=255)                                   ccATTTAACCCGGTGTATCGGCATTTTGACAACGGTCTTAAGGAAATCACCGACTGGTTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCACTGACGTTACCGAGTGGAC‑ACCCGCAGAGCAGATTGATCAAAAAATTTACCGCCCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc                                                                                                                                                                                                         >  2:707056‑M1/1‑214 (MQ=255)                                        tAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                                                                                                                    >  2:251399‑M1/1‑209 (MQ=255)                                        tAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                                                                                                                    >  1:515072‑M1/1‑209 (MQ=255)                                           cccagtgtagcgGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATGGGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑CACCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGCGAAGGTAagtgcaaagataatcgattctttttcgagtgtctggctgtatgcg                                                                                                                                                                                                 <  2:239750‑M1/247‑46 (MQ=255)                                                           tttGACAAAAGTCATCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAAATTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                                                                                                                                                                 >  1:178950‑M1/1‑190 (MQ=255)                                                                         tcaggtcttaaacgactggatagaaAACCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAAAAGTTGAATCAACACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGACGAGCAGATTGGTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagagcgtg                                                                                                                                                                   >  2:645982‑M1/13‑176 (MQ=255)                                                                                aTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctg                                                                                                                                                            >  1:333815‑M1/1‑169 (MQ=255)                                                                                       cctGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagggcgtgacttctggggctgg                                                                                                                                                     >  2:28729‑M1/2‑162 (MQ=255)                                                                                       aCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgg                                                                                                                                                     >  1:225788‑M1/1‑162 (MQ=255)                                                                                        cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggg                                                                                                                                                    >  2:254421‑M1/1‑161 (MQ=255)                                                                                        cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                                                                                                                                                                          >  1:412115‑M1/1‑161 (MQ=255)                                                                                        cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                                                                                                                                                                          <  2:412115‑M1/181‑21 (MQ=255)                                                                                          ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgg                                                                                                                                                  >  2:509138‑M1/1‑159 (MQ=255)                                                                                                  aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagct                                                                                                                                          >  1:684550‑M1/1‑151 (MQ=255)                                                                                                          gtgtTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcg                                                                                                                                   <  1:691758‑M1/250‑108 (MQ=255)                                                                                                          gtgtTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcg                                                                                                                                   <  1:251399‑M1/250‑108 (MQ=255)                                                                                                              tAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGCAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggagggagctggaagcgcagg                                                                                                                               >  2:591269‑M1/1‑139 (MQ=255)                                                                                                                            aTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacct                                                                                                                >  2:136559‑M1/1‑125 (MQ=255)                                                                                                                              tttGCCAATTTGGTCAAACAGTTGCATCACAACAGGCGATAGCAATGACGTTACAGAGTGGAA‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggcggatggagctggaagcgcaggaatcccgttttccctac                                                                                                              >  2:270274‑M1/1‑123 (MQ=255)                                                                                                                                tGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggcatcccgttttacctac                                                                                                              >  1:443377‑M1/1‑121 (MQ=255)                                                                                                                                    aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttac                                                                                                        >  1:518656‑M1/1‑117 (MQ=255)                                                                                                                                       ttGGTAAAACAGTTGCATCACAACAGGAGCTAGAAATAACGTTACCGAGAGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATCTTCGTGAAGGTCagtgcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaca                                                                                                       >  2:497863‑M1/1‑114 (MQ=255)                                                                                                                                          gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttt                                                                                                  >  2:324695‑M1/1‑111 (MQ=255)                                                                                                                                          gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttt                                                                                                  >  2:137195‑M1/1‑111 (MQ=255)                                                                                                                                           tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcttgaattctgggg                                                                                                                                                         <  1:254717‑M1/195‑86 (MQ=255)                                                                                                                                           tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggg                                                                                                                                                         >  2:254717‑M1/1‑110 (MQ=255)                                                                                                                                                  gTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAACATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgcttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagagcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtt                                                                                          >  1:628927‑M1/1‑103 (MQ=255)                                                                                                                                                         cacaACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataa                                                                                    >  2:332577‑M1/1‑96 (MQ=255)                                                                                                                                                           caacaGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                 >  1:162526‑M1/1‑94 (MQ=255)                                                                                                                                                            aacaGGAGATAGCACTGCCGTTACCGAGTGGCC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGCAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaccgttaccagtttggtctgttcgataaagc                                                                                 >  2:671103‑M1/1‑93 (MQ=255)                                                                                                                                                              caGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcagtaatcccgttttacctacagttaccagtttggtctgtgcgataaagcagg                                                                              <  2:683749‑M1/251‑161 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGCGTGGAC‑AACCGAAGCGAAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGAGCAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgccggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggtgact                                                                         >  2:464401‑M1/1‑86 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccag                                                                                                     <  2:568930‑M1/223‑138 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccag                                                                                                     >  1:568930‑M1/1‑86 (MQ=255)                                                                                                                                                                   gaTAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtgttacctacagttaccagtttggtctgttcgataaagcaggcgact                                                                         <  2:684550‑M1/251‑166 (MQ=255)                                                                                                                                                                      aGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactgga                                                                      <  1:509138‑M1/251‑169 (MQ=255)                                                                                                                                                                      aGCAATGACGTGACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAAGTTACCTCCCTAGGCCCGTATAGTCGTGCAGGTAagtgcaaagataatcgagtctttttctattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagttcgtgaattctgtggctggtggaaggagctggaagcgcatgaatcccgttttacccacagttaccagtttggtcgtttcgataaagcaggcgactgg                                                                       <  2:326454‑M1/250‑168 (MQ=255)                                                                                                                                                                             aCGTTACCGAGTGGACGACCC‑AAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgttaattctgggtctggtggatggagctggaagcgcaggaatcccgttgtacctacagtgaccagtttggtctgttcgataaagcaggccactggaagagtgt                                                               <  2:128650‑M1/251‑176 (MQ=255)                                                                                                                                                                             aCGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGAGCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctgtttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccggtttgcctacagtgatcagtttggtctgttcgataaagcagtcgactggaagagtgt                                                               <  2:65140‑M1/251‑176 (MQ=255)                                                                                                                                                                              cGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaagccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgtt                                                              <  2:443377‑M1/251‑177 (MQ=255)                                                                                                                                                                                 tACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtgcgataaagcaggcgactggaagagtgttcc                                                            <  2:582030‑M1/250‑179 (MQ=255)                                                                                                                                                                                     gAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaa                                                       <  1:671103‑M1/251‑184 (MQ=255)                                                                                                                                                                                      aGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggagaagtgtccggttaa                                                       >  1:716312‑M1/1‑67 (MQ=255)                                                                                                                                                                                          gAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcagtaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagac                                                  <  1:497863‑M1/251‑189 (MQ=255)                                                                                                                                                                                              acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggtgctggtggatgtagctggaagctcaggaatcccgttttaccttcagttaccagttgggtcgggtcgataaagcaggcgactggaagagtgttccggtaaaagacact                                               <  1:464401‑M1/251‑192 (MQ=255)                                                                                                                                                                                              acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctgcagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacact                                               >  1:149896‑M1/1‑60 (MQ=255)                                                                                                                                                                                               cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagccactg                                              >  2:209055‑M1/1‑59 (MQ=255)                                                                                                                                                                                                         aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttg                                     >  1:567039‑M1/1‑49 (MQ=255)                                                                                                                                                                                                         aGATTGATCAAAAAATTTACCGAACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttg                                     >  1:198121‑M1/1‑49 (MQ=255)                                                                                                                                                                                                              gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacttgaaaccggcaccggaagtgcgtgaattctggggctggtggagtgagctggaagcgcaggaatccctttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgac                               <  1:707056‑M1/251‑208 (MQ=255)                                                                                                                                                                                                                 cAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactg                             <  2:198121‑M1/250‑210 (MQ=255)                                                                                                                                                                                                                   aaaaaTTTACCGCACTAGGCACGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgacgtggttgaacgactgg                            >  2:722666‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                             cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgc                <  1:318788‑M1/251‑223 (MQ=255)                                                                                                                                                                                                                              gCACTCGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagtgcgtgacttctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgatcaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcg               >  2:45605‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                              gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtgcgataaagcaggcgactggaagagtgttccggtaaaagacactgatttggttgaacgactggagcataccctgcg               <  1:254421‑M1/251‑224 (MQ=255)                                                                                                                                                                                                                               cACTAGGCCCGTATATTTGAGAAGGTAagtgcaaagatcctcgattctttttcgattgtctggctgtatgagtcaaagtgaaaccggaaccggaagtgcgtgaattatggggctggtggatggcgctggaagcgcaggaatccagttttacatacagttaccagtttggactgttcgataaagcaggcgactggaagcgtgttccggaacatgacactgaagtggttgaacgactggagcacaccctgcgt              >  2:650155‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                   aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgcagtggttgaacgactggagcacaccctgcgtgagt          >  2:240516‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                   aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagt          >  1:443234‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                     gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagttt        >  1:586638‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                        cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcac     <  1:28729‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                                           atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccggttgacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgtgccggtaaaagacactgaagtgggtgaacgactggagcacaccctgcgtgagtttcacgag  <  1:681114‑M1/251‑237 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             CTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGAC‑ACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCA  >  NC_000913/257660‑258143

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