Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_bme_pgi	2,911,480	Δ12 bp	coding (2161‑2172/2235 nt)	relA ←	(p)ppGpp synthetase I/GTP pyrophosphokinase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	2911480	2911491	12	77 [0]	[0] 74	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 2911479	0 (0.000)	70 (1.230)	49/282	0.1	100%	coding (2173/2235 nt)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
	?	NC_000913_3_bme_pgi	2911492 =	0 (0.000)					coding (2160/2235 nt)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase

TTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCT  >  NC_000913_3_bme_pgi/2911345‑2911636                                                                                                                                                    |                                                                                                                                                 ggatcataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTg                                                                                                                                                <  2:446046‑M2/3‑1 (MQ=255) ggatcataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttaggcctaactcccgtgcaaccgacgcgcctcgataacatccggcacctggttgtgtttacccagcacTTg                                                                                                                                                >  1:879109‑M2/148‑150 (MQ=255)     cgtaaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacagtgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgaaaacatccggaacctggtttagtttacacagcacTTGCAgg                                                                                                                                            >  2:531773‑M2/144‑150 (MQ=255)     cataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacctccggcacctggttgagtttacccagcacTTGCAgg                                                                                                                                            >  2:149420‑M2/144‑150 (MQ=255)     cataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAgg                                                                                                                                            >  2:755921‑M2/144‑150 (MQ=255)     cataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAgg                                                                                                                                            >  2:475999‑M2/144‑150 (MQ=255)     cataaccctttcctcaaaccgctatcatatgtagataaagtatatatcaaactacattgtagatacgatcaaatttcggactaactcacgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagtacTTGCAgg                                                                                                                                            <  1:724933‑M2/7‑1 (MQ=255)      ataaccctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattatagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGt                                                                                                                                           <  2:631374‑M2/8‑1 (MQ=255)         accctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGt                                                                                                                                        >  2:245114‑M2/140‑150 (MQ=255)           cctttcctcaaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTaa                                                                                                                                      <  1:379745‑M2/13‑1 (MQ=255)                cctcaaaccgctatcatatgtagatactgtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTGGTAAAtctc                                                                                                                                 <  1:26729‑M2/18‑1 (MQ=255)                   caaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcttcgataacatccggcaccgggtggaggttaccccgcacGTGCAGGGTGTAAATCtcact                                                                                                                              >  2:251318‑M2/130‑148 (MQ=255)                   caaaccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagaaaatttcggcctacctcccgtgcaaccgacgcgcgtcgataacatccggcacatggttgagttaacccagcacTTGCAGGTTGTAAATCTCAAt                                                                                                                              >  2:645950‑M2/130‑150 (MQ=255)                       ccgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGCAAATCTCAATGGTc                                                                                                                          >  1:1735‑M2/126‑150 (MQ=255)                       ccgctatcatatgtagatacagtatatatcaatctacattgtagataagagcaaatttcggcctaactcccgtgcaacagacgcgcgacgataacatccggcccctggttgagtttacccagcccTTGGAGGTTGTTCATCTCAATGGTc                                                                                                                          >  2:656815‑M2/126‑150 (MQ=255)                        cgctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCa                                                                                                                         <  1:54828‑M2/26‑1 (MQ=255)                         gctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCAt                                                                                                                        >  1:744386‑M2/124‑150 (MQ=255)                         gctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgaagcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCAt                                                                                                                        >  1:563648‑M2/124‑150 (MQ=255)                          ctatcatatgtagatacagtatatatcaatctacattgtagatacgagcaattttcggcctaactcccgtgcaaccgacgcgcgtcgatatcatccggcagcgggttgagtttacgcagcacTGGCAGGTTGTAATTCTCAATGGTCATg                                                                                                                       >  2:553452‑M2/123‑150 (MQ=255)                              catatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcgtcctaactcccttgcaaccgacgcgcgtcgataacatccggcacctgattgagtttacccagcacTTGCAGGTAGTAAATCACAATGGTCAGGACGa                                                                                                                   >  2:186526‑M2/119‑150 (MQ=255)                                tatgtagatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATg                                                                                                                 <  2:147274‑M2/34‑1 (MQ=255)                                      gatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGcc                                                                                                           <  1:770313‑M2/40‑1 (MQ=255)                                      gatacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGcc                                                                                                           <  1:684679‑M2/40‑1 (MQ=255)                                        tacagtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatacggcacctggttgagtttacccagcacTTGCAGGTTGTAATTCTCAATGGTCATGTCGATGGTCGCCAg                                                                                                         >  1:235567‑M2/109‑150 (MQ=255)                                           agtatatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTg                                                                                                      <  2:253520‑M2/45‑1 (MQ=255)                                              atatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgaggcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTg                                                                                                   <  1:951643‑M2/48‑1 (MQ=255)                                              atatatcaatctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTg                                                                                                   <  1:503417‑M2/48‑1 (MQ=255)                                              ataaataaatctacattgtatatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcggcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTtgccg                                                                                                   <  2:533116‑M2/48‑3 (MQ=255)                                                       tctacattgtagatacgagcaaatatcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTc                                                                                          >  1:497063‑M2/94‑150 (MQ=255)                                                       tcaacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTc                                                                                          >  1:397876‑M2/94‑150 (MQ=255)                                                        ctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcccTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCg                                                                                         >  2:690574‑M2/93‑150 (MQ=255)                                                        ctacattgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGCAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCg                                                                                         <  2:1735‑M2/58‑1 (MQ=255)                                                          acatcgtagatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCt                                                                                       <  2:715820‑M2/60‑1 (MQ=255)                                                          aaattgtagatactagcaaatatcggcctatctcccgtgctaccgacgcgagtcgataacatccggcacctggttgagtttacccagctcTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCt                                                                                       <  1:889777‑M2/60‑1 (MQ=255)                                                                 agatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTg                                                                                >  1:164590‑M2/84‑150 (MQ=255)                                                                 agatacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTg                                                                                <  2:297691‑M2/67‑1 (MQ=255)                                                                    tacgagcaaatttcggcctaactcccgtgcaaccgccgcgagtcgataacatccggcacctgtttgagtttacccagcacTTGCAGGTTGTAAATCTCAAGGGTCATGTCGATGGTCGCCAGGTGCTGGTTGGTGTCGCTACGGCTGGCa                                                                             >  1:176113‑M2/81‑150 (MQ=255)                                                                    tacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCa                                                                             >  1:917028‑M2/81‑150 (MQ=255)                                                                    tacgagcaaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCa                                                                             >  1:130746‑M2/81‑150 (MQ=255)                                                                          caaatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCa                                                                       >  2:642112‑M2/75‑150 (MQ=255)                                                                            aatttcggcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAg                                                                     >  1:878434‑M2/73‑150 (MQ=255)                                                                                   gcctaactcccgtgcacccgccgcgcgtcgataacatcacgcacctggttgagtttacccagcacTTGCAGGTTGTAACTCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTc                                                              <  1:251318‑M2/85‑1 (MQ=255)                                                                                   gcctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTc                                                              <  1:776148‑M2/85‑1 (MQ=255)                                                                                   gccgaactcccgtgcaaccgacgcgcgtcgagaacatccggcaccgggttgaggttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTc                                                              <  1:645950‑M2/85‑1 (MQ=255)                                                                                    cctaacttccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAACTCTCAATGGTCATGTCGCTGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTtcc                                                             >  1:424862‑M2/65‑149 (MQ=255)                                                                                    cctaactcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCa                                                             <  1:245114‑M2/86‑1 (MQ=255)                                                                                      aaacgcccgtgcaaccgacgcccgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCAcc                                                           <  1:529265‑M2/88‑1 (MQ=255)                                                                                         ctcccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTAACcttg                                                        >  2:439006‑M2/60‑149 (MQ=255)                                                                                         ctcccgtgcaaccgacgcgcgtcgataacatccggcacctgggtgaggttacccagcacTTGCAGGGTGTAAAGCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTtc                                                        <  1:209673‑M2/91‑1 (MQ=255)                                                                                            ccgtgcaaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCg                                                     <  2:274704‑M2/94‑1 (MQ=255)                                                                                                 caaccgacgcgcgtcgatatcatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCTAGCACGTTCACCTTCTCGTTGGc                                                >  2:426047‑M2/52‑150 (MQ=255)                                                                                                 caaacgacgagcgtcgataacatccggcacctggttgagtatacccagcacTTTCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCTAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGc                                                <  1:531773‑M2/99‑1 (MQ=255)                                                                                                  aaccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTGTGGTGTCGCTTCGGCTGGCAACGCCACGCACGTTCACCTTCTCGTTGGCg                                               >  1:446966‑M2/51‑150 (MQ=255)                                                                                                    ccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCgag                                             <  1:755921‑M2/102‑1 (MQ=255)                                                                                                    ccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCgag                                             <  1:475999‑M2/102‑1 (MQ=255)                                                                                                    ccgacgcgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCgag                                             >  2:410102‑M2/49‑150 (MQ=255)                                                                                                    ccgacgcgcgtcgataacagccggcacctggttgagattacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGGCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCgag                                             <  1:149420‑M2/102‑1 (MQ=255)                                                                                                          cgcgtcgataacatccggcccctggttgagtttacccaccacTTGCAGGGTGTAAATCTCACTGGGCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGt                                       <  2:559632‑M2/108‑1 (MQ=255)                                                                                                          cgcgtcgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGGTGGCGAGAATGGt                                       >  1:342232‑M2/43‑150 (MQ=255)                                                                                                               cgataacatccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGa                                  <  1:289049‑M2/113‑1 (MQ=255)                                                                                                                       tccgtcacctgtttgagtttacccagctcTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCTGTTTCTGTTTGTTGTCGCTACGGCTGGCAACGCCAAGCACGTGCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGt                          <  2:397876‑M2/121‑1 (MQ=255)                                                                                                                       tccggcacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGt                          <  2:497063‑M2/121‑1 (MQ=255)                                                                                                                         ctgccactggttgagtttacccagcacTTGCAGGTTTTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGGGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTaa                        <  2:82667‑M2/123‑1 (MQ=255)                                                                                                                         cggcacctggttgagtttacccatcacTTGCAGGTTGGAAATCTCAATGGTCATGTCGATGGGCGCCAGTTGCTGTTTGGTGGCTCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTaa                        <  1:416508‑M2/123‑1 (MQ=255)                                                                                                                            cacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGGTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTaacaa                     >  2:411170‑M2/25‑150 (MQ=255)                                                                                                                            cacctggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTaacaa                     <  1:855807‑M2/126‑1 (MQ=255)                                                                                                                                tggttgagtttacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCa                 >  1:186097‑M2/21‑150 (MQ=255)                                                                                                                                          tacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTAGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCAt       >  1:946177‑M2/11‑150 (MQ=255)                                                                                                                                          tacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGACAAGCACTTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCAt       >  1:638397‑M2/11‑150 (MQ=255)                                                                                                                                          tacccagcacTTGCAGGTTGTAAATCTCAATGGTCATGACGATGGTCGCCAGTTGCTGTTTCGTGTCGCTACGCCTGGCAACGCCAAGCACGATCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCAt       >  1:826510‑M2/11‑150 (MQ=255)                                                                                                                                            cccagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACACCCCACTACGATCATTa     >  1:603368‑M2/9‑150 (MQ=255)                                                                                                                                              cagcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATGACGTAACAACCCACTACGATCATTAGc   >  1:30106‑M2/7‑150 (MQ=255)                                                                                                                                               agcacTTGCAGGTTGTAAATCTCAATGGTCCTGGCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCt  <  2:563648‑M2/145‑1 (MQ=255)                                                                                                                                               agcacTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGGCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGAGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCt  <  2:744386‑M2/145‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCT  >  NC_000913_3_bme_pgi/2911345‑2911636

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