Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	AE015451	1,709,370	2640 bp→CGGCCGC		[ppc]	[ppc]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	AE015451	1709370	1712009	2640	38 [0]	[0] 38	[ppc]	[ppc]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1709369	0 (0.000)	38 (0.550) +CGGCCGC	38/586	0.6	100%	intergenic (‑136/‑13)	PP_1504/ppc	conserved protein of unknown function/Phosphoenolpyruvate carboxylase
	?	AE015451	1712010 =	0 (0.000)					intergenic (+1/‑257)	ppc/adk	Phosphoenolpyruvate carboxylase/Adenylate kinase

GAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGCACAAGATGCGCGTAATTCTGCTGGGAGCTCCCGG  >  AE015451/1711716‑1712295                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                              cgatcagcgcggccaggcgcgcttcggactccaggtccaggtgcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAg                                                                                                                                                                                                                                                                                         >  3:482569‑M2/295‑301 (MQ=255)   atcagcgcggccaggcgcgcttcggactccaggtccaggtgcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGcc                                                                                                                                                                                                                                                                                       <  1:482567‑M2/9‑1 (MQ=255)     cagcgcggccaggcgcgcttcggactccaggtccaggtgcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCg                                                                                                                                                                                                                                                                                     <  2:423032‑M2/11‑1 (MQ=255)        gggggccaggcgcgcgtcggattcaaggtcaaggggaaggacaggacaggtacggctgagggtaccggcggcattctgttcccagggctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGt                                                                                                                                                                                                                                                                                  <  2:391907‑M2/14‑1 (MQ=255)                     gcgtgtcgggttgggtccaggggcaggacaggaccggtacggctgagggtaccggcgggattctgtgcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGAc                                                                                                                                                                                                                                                                     <  2:290651‑M2/27‑1 (MQ=255)                             ctccaggtccaggtgcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggagaccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCAtaa                                                                                                                                                                                                                                                             >  2:504451‑M2/267‑300 (MQ=255)                                  ggtccaggtgcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgcggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGtcatgtctggtgacagatcggaagatc                                                                                                                                                                                                                                                        >  4:290719‑M2/262‑288 (MQ=255)                                           tcaggacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgattacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAg                                                                                                                                                                                                                                               <  4:504593‑M2/49‑1 (MQ=255)                                               gacatgaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACg                                                                                                                                                                                                                                           >  3:67241‑M2/249‑301 (MQ=255)                                                    gagcggtacgggtgagggtatcggcggcatgctgttcccagggctgggcgagttctaggatgttcatgatgacctctgttgacggctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAg                                                                                                                                                                                                                                      <  4:598398‑M2/58‑1 (MQ=255)                                                    gaccggtacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTTACGGCCATACAGGCCCGGCACAAGGACGGGGAg                                                                                                                                                                                                                                      >  3:25765‑M2/244‑301 (MQ=255)                                                          tacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATTTCTGGTGACGGCCATACAGGCCCTGCCCCAGGACGGGGacggcggg                                                                                                                                                                                                                                >  2:10457‑M2/238‑300 (MQ=255)                                                          tacggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGt                                                                                                                                                                                                                                <  3:254971‑M2/64‑1 (MQ=255)                                                           acggctgagggtaccggcggcattctgttcccagtgctgggcgagttcacggattttcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGtt                                                                                                                                                                                                                               <  1:67239‑M2/65‑1 (MQ=255)                                                                                                              gattttcatgatgacctctgttgacgcctggaaaggctggattgaaggggagtgctgtgttttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCa                                                                                                                                                                            <  4:10457‑M2/116‑1 (MQ=255)                                                                                                                   tcatgatgacctctgttgacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACttt                                                                                                                                                                       >  3:625064‑M2/181‑301 (MQ=255)                                                                                                                                tgttgacgcctggaaaggctgcactgcatgggtgtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCt                                                                                                                                                          <  4:203739‑M2/134‑1 (MQ=255)                                                                                                                                    gacgcctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCt                                                                                                                                                      <  1:625062‑M2/138‑1 (MQ=255)                                                                                                                                        cctggaaaggctgcactgcatgggagtgctgtgctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCg                                                                                                                                                  >  3:244273‑M2/160‑301 (MQ=255)                                                                                                                                                               gagtgctgtgctgtcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgtgccggcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGcc                                                                                                                           <  1:25765‑M2/165‑1 (MQ=255)                                                                                                                                                                        ctcttccgatctgtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCg                                                                                                                  <  1:11115‑M2/174‑1 (MQ=255)                                                                                                                                                                         ctttcactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGa                                                                                                                 >  4:340825‑M2/127‑301 (MQ=255)                                                                                                                                                                             cactttagtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGatttt                                                                                                             >  3:41353‑M2/123‑300 (MQ=255)                                                                                                                                                                                    gtttggtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGccttctgaccgagatcggaagag                                                                                                      >  3:11115‑M2/116‑290 (MQ=255)                                                                                                                                                                                        ggtttggccgatcacggtgggcctttggcgctggggtgtgggcactcttgggggataccgctcgctcatgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATAccc                                                                                                  <  2:340746‑M2/190‑1 (MQ=255)                                                                                                                                                                                         gtttggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATacccc                                                                                                 >  3:580929‑M2/111‑300 (MQ=255)                                                                                                                                                                                             ggccgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGAttt                                                                                             >  3:420228‑M2/107‑301 (MQ=255)                                                                                                                                                                                                cgatcacggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgcgcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTcc                                                                                          <  1:580927‑M2/198‑1 (MQ=255)                                                                                                                                                                                                      cggcttgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTTCGGTCtt                                                                                    >  4:780179‑M2/98‑301 (MQ=255)                                                                                                                                                                                                           tgccatcgtcgctcggctgtgggcactcttgggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGccc                                                                               >  3:345175‑M2/93‑301 (MQ=255)                                                                                                                                                                                                                                gcgacactggggggataccgtgcgctcgtgcacaaggcgcgctaggactgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATCCAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAAt                                                          <  2:779899‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                          ggggataccgctcgctcctgcacagggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACaaa                                                >  4:783584‑M2/62‑301 (MQ=255)                                                                                                                                                                                                                                                       ctgcggatctagggcgggctgaggctgatagggatgaaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGca                                   <  2:395977‑M2/253‑1 (MQ=255)                                                                                                                                                                                                                                                                 agggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAatttgaggagcaagatcggaag                         >  4:396074‑M2/39‑291 (MQ=255)                                                                                                                                                                                                                                                                   ggcgcgctgaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGCACAAGATGCGCGt                       <  1:41352‑M2/265‑1 (MQ=255)                                                                                                                                                                                                                                                                          ggaggctgatagggaggcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGCACAAGATGCGCGTAATTctg                <  2:783304‑M2/272‑1 (MQ=255)                                                                                                                                                                                                                                                                           gaggctgatagggatgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATGTAATTTAATACAAAAATTTGAGGAGCACAAGGTGCGGGTAattcagg               >  2:264934‑M2/29‑298 (MQ=255)                                                                                                                                                                                                                                                                                        gtgcaacgcggccgcGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGCACAAGATGCGCGTAATTCTGCTGGGAGCTCCCgg  <  4:265000‑M2/286‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                              GAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCGGGCTGCGCAACACCGGCTAAGGGCAAGCCCGGGTCATGTCTGGTGACGGCCATACAGGCCCGGCACAAGGACGGGGAGGGCGGTTGCGCCAGGCGCAACCGCCTGCCCGGTCAGCCGCAAGTGGCGCATTCCTGCAACTTTGGGACGCTTGTCTGGCTGCCGGGCGCTGTGTATCTTGAGCAGCCTTCTGACCGATTTATGGTCATACCCGATTTTCCGGACTTGGCCCTGGTCGCCGAATCCATTGAATTTCATACAAAAATTTGAGGAGCACAAGATGCGCGTAATTCTGCTGGGAGCTCCCGG  >  AE015451/1711716‑1712295

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