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	87 [3]	[2] 87	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	3 (0.030)	84 (0.900)	73/478	0.2	98.3%	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

CCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGAC  >  NC_000913/257659‑258157                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           ccTCATCAGCTGCTGACGGGAAATTATTAGAGAGCCCATTTAACCCGGTGTATCGGAATGTTGACAAAGGTCATCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATATGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGTAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGAATGATCAAAAAAGTTACCGCTCTAGGCCCGTATATTCGTGAAGGTAag                                                                                                                                                                                                                                                          <  2:492600‑M1/251‑3 (MQ=255)    catcatCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                                         <  2:666479‑M1/249‑4 (MQ=255)       catcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                                     <  1:374021‑M1/250‑8 (MQ=255)         tcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGCATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                                     >  1:602794‑M1/1‑241 (MQ=255)            gctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTGCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAAGTTGCTAAAGTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTGACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagata                                                                                                                                                                                                                                                <  2:970368‑M1/250‑13 (MQ=255)                   gTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                                                                                                                                                                                        <  2:977116‑M1/251‑21 (MQ=255)                     aaaTTATTTGAGATCCCATTTAACCCGGTGTATTGTAATTTTGACTAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTTTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGATCAGATTGATCAAAAAATTTCCCGCACTAGGCCCGTATATTCGTGAATGTAagtgcaaagataatcgattctt                                                                                                                                                                                                                                      <  2:496750‑M1/251‑23 (MQ=255)                     aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                                                                                                                                                                                                                      <  2:369344‑M1/251‑23 (MQ=255)                       attattAGAGATCCCATTTACCCCGGTGTCTAGGACTTTTGACAAAGGACTTCAGGCAATCACCGACTGGATCGAAAAACGCTCGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATAACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGCGCAGATTGATCAACAAATTTACCGCACTAGGCCCGAATATCTGTGCAGGTAagagcaaagataatcgattctttt                                                                                                                                                                                                                                    >  2:68840‑M1/1‑227 (MQ=255)                          attaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcg                                                                                                                                                                                                                                 <  1:155059‑M1/251‑28 (MQ=255)                             agagAACCCAGTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATTACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAAGTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTAGATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                                                                                                                                                                              <  1:817237‑M1/251‑31 (MQ=255)                              gagaTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattg                                                                                                                                                                                                                             <  2:195767‑M1/251‑32 (MQ=255)                                   cccATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTCCCGAGTGGACACCCGAAGAGCAGATTGATCAAGAAATTTACCGCACTCGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                                                                                                                        >  1:958258‑M1/1‑215 (MQ=255)                                     cATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggct                                                                                                                                                                                                                      <  2:756219‑M1/251‑39 (MQ=255)                                           aCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATTGAAAAACGCTTGTGTTAAAAATTTGCTAACTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGCTCAAAAAATTTACCGCACTAGGCCAGTATATTCGTGCAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                                                                                                                                                                                >  1:907546‑M1/1‑207 (MQ=255)                                                          aTTTTGACACAGGTCTCCAGGAAATCACCGCCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGCTGCATCACAACAGGGGATAGCAATGCCGTTACCGAGGGGACACCCTAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaacgataatcgattctttttcgattgtctgtctgtatgcgtccactggaaac                                                                                                                                                                                                   >  2:255807‑M1/1‑192 (MQ=255)                                                                    aGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgga                                                                                                                                                                                         >  1:680215‑M1/1‑182 (MQ=255)                                                                       tCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTCAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGCTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgacaccggcaccggaagtgc                                                                                                                                                                                    >  1:1131‑M1/1‑179 (MQ=255)                                                                          tCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                                                                                                                                  <  1:667370‑M1/250‑75 (MQ=255)                                                                            aGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaa                                                                                                                                                                               <  2:978467‑M1/251‑78 (MQ=255)                                                                              gAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaatt                                                                                                                                                                             <  2:699841‑M1/251‑80 (MQ=255)                                                                                  tCACCGACTGGATCGACAAACGCTTGTGTTAAAAATTTGCTAAATTTTACCAATTTGGTAAAACAGTTGCCTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCCGATTGCTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGCAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacttgcaaccggcaccggaagtgcgtgatttctgg                                                                                                                                                                         >  2:939361‑M1/1‑168 (MQ=255)                                                                                   cACCGACTGGATTGAAAAACGCTTGGGTTAAAAAGTTGCTAAATTTTGCCAAGTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAGTACGGTCCCGAGTGGACACCCGAAGAGCAGATTTATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgagtcttttgcgattgtctggctgtatgcggcaacgtgaaaccggcaccggaagtgcgtgaatgctgg                                                                                                                                                                         <  2:433067‑M1/250‑84 (MQ=255)                                                                                      cGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                                                                                                                        <  1:158231‑M1/216‑53 (MQ=255)                                                                                      cGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                                                                                                                                        >  2:158231‑M1/1‑164 (MQ=255)                                                                                       gACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctg                                                                                                                                                                    >  2:175676‑M1/1‑163 (MQ=255)                                                                                         cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGATAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaacccggcaccggaagtgcgtgtattctggggctggg                                                                                                                                                                  >  2:567421‑M1/1‑161 (MQ=255)                                                                                           ggATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgg                                                                                                                                                                <  1:980616‑M1/251‑93 (MQ=255)                                                                                            gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggg                                                                                                                                                               >  2:590979‑M1/1‑158 (MQ=255)                                                                                            gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgga                                                                                                                                                               <  1:90514‑M1/251‑94 (MQ=255)                                                                                             aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGTACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggcgggtggat                                                                                                                                                              <  2:368823‑M1/251‑95 (MQ=255)                                                                                              acgAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagtgcgtgaattctggggctggtggatg                                                                                                                                                             <  2:940723‑M1/250‑96 (MQ=255)                                                                                                       gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgga                                                                                                                                                     >  1:859227‑M1/1‑147 (MQ=255)                                                                                                          tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcg                                                                                                                                                 >  1:788126‑M1/1‑144 (MQ=255)                                                                                                               tAAAAATTTGCTAAATTTTGCCCATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGCAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattatggggctggtggatggagctggaagcgcagg                                                                                                                                             >  1:953026‑M1/1‑139 (MQ=255)                                                                                                                aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcagga                                                                                                                                            >  1:972326‑M1/1‑138 (MQ=255)                                                                                                                 aaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggatggtggatggagctgggagcgcaggaat                                                                                                                                          >  1:109711‑M1/1‑137 (MQ=255)                                                                                                                       tGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgcaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtt                                                                                                                                    >  1:29‑M1/1‑131 (MQ=255)                                                                                                                        gCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctcgaagcgcaggaatcccgttt                                                                                                                                   >  2:713284‑M1/1‑130 (MQ=255)                                                                                                                           aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACATGAGATAGCAATGACGTTACCGAGTGTACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                                  <  2:952898‑M1/249‑123 (MQ=255)                                                                                                                            aaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACAGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacc                                                                                                                               >  2:49285‑M1/1‑126 (MQ=255)                                                                                                                               tttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggcagcgcgggaatcccgttttacctac                                                                                                                            >  2:201126‑M1/1‑123 (MQ=255)                                                                                                                                ttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaacccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaataccgttttacctac                                                                                                                            >  2:628885‑M1/1‑122 (MQ=255)                                                                                                                                ttGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac                                                                                                                            <  2:773075‑M1/250‑129 (MQ=255)                                                                                                                                 tGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac                                                                                                                            >  1:837573‑M1/1‑121 (MQ=255)                                                                                                                                  gCCAATTTGGTAAATCAGTTGCATCACAACAGGAGATAGCAATGACGTTTCCGATTGGCCACCCGAATAGCAGATTGATCAAAGAATTTACCGCACTCGGCCCGTATATTCGTGAAGGTAagtgcaaagatgatcgattttatgtcgagtgtctggctgtaggcgtcaacgtgaatccggcaccggaattgcgtgaattctggggctgggggatggagctggaagcgcaggaatctcgttttacctacag                                                                                                                          <  2:812121‑M1/250‑131 (MQ=255)                                                                                                                                         tGGTAAAACAGTTTCATCACAACAGGAGATAGCAAGGACGTTACCGAGTGGACACCCGAAGAGCATATTGATCAAAAAATTTACCGCACTAGGCCCGTAGATTCGTGAATGTAagtgcaaagataatcgagtctttttcgattgtctggcggtatgcggcaacgtgaaaccggccccggaagtgcgtgaattctttggctggtggatggagctggaagcgcaggaatcccgtgttacctacagtaaccag                                                                                                                   <  2:913684‑M1/250‑138 (MQ=255)                                                                                                                                            tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttg                                                                                                               >  1:690075‑M1/1‑110 (MQ=255)                                                                                                                                                      gCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcg                                                                                                      <  1:555558‑M1/250‑151 (MQ=255)                                                                                                                                                          cacaACAGGAGATAGCAAGGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaa                                                                                                 <  2:696006‑M1/251‑156 (MQ=255)                                                                                                                                                             aacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                               >  1:870598‑M1/1‑93 (MQ=255)                                                                                                                                                             aacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                               >  1:683191‑M1/1‑93 (MQ=255)                                                                                                                                                                aGGAGATAGCAATGACGTTACCGAGTGGACGCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAATGTAagtgcaaagataatcgattcattttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgggggatggagctggaagcgcatgaatcccgttttacctacagttaccagtttgggctgttcgataaagcaggc                                                                                           <  1:499493‑M1/251‑162 (MQ=255)                                                                                                                                                                aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggc                                                                                           >  1:16855‑M1/1‑90 (MQ=255)                                                                                                                                                                 ggAGATAGCAATGACGTTAACGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcatcgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttttcctacagttaccagtttggtcttttcgataaagcaggcg                                                                                          <  2:852810‑M1/251‑163 (MQ=255)                                                                                                                                                                   agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgac                                                                                        <  1:119208‑M1/251‑165 (MQ=255)                                                                                                                                                                   agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgg                                                                                                              >  2:175799‑M1/1‑87 (MQ=255)                                                                                                                                                                   agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgg                                                                                                              <  1:175799‑M1/229‑143 (MQ=255)                                                                                                                                                                   aaataGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAACAAATTTACCGCACTAGGCCCGTAAATTCGTGAAGGTAagtgcaaagataatcgattctttttagattgtctggctgtatgcgtcaaagtgaaaccggcaccggaagtgcgtgaattctggggctggtggagggagctggaagcgcaggaatcccgttttaactacagttaccaatttggtctgttcgataaagcaggcgac                                                                                        >  2:656906‑M1/3‑87 (MQ=255)                                                                                                                                                                    gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgacaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgact                                                                                       >  2:777170‑M1/1‑86 (MQ=255)                                                                                                                                                                              aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgctttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgt                                                                             >  1:876560‑M1/1‑76 (MQ=255)                                                                                                                                                                               cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgtt                                                                            <  1:30840‑M1/251‑177 (MQ=255)                                                                                                                                                                                gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttc                                                                           >  1:535529‑M1/1‑74 (MQ=255)                                                                                                                                                                                  tACCCAGTTGACACCAGAAGTGCAGATTGATCAAAAAAATTACCGCAATAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtaaacgtgaaaccggcaccggaagtgcgtgaattctggggcgggtggatggagctggaagtgcaggaatcccgttttacctacagtgaccagtttggtctgttcgataaagcaggcgactggtagagtctgcc                                                                          <  2:81654‑M1/250‑179 (MQ=255)                                                                                                                                                                                       aGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactgtaagagtgttccggtaa                                                                     >  1:169951‑M1/1‑67 (MQ=255)                                                                                                                                                                                           gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagac                                                                <  1:174371‑M1/251‑189 (MQ=255)                                                                                                                                                                                             cacCCGAATAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgttaattctggtgctggtggatggagctggaagcgcaggaatcccgttttacctacagtgaccagtttggtctgtgcgataaagcaggcgactggaagagtgttccggtaaaagacac                                                              <  1:403610‑M1/251‑191 (MQ=255)                                                                                                                                                                                               cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgcctggaagagtgttccggtaaaagacactg                                                            >  1:193488‑M1/1‑59 (MQ=255)                                                                                                                                                                                                   aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagggttccggtaaaagacactgaagg                                                        >  2:177026‑M1/1‑55 (MQ=255)                                                                                                                                                                                                   aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtattcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataatgcaggcgactggaagcgtgttccggtaaaagacactgaagg                                                        >  1:473006‑M1/1‑55 (MQ=255)                                                                                                                                                                                                         atattGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataagctattctttttcgatggtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggtctggtggattgagctggaagcgcaggaatcccgttttacctacagttaccagtgtggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttga                                                  <  2:602794‑M1/249‑203 (MQ=255)                                                                                                                                                                                                         aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcgttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattcttgggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttga                                                  <  1:993100‑M1/251‑203 (MQ=255)                                                                                                                                                                                                           aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaac                                                >  1:429379‑M1/1‑47 (MQ=255)                                                                                                                                                                                                           aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaac                                                <  1:506868‑M1/251‑205 (MQ=255)                                                                                                                                                                                                              gttcaAAAAAATTACCGCACTAGGCCCGTATAGTGGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgtagctggaagcgcatgaatctcgttttacctacagtttccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgatgtggttgaactac                                             <  2:620431‑M1/249‑208 (MQ=255)                                                                                                                                                                                                                  aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGACGGTAagggcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttaaacgactgg                                          >  1:284101‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                  aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggg                                         >  1:247136‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                   aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggttggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaagggcactgaagtggttgaacgactgg                                          >  2:73588‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                        attaccgacctaggcccgtATATTCGTGAAGGTAagttcacagataaacgattctttttcgtttgtctggcggtctgcgacaacgtgaaaccggcaccggaagtgcgtgaattctggggcttgtggatggagctggaagcgcaggaagcccgtttttcctacagttaccagtttggtctgttcgagaaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgaatggagcatac                                   <  2:569526‑M1/242‑218 (MQ=255)                                                                                                                                                                                                                               cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgt                            <  1:255807‑M1/251‑225 (MQ=255)                                                                                                                                                                                                                                aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgt                           <  1:518542‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                   aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagt                        >  2:807165‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                        cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcac                   >  1:862509‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                        cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcac                   >  1:90459‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                             aatacgcatcccaggtaATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCGCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATcaccac               <  2:973314/238‑1 (MQ=2)                                                                                                                                                                                                                                                        taggtaATGGCTCCAACTTAGTGATAGTGTTTTATGTTCAGATAATGCCCGATGTCTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTg    <  1:962137/250‑1 (MQ=11)                                                                                                                                                                                                                                                         aggtaATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGGGAc  >  1:409323/1‑251 (MQ=11)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           CCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGAC  >  NC_000913/257659‑258157

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