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	59 [0]	[3] 62	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)	59 (0.820)	56/472	0.2	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

CATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257665‑258141                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           catcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAACCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                     >  1:26305‑M1/1‑243 (MQ=255)         gACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                                                                                                                                                                                              <  2:178407‑M1/249‑15 (MQ=255)                             cccATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                                                                                                        >  1:113003‑M1/1‑215 (MQ=255)                             cacatTGAACCCGGTGTATCGGAATTTTGACAAAGGGCTTCAAGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGGGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                                                                                                        <  2:760‑M1/249‑37 (MQ=255)                                  gtaACCCGGTGTATCGGAATTTTGACAAAGGTCGTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGGAagtgcaaagataatcgattctttttcgattgtctggctgt                                                                                                                                                                                                    <  2:168261‑M1/249‑41 (MQ=255)                                    aaCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatg                                                                                                                                                                                                 <  2:59204‑M1/251‑44 (MQ=255)                                                   aaTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc                                                                                                                                                                                  >  1:420043‑M1/1‑193 (MQ=255)                                                       ttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACCCTTGTATTACAAATTTGCTAAATTTTGCCCATTTGGTAAAAAAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcc                                                                                                                                                                              >  1:319301‑M1/1‑189 (MQ=255)                                                         gACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacc                                                                                                                                                                            >  2:453293‑M1/1‑187 (MQ=255)                                                                gTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtg                                                                                                                                                                     <  2:263380‑M1/251‑72 (MQ=255)                                                                    tCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTTGTAAAACAGTTGCATCACAACAGGAGATAGCAATGGCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGGATATTCGTGAAGGTAagggcaaagataatcgattcttggtcgagtgtctggctgtatgcgtcaacgtgagtccggcactggaagtgcgt                                                                                                                                                                  <  2:598971‑M1/250‑75 (MQ=255)                                                                    tCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtg                                                                                                                                                                 >  1:379552‑M1/1‑176 (MQ=255)                                                                       ggAAATCACCGACTGGATCGACAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCAAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACTGAGTGGAACCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcagcgtgaaaacggcaccggaagtgcgtgaat                                                                                                                                                              >  2:343254‑M1/1‑173 (MQ=255)                                                                       ggAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaat                                                                                                                                                              >  1:571937‑M1/1‑173 (MQ=255)                                                                                       aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                                                                                                                              <  2:652441‑M1/251‑95 (MQ=255)                                                                                          gAAAAACGCTTATGTTAAAAATTTGCTAAATTTTCCCCATTTGGTAAAACAGTCGCATCACAACAGGAGATAGCAATGACGATACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATCTTCGAGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaatggcgtgaattctggggctggtggatgg                                                                                                                                            >  2:291683‑M1/1‑154 (MQ=255)                                                                                             aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAACACAGTTGCATAAAAAAAGGAGATAACAATGACGTTACCGAGGGGACACCAGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATCCGTGAGGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggacgtgcgtgaattctggggctggcggatggagcg                                                                                                                                        >  2:324466‑M1/1‑151 (MQ=255)                                                                                                  cTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgga                                                                                                                                     >  2:525217‑M1/1‑146 (MQ=255)                                                                                                           aaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATAACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGAAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaa                                                                                                                           >  2:446599‑M1/1‑137 (MQ=255)                                                                                                               tttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAAACGGAGATAGCACAAACGTTACCGAGTGGACACCAGAAGAGCAGATTGATCAACAAATTTACCGCACTAGGCCCGTATATTCGTGAAGATAagtgcaaagataatcgattctttttcgattgtctggctgtatgagtcaacgggaaaccggcaccggaagtgcgtgaattctggggctggtgtatggagctggccgcgcaggaaacacg                                                                                                                      >  2:586916‑M1/1‑133 (MQ=255)                                                                                                                     aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                  <  2:376778‑M1/249‑123 (MQ=255)                                                                                                                      aaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggctggtggatggagctggaagcgcaggaatcccgttttacc                                                                                                               >  2:135921‑M1/1‑126 (MQ=255)                                                                                                                            gCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacag                                                                                                          <  1:453293‑M1/250‑131 (MQ=255)                                                                                                                             ccAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatccctttatacctacagtt                                                                                                        <  2:66362‑M1/251‑133 (MQ=255)                                                                                                                                aTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttacc                                                                                                     <  2:113003‑M1/251‑136 (MQ=255)                                                                                                                                   tGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                  >  1:454378‑M1/1‑113 (MQ=255)                                                                                                                                     gTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttt                                                                                                >  1:226474‑M1/1‑111 (MQ=255)                                                                                                                                       aaaaCAGTTGCATCACACCAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGCGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgg                                                                                              >  2:673728‑M1/1‑109 (MQ=255)                                                                                                                                           cAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctg                                                                                          >  2:177947‑M1/1‑105 (MQ=255)                                                                                                                                            aGTTGCATCACAACAGGAGCTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgt                                                                                         >  2:75954‑M1/1‑104 (MQ=255)                                                                                                                                            aGTTGCATCACAACAGGAGATAGCAATGACGTTACCGCGTGGACACCCGAAGCGCAGATTGCTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgt                                                                                         >  1:638440‑M1/1‑104 (MQ=255)                                                                                                                                                 cATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgat                                                                                    >  2:230650‑M1/1‑99 (MQ=255)                                                                                                                                                     acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatccagttttacctacagttaccagtttggtctgttagataa                                                                                  >  1:108225‑M1/1‑95 (MQ=255)                                                                                                                                                        acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcatgaatcccgttttacctacagttaccagtttggtctgttcgataaagcag                                                                             <  1:673728‑M1/251‑160 (MQ=255)                                                                                                                                                        acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                               >  2:531383‑M1/1‑92 (MQ=255)                                                                                                                                                             agaTAGCAATGACGATACCGGGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGGAagtgcaaagatactcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggcggctggagctggaagcgcaggaatcccgttttacctgcagttaccagcttggtctgttcgataaagcagtcgac                                                                        >  2:62963‑M1/1‑87 (MQ=255)                                                                                                                                                              gaTAGCAATGACGGCACCGAGTGGACACCCGAAGCGCAGATTGATCAAAAAATTTACCGCACTAGGCCAGCATATTCGTGAAGGTAagtgcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcggaggaatcccgttttacctacagttaccagtttggtctgtccgataaagcaggcgatt                                                                       >  2:642475‑M1/1‑86 (MQ=255)                                                                                                                                                                         cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgtt                                                            <  1:177947‑M1/251‑177 (MQ=255)                                                                                                                                                                          gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttc                                                           >  1:177783‑M1/1‑74 (MQ=255)                                                                                                                                                                            tACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacctgcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttcc                                                          <  1:291683‑M1/250‑179 (MQ=255)                                                                                                                                                                             aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCATGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccgg                                                        <  1:518687‑M1/251‑181 (MQ=255)                                                                                                                                                                                       cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtcggttcgataaagcaggcgactggaagagtgttccggtaaaagacac                                              <  2:453123‑M1/251‑191 (MQ=255)                                                                                                                                                                                        acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacact                                             >  1:593597‑M1/1‑60 (MQ=255)                                                                                                                                                                                         cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactg                                            >  2:564560‑M1/1‑59 (MQ=255)                                                                                                                                                                                          ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactg                                            >  1:245182‑M1/1‑58 (MQ=255)                                                                                                                                                                                          ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGAAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcaggaatcccggtttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtaatcgcggtaagctgttcct                                           >  1:248374‑M1/1‑58 (MQ=255)                                                                                                                                                                                            gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                           >  1:652572‑M1/1‑56 (MQ=255)                                                                                                                                                                                             aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgtgcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagt                                        <  1:164091‑M1/251‑197 (MQ=255)                                                                                                                                                                                              agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgcttgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcacgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtg                                       >  1:657893‑M1/1‑54 (MQ=255)                                                                                                                                                                                                   aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttg                                   >  2:513512‑M1/1‑49 (MQ=255)                                                                                                                                                                                                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaag                                                  >  1:434558‑M1/1‑48 (MQ=255)                                                                                                                                                                                                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaagcccgttttaccgacagttaccagtttggtctgttcgataaagcaggctactggaagagtgttccggtaaaag                                                  <  2:434558‑M1/234‑187 (MQ=255)                                                                                                                                                                                                     aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccggtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaac                                >  2:210878‑M1/1‑47 (MQ=255)                                                                                                                                                                                                         aTCAAAAAACTTACCGCACTAGGCCCGTATAGTCGTGAATGTAagtgcaaagataatcgattcttattcgattgtctggatgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaagtcgggggctggtggatggagctggaagcgcagcactcccgttgtacctacagttgccagtttggtctgtgcgataaagcaggcgactggaagagggatccggtaaaagacactgaagttgtttaacgact                            <  2:56897‑M1/251‑209 (MQ=255)                                                                                                                                                                                                             aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctg                                                                                                                                       >  1:97380‑M1/1‑39 (MQ=255)                                                                                                                                                                                                             aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctg                                                                                                                                       <  2:97380‑M1/140‑102 (MQ=255)                                                                                                                                                                                                              aaaaTTTACCGCACTAGGCCCGTATATTCGTGTAGGTAagtgcaatgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaggtgcgtgaattctggggctggtggatggagctggcaacgcaggcatcccgttttacctacagttaccagtttggtctgttcgataaagcagtcgactggaagagtgttccgttaaaagacactgaagaggttgaacgactggagc                       >  2:29599‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                         cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggcacaagacactgaagtggttgaacgactggagcacatcctgcgt            >  2:566592‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                    tataTTCGTGAAGGTAcgtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacg  >  2:573103‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCAC  >  NC_000913/257665‑258141

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: ‑