Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,233,804	1557 bp→GGTGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT	coding (47‑1603/1650 nt)	pgi →	glucosephosphate isomerase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4233804	4235360	1557	73 [0]	[0] 65	pgi	glucosephosphate isomerase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4233803	0 (0.000)	40 (0.600) +GGTGTAGGCTGGAGCTGCTTCGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGAACTGCAGGTCGACGGATCCCCGGAAT	30/332	1.0	100%	coding (46/1650 nt)	pgi	glucosephosphate isomerase
	?	NC_000913	4235361 =	0 (0.000)					coding (1604/1650 nt)	pgi	glucosephosphate isomerase

CCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTAGCGCTGCGTAACCGTAGCAATACCCCGATTTTGGTTGATGG  >  NC_000913/4233508‑4234095                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     cttcttctggACTGCGCTAGCGCAGGTAGTACATTTATAAACAAAGGGTGAGCGGGGCGGTTGTCAACGATAGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTGGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgt                                                                                                                                                                                                                                                                                                 <  2:358593‑M1/299‑6 (MQ=255)                         ggTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggt                                                                                                                                                                                                                                                                                                   <  2:656606‑M1/275‑4 (MQ=255)                          gtagtaCATTTAGAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTGACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggct                                                                                                                                                                                                                                                                                            <  2:379917‑M1/281‑11 (MQ=255)                          atagtacATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcc                                                                                                                                                                                                                                                                         >  2:121183‑M1/2‑271 (MQ=255)                             gtaCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctata                                                                                                                                                                                                                                                                     >  2:230228‑M1/1‑268 (MQ=255)                              taCATTTAGAAGTAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGGCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                    <  2:839186‑M1/301‑35 (MQ=255)                              taCAATTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                    <  2:13962‑M1/301‑35 (MQ=255)                              aacATTGATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                    <  2:939226‑M1/300‑35 (MQ=255)                                cATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactt                                                                                                                                                                                                                                                                  >  2:462113‑M1/1‑265 (MQ=255)                                   ttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGCGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttct                                                                                                                                                                                                                                                               <  2:252274‑M1/301‑40 (MQ=255)                                        aataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagaga                                                                                                                                                                                                                                                          >  2:640416‑M1/1‑257 (MQ=255)                                        aataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagaga                                                                                                                                                                                                                                                          >  2:161082‑M1/1‑257 (MQ=255)                                          aaaagGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaat                                                                                                                                                                                                                                                        <  2:818995‑M1/300‑47 (MQ=255)                                              gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaatagga                                                                                                                                                                                                                                                    >  2:646627‑M1/1‑251 (MQ=255)                                                                            tcaggtggatttttGATCCAATTCTGGCGGGCAGTAGTTCATCTAATAAATGCTTCAGTTCGAAAAGGGGTTAAACTCAAAATTACGCTAACGGCAATAAAACCATATCATTTTGCTGTGACTGGCGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAGCATCAATCCAACGCAGACCGCTGCCTTGCAGGCACTATggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttagaactgcaggtcgacggatccccgg                                                                                                                                                                                                                      <  2:490991‑M1/295‑81 (MQ=255)                                                                                 cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACgg                                                                                                                                                                                                                                                                                                    <  1:539682‑M1/218‑3 (MQ=255)                                                                                        ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggc                                                                                                                                                                                                                                                                                             <  1:94208‑M1/218‑10 (MQ=255)                                                                                        ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggc                                                                                                                                                                                                                                                                                             <  1:131859‑M1/218‑10 (MQ=255)                                                                                        ntcATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggc                                                                                                                                                                                                                                                                                             <  1:685466‑M1/217‑10 (MQ=255)                                                                                             ccACTCCTGGCGGTAAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACAATACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggt                                                                                                                                                                                                     <  2:103608‑M1/301‑98 (MQ=255)                                                                                                                 ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGAACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaagctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                    <  1:274195‑M1/218‑35 (MQ=255)                                                                                                                 ntcATCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGATAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatac                                                                                                                                                                                                                                                                    <  1:833526‑M1/217‑35 (MQ=255)                                                                                                                  tCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcg                                                                                                                                                                                <  2:589528‑M1/301‑119 (MQ=255)                                                                                                                        aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgcgccacgctttatagc                                                                                                                                                                          <  2:334404‑M1/301‑125 (MQ=255)                                                                                                                        aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctag                                                                                                                                                                                                                                                             <  1:462110‑M1/218‑42 (MQ=255)                                                                                                                        aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttct                                                                                                                                                                                                                                                               >  1:800701‑M1/1‑177 (MQ=255)                                                                                                                        aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTGTCTGTGACGGGCGCTACAATCGTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctag                                                                                                                                                                                                                                                             <  1:101475‑M1/218‑42 (MQ=255)                                                                                                                        aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAAAGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTATCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggagtaggctggagatgcctcgagatttctataatgtctan                                                                                                                                                                                                                                                             >  1:662532‑M1/1‑177 (MQ=255)                                                                                                                         ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTTGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctaga                                                                                                                                                                                                                                                            >  1:453297‑M1/1‑176 (MQ=255)                                                                                                                                cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcacggttaatc                                                                                                                                                                  <  2:623171‑M1/301‑133 (MQ=255)                                                                                                                                  tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTTCggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaa                                                                                                                                                                                                                                                   <  1:569704‑M1/218‑52 (MQ=255)                                                                                                                                  ncaCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaa                                                                                                                                                                                                                                                   <  1:934431‑M1/217‑52 (MQ=255)                                                                                                                                       gcgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcg                                                                                                                                                                                                                                              >  1:349262‑M1/1‑162 (MQ=255)                                                                                                                                               gttttacACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaa                                                                                                                                                                                                                                                   <  1:542832‑M1/203‑52 (MQ=255)                                                                                                                                               ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgtcggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcggtatgtag                                                                                                                                                   >  2:459284‑M1/1‑154 (MQ=255)                                                                                                                                               ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaa                                                                                                                                                                                                                                                   >  2:542839‑M1/1‑154 (MQ=255)                                                                                                                                               ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACAGCTGCCTGGCAGGCACTACggtgtaggctggagctgctccgaagttcctatactttctagaaaataggaacttcgaactgcaggtagacggatccccggaatgcacgaccaatggtctgattaaccgatataaagcgtggcgccgttaatcatcttcgttatgtaa                                                                                                                                                   >  2:432159‑M1/1‑154 (MQ=255)                                                                                                                                                  ntaCACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtc                                                                                                                                                                                                                                   <  1:823706‑M1/217‑68 (MQ=255)                                                                                                                                                     cacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgan                                                                                                                                                                                                                                >  1:799790‑M1/1‑148 (MQ=255)                                                                                                                                                       cgcaaCATTACGCTAACGGCACTAAAACCAACACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAAGCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggcaggata                                                                                                                                           <  2:724035‑M1/299‑156 (MQ=255)                                                                                                                                                           aCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgcactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatctgtaggccggataagg                                                                                                                                        >  2:929177‑M1/1‑142 (MQ=255)                                                                                                                                                                            cTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccgg                                                                                                                      >  2:206362‑M1/1‑125 (MQ=255)                                                                                                                                                                                    aTCACATTATTCTGTGACTGGCGCTACAATCGTCCAAAGTCACAATTCTCAAAATCAGAATAGTATTGCTAATGAAAAACATCAATCCAACGCTGCCCGCTGCCTGGCAGGCAATACtgtgtaggcaggagctgcttcgaagttcctatactttctagagaataggaacttcgaactggaggtcgacggatccccggaatgctcgaccaatggtctga                                                                                                                                                                                                 <  1:429501‑M1/218‑102 (MQ=255)                                                                                                                                                                                                        gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtgg                                                                                                                                                                             >  1:191543‑M1/1‑97 (MQ=255)                                                                                                                                                                                                                accttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGACAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatc                                                                                  <  2:336393‑M1/299‑213 (MQ=255)                                                                                                                                                                                                                    cacacagtcacAATTCTCAAAATCAGACGAGTATTTCTATGTATAACCATCACTCCAACGCGAACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatgctttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaattgtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggc                                                                              <  2:209323‑M1/296‑217 (MQ=255)                                                                                                                                                                                                                          gTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatttaggccggctaaggcgttcacgccgcatccggcaacctatgccttatgccacgcggtcgcgtcttctcaggcctacag                                                                        >  2:217639‑M1/1‑79 (MQ=255)                                                                                                                                                                                                                                nttcTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaggcgtggcgcggttaatcatcgtcgatatgt                                                                                                                                                     <  1:926972‑M1/217‑146 (MQ=255)                                                                                                                                                                                                                                aTTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgt                                                                                                                                                     <  1:230224‑M1/218‑146 (MQ=255)                                                                                                                                                                                                                                  tctcAACATCAGAAGCGTATTGCTAATGCAAAACATCACTCCTCCCCCGACCGCTGCATGGCCGGCACTTCtgtgtaggctggagctgattcgaagttccgatactttctagagaataggatcttcgaactgcaggtcgacggatccccggaatgctcgaccactgttctgattaaccgctataaagcgtggcgcggttaatcatcgtcgctatgtaggccggataaggcgttcccgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgc                                                                <  2:326839‑M1/301‑231 (MQ=255)                                                                                                                                                                                                                                  tctcAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgc                                                                <  2:191546‑M1/301‑231 (MQ=255)                                                                                                                                                                                                                                    tcAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggc                                                                                                                                                 <  1:121182‑M1/218‑150 (MQ=255)                                                                                                                                                                                                                                        aaTCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccgga                                                                                                                                             >  1:19091‑M1/1‑65 (MQ=255)                                                                                                                                                                                                                                              aagtGTATGACTAATGAACAACATCCATCCAACGCAGCCCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacat                                                    <  2:928220‑M1/301‑243 (MQ=255)                                                                                                                                                                                                                                              aagaGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgtcggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataagg                                                                                                                                        <  1:459281‑M1/217‑159 (MQ=255)                                                                                                                                                                                                                                              aagaGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacat                                                    <  2:19091‑M1/301‑243 (MQ=255)                                                                                                                                                                                                                                              aagaGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtac                                                      <  2:816250‑M1/299‑241 (MQ=255)                                                                                                                                                                                                                                              aagaGTATTGCTAATGAAAAACAGCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtac                                                      <  2:668813‑M1/299‑241 (MQ=255)                                                                                                                                                                                                                                                   taatgctAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaatatgaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtat                                               <  2:800727‑M1/298‑248 (MQ=255)                                                                                                                                                                                                                                                        cTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggc                                          <  2:651860‑M1/301‑253 (MQ=255)                                                                                                                                                                                                                                                              aaaaaCATCAAGCCAACGCAGACCGCTGCCGGGCAGGCACTACggtgcagtctggagatgctgcgcagttcctatactttcaagagaataggtagttcgaaccgcaggtcgacgcagccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggtgaatcatcgtggatatggaggtcggataaggcgtccacgccgcctc                                                                                                                         <  1:410435‑M1/216‑174 (MQ=255)                                                                                                                                                                                                                                                                   tatCAATCCAGCGCAGACCGCTGCTTGGCATGCACCACgtttggcgctggagttgcttagaagttattatacttgttatagaataggaacttcgaactgcaggtcgacggatccccgtaatgctcgaccaatggtctgattaaccgctaaaaagagtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacat                               <  2:662543‑M1/300‑264 (MQ=255)                                                                                                                                                                                                                                                                   cATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaa                                                                                                                   >  1:490116‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                                                   cATCAAGCCAACGCAGACCGCTGCCGGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttatagagaataggaacgtcgaacggcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggc                                          <  2:490119‑M1/290‑253 (MQ=255)                                                                                                                                                                                                                                                                        aTCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcccgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggtaattacttcct                               >  2:766784‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                                                        aTCCAACGCAGACCGCTGCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggccattaatacatagca                           >  2:140856‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                                                           cAACGCAGACCGCTGCCTGGCAGGCACGACggtgtaggctggagatgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgacaaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacatagcacgat                       <  2:281744‑M1/301‑272 (MQ=255)                                                                                                                                                                                                                                                                           cAACGCAGACCGCTGCCTGGCAGGCAATACgaggtaggcgtgagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacatagcacgat                       <  2:490032‑M1/301‑272 (MQ=255)                                                                                                                                                                                                                                                                             aCGCAGACCGCGTCCTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacatagcacgattg                     <  2:453300‑M1/301‑274 (MQ=255)                                                                                                                                                                                                                                                                                          ccTGGCAGGCACTACggtgtaggctggagctgcttcgaagttcctatacttcttagagaataggaacttcgaactgaaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacatagcacgattgattaaataacct         <  2:799815‑M1/300‑286 (MQ=255)                                                                                                                                                                                                                                                                                              gcaggcaCTACggtgtaggctggagctgcttcgaagttcctatactttctagagaatagcaacttcgaactgcaggtcgacggatccccggaatgctctaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgt                                                                                        <  1:439316‑M1/217‑207 (MQ=255)                                                                                                                                                                                                                                                                                               caggcaCTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgt                                                                                        <  1:640405‑M1/216‑207 (MQ=255)                                                                                                                                                                                                                                                                                                 ggcaCTACggtgtaggctggagctgcttcgaagttcctatactttctagagaataggaacttcgaactgcaggtcgacggatccccggaatgctcgaccaatggtctgattaaccgctataaagcgtggcgcggttaatcatcgtcgatatgtaggccggataaggcgttcacgccgcatccggcaaccgatgcctgatgcgacgcggtcgcgtcttatcaggcctacaggtcgatgccgatatgtacatcgtattcggcaattaatacatagcacgattgattaaataaccttaataac  <  2:767705‑M1/300‑293 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     CCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCCTGGCAGGCACTACAGAAACACTTCGATGAAATGAAAGACGTTACGATCGCCGATCTTTTTGCTAAAGACGGCGATCGTTTTTCTAAGTTCTCCGCAACCTTCGACGATCAGATGCTGGTGGATTACTCCAAAAACCGCATCACTGAAGAGACGCTGGCGAAATTACAGGATCTGGCGAAAGAGTGCGATCTGGCGGGCGCGATTAAGTCGATGTTCTCTGGCGAGAAGATCAACCGCACTGAAAACCGCGCCGTGCTGCACGTAGCGCTGCGTAACCGTAGCAATACCCCGATTTTGGTTGATGG  >  NC_000913/4233508‑4234095

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑