Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	402,391	Δ723 bp	coding (645‑1367/1416 nt)	phoA →	alkaline phosphatase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	402391	403113	723	43 [0]	[1] 44	phoA	alkaline phosphatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 402390	0 (0.000)	43 (0.760)	40/480	0.4	100%	coding (644/1416 nt)	phoA	alkaline phosphatase
	?	NC_000913	403114 =	0 (0.000)					coding (1368/1416 nt)	phoA	alkaline phosphatase

CGGTGTCAAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAG  >  NC_000913/402145‑402620                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       cGGTGTCAAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCCGAAAAGGATCGATTACCGAccaga                                                                                                                                                                                                                                   <  2:62825‑M1/251‑6 (MQ=255)  ggTGTCAAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATTCTACGGTCCGAGCGCTACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagac                                                                                                                                                                                                                                  <  1:526225‑M1/251‑7 (MQ=255)        aaaaCCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCCCCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatct                                                                                                                                                                                                                            <  1:409910‑M1/251‑13 (MQ=255)         aaaCCTATAACGGCGCGCTGGGCGGCGATCTTCACGAAAAAGATCACCCAACGCTTCTGGCAATGGCAAAAGCCGCAGGTATGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAAAGTCCGGGTAACGCTATGGAAAAAGGCGGAAAAGGATCGATTACAGAccagaccgatctc                                                                                                                                                                                                                           >  2:588947‑M1/1‑238 (MQ=255)               ataACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctac                                                                                                                                                                                                                     >  1:315931‑M1/1‑232 (MQ=255)                    ggCGCGCTGGGCGTCGATATTCACGCAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACATCGCGCAACTGCTACGGTCCGACCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGCAAAAGGCGGAAAAGGATCGATTACAGAccagaccgatctattatacaccat                                                                                                                                                                                                                >  2:26091‑M1/1‑227 (MQ=255)                    ggCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGTATCGATTACCGAccagaccgatctcttctacaccat                                                                                                                                                                                                                <  2:358530‑M1/251‑25 (MQ=255)                         gcTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAACAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaa                                                                                                                                                                                                             >  2:404092‑M1/1‑222 (MQ=255)                           tGGGCGTCGATATTCACGAAAAAGACCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGCAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAACTGCTACGGTCCGCGCGCGCCCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagcc                                                                                                                                                                                                         >  2:518801‑M1/1‑220 (MQ=255)                                  cGATATTCACGAAAAAGATCACCCAACGATTCGGGAAATGGCAAAAGCCGCAGGTCTGGCGACCTGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctgg                                                                                                                                                                                                  <  2:17541‑M1/251‑39 (MQ=255)                                     tatTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggc                                                                                                                                                                                               >  1:29435‑M1/1‑210 (MQ=255)                                       ttCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctg                                                                                                                                                                                             >  1:533599‑M1/1‑208 (MQ=255)                                                        ccAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAAGGTTTCTACCGCAGAGTTGCAGGTTGCCCCGCCCGCTGCGCTGGGGGCACATGTGCCCTCGCGCAAATGCTAGGGTCCGAGCGCGACCAGGGAAAAAGGCCCGGGAACCGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagccggttctcttctacaccttaaaagccgctctggggctgaaataaaaccgcccccg                                                                                                                                                                            >  1:433939‑M1/1‑191 (MQ=255)                                                         cAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccgg                                                                                                                                                                           >  1:133433‑M1/1‑190 (MQ=255)                                                                tcttgaaatGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtga                                                                                                                                                                     <  1:34863‑M1/246‑68 (MQ=255)                                                                  tGGAAATGGCAAACGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggccgtgaatt                                                                                                                                                                  >  2:307415‑M1/1‑181 (MQ=255)                                                                   ggAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaa                                                                                                                                                                    <  1:129595‑M1/248‑69 (MQ=255)                                                                            caaagACGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGCTGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAACTGCTACGGTCCGAGCGCGCCCAGTGAAAAATGTCCGGGTAACGCTCTGGCAAAAGGCGGACAAGGCTCGCTTACCGAccagaccgatctcttctacaccatacaagccgctctggggctgaaataaaaccgcgcccggcagtgacttttcgctgccg                                                                                                                                                        >  2:540342‑M1/2‑171 (MQ=255)                                                                               aGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGCACGCCCCGCCCGCTGCGCTGGTGGCACATATGACCTCGCGCAACTGCTACGGTCCGAGCGCGACCCGCGAAAATTGTCCCGGTGAAGATCTGGAAGAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagcagctctggggctgaaataaacccgcgcccggcagtgaattttcgctgccgggg                                                                                                                                                     >  1:270806‑M1/1‑168 (MQ=255)                                                                                   gCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCTCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtt                                                                                                                                                 <  2:410780‑M1/251‑88 (MQ=255)                                                                                   gCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtt                                                                                                                                                 >  1:582761‑M1/1‑164 (MQ=255)                                                                                     aGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttt                                                                                                                                               >  1:569063‑M1/1‑162 (MQ=255)                                                                                                       cGTTTCTACCGCAGAGTTGCATGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAAGGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcgtctacaccatgaaagccgctcgggggctgaaataaggccgcgcctggcagtgaatgttcgctgccgggtggtttttttgctgttagcaaccag                                                                                                                             <  2:587212‑M1/251‑108 (MQ=255)                                                                                                              aCCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaatgaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagactta                                                                                                                        <  1:439818‑M1/249‑113 (MQ=255)                                                                                                                    gagTTGCAGGATGCCACGCCCGCTGAGCTGGGGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaacctcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggaaga                                                                                                                <  2:374182‑M1/251‑121 (MQ=255)                                                                                                                                      cccGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcg                                                                                                       >  1:228648‑M1/1‑113 (MQ=255)                                                                                                                                      cccGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcg                                                                                                       <  2:228648‑M1/242‑130 (MQ=255)                                                                                                                                                       acaTGTGACAGCGCGCAAATGCTATGGTCCGAGCGCGACCAGTTAAAAATGTCCGAGTAACGCTCGGGAAAAAGGCGGAAAAGGATCGATTTCCGTccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaacccgcgcccggcagtgaattttcgctgctgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggctcatacgctcatggttaaaacatgaag                                                                             <  2:542252‑M1/251‑156 (MQ=255)                                                                                                                                                          tgtgACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagg                                                                          >  1:112037‑M1/1‑93 (MQ=255)                                                                                                                                                           gtgACCTCGCGCAAATGCTACGGGCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagg                                                                          <  2:299331‑M1/250‑159 (MQ=255)                                                                                                                                                                   gcgcAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgct                                                                 <  1:177819‑M1/251‑168 (MQ=255)                                                                                                                                                                     gcAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTACCGCTCTGGAAAAAGGCGGAAAAGGATCGCTTACCGAcaagaccgatctcttctacaccatgaaagccgctctggggatgaactaaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaacccgacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggagctat                                                               >  1:76781‑M1/1‑82 (MQ=255)                                                                                                                                                                     gcAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacacaatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagaattaatggcagatcacgggcgaatacgctcatggttaaaacatgaagagggatggggctat                                                               >  1:253264‑M1/1‑82 (MQ=255)                                                                                                                                                                      cAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatg                                                              <  1:404092‑M1/251‑171 (MQ=255)                                                                                                                                                                           gCTACGGTCCGTGCGCGACCAGTAAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcgtctgctctatgtaagccgctctggggcttaagtaaaaccgcgccctgcagtgaattgtcgctgccgggttgttgttttgctgttagcgaccaggcttaatggcagatcgcgggggcatacgctcatggttaaaacatgaagagggatggtgctatgaaaat                                                         <  1:578859‑M1/251‑176 (MQ=255)                                                                                                                                                                              aCGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaata                                                        <  2:166032‑M1/249‑177 (MQ=255)                                                                                                                                                                              aCGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaagggcagatcacgggcgcatacgctcatggttaaaacatgaagatggatggtgctatgaaaagaa                                                       <  2:129057‑M1/250‑178 (MQ=255)                                                                                                                                                                                             ccAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttataacatgaagagggatggtgctatgaaaataacattactggttccctt                                       >  1:497114‑M1/1‑58 (MQ=255)                                                                                                                                                                                                    aaaaTGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaataacattactggttaccttgcttttc                                >  1:393472‑M1/1‑51 (MQ=255)                                                                                                                                                                                                           ccGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaatccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaattgcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaataacattactggttaccttgcttttcggtctgg                         <  1:49814‑M1/251‑208 (MQ=255)                                                                                                                                                                                                                         ggAAAAAGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaataacattactggttaccttgcttttcggtctggtttttttaacc              <  1:260456‑M1/248‑219 (MQ=255)                                                                                                                                                                                                                             aaaGGCGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgctgttagcaaccagacttaatggcagatcacgggcgcatacgctcatggttaaaacatgaagagggatggtgctatgaaaataacattactggttaccttgcttttcggtctggtttttttaaccaccgtcg       <  2:315931‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                  cGGAAAAGGATCGATTACCGAccagaccgatctcttctacaccatgaaagccgctctggggctgaaataaaaccgcgcccggcagtgaattttcgctgccgggtggtttttttgcggttagcaaccagacttaatgacagatcaagggcgaatacgcacatggttaaaacatgaagagggatggagctaggaaaataacattactggttaccttgcttttcggtctggtttttttaaccaccgacggcgct  >  2:114666‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       CGGTGTCAAAACCTATAACGGCGCGCTGGGCGTCGATATTCACGAAAAAGATCACCCAACGATTCTGGAAATGGCAAAAGCCGCAGGTCTGGCGACCGGTAACGTTTCTACCGCAGAGTTGCAGGATGCCACGCCCGCTGCGCTGGTGGCACATGTGACCTCGCGCAAATGCTACGGTCCGAGCGCGACCAGTGAAAAATGTCCGGGTAACGCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAG  >  NC_000913/402145‑402620

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