Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_bme_pgi	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1978503–1979270	1979278	9–776	68 [2]	[0] 64	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1978502	2 (0.030)	60 (1.080)	48/282	0.1	98.5%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_bme_pgi	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

ATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  NC_000913_3_bme_pgi/1978353‑1978651                                                                                                                                                      |                                                                                                                                                      atattaatgggtCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACAAGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                       <  2:918139‑M1/145‑1 (MQ=255) aTAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAATGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                       <  1:908742‑M1/150‑1 (MQ=255)  tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                                      >  1:200112‑M1/1‑149 (MQ=255)       aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                                 <  2:494483‑M1/150‑7 (MQ=255)         gggTCTTTACACTTATCTAAGATTTTTCCTACATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                               <  2:809307‑M1/150‑9 (MQ=255)          ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                              >  1:413412‑M1/1‑141 (MQ=255)           gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                             >  2:372651‑M1/1‑140 (MQ=255)                  cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGAGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                                      >  1:681160‑M1/1‑133 (MQ=255)                        tCTAAGATTTTTCCTAAATCGACGCAACAGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                    <  2:818756‑M1/130‑4 (MQ=255)                        tCTAAGATTTTTCCTAAATCGACGCAACAGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                    >  1:818756‑M1/1‑127 (MQ=255)                          tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                              >  1:174788‑M1/1‑125 (MQ=255)                              aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                          >  1:754099‑M1/1‑121 (MQ=255)                              aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                          <  2:200112‑M1/150‑30 (MQ=255)                                     cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                   >  2:787972‑M1/1‑114 (MQ=255)                                       aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                                 <  1:441455‑M1/150‑39 (MQ=255)                                           cgccgcaACTGTACTCGTCAATACACGCACATACAACGTAGCGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcaaagttttgttttgggtgtatttcaatcagtaaa                                                                                                             >  2:27513‑M1/4‑108 (MQ=255)                                           cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                                                                             >  2:92628‑M1/1‑108 (MQ=255)                                              cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                          <  1:897660‑M1/150‑46 (MQ=255)                                              agcAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                          <  1:722005‑M1/149‑46 (MQ=255)                                                  aCTGTACTCGTCACTACACGCACATACAACGGATTGGGTCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                      <  1:691838‑M1/150‑50 (MQ=255)                                                    tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                    >  1:818362‑M1/1‑99 (MQ=255)                                                     gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                   <  1:92628‑M1/150‑53 (MQ=255)                                                          cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaa                                                                                              >  1:925086‑M1/1‑93 (MQ=255)                                                          cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTACAATAAAGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagacatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctac                                                                                              >  1:477761‑M1/1‑93 (MQ=255)                                                               cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                         <  1:599596‑M1/150‑63 (MQ=255)                                                                tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                        <  1:93954‑M1/150‑64 (MQ=255)                                                                    cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                    <  2:709512‑M1/150‑68 (MQ=255)                                                                          tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                              >  2:809252‑M1/1‑77 (MQ=255)                                                                           aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatt                                                                             >  2:604746‑M1/1‑76 (MQ=255)                                                                             aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATCAATTACTTAACATAAATGtattagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                           >  1:396951‑M1/1‑74 (MQ=255)                                                                              aCGGAGGGGGGCTGCGTTTTTCAATAATGCTTGATGCAGATCACACAAAACAATCAATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaattttttcgctacgctagcaagcaaaaatgaaa                                                                          <  2:477761‑M1/150‑78 (MQ=255)                                                                                 gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                       >  1:471450‑M1/1‑70 (MQ=255)                                                                                 gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                       >  1:278690‑M1/1‑70 (MQ=255)                                                                                  aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                                      <  1:860200‑M1/150‑82 (MQ=255)                                                                                         cTGCGATTTTCAAAAATGCGTGATGCAGATCACACATAACACTCAATTACTTAACATAAATGtataagtcatacttttgttatgggtgtatttcaatctgttaaaaagtttttcgctacgctagctagcaaaaatgaaacaggaataatc                                                               <  2:733168‑M1/150‑89 (MQ=255)                                                                                          tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              <  2:681160‑M1/150‑90 (MQ=255)                                                                                          tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              >  2:309264‑M1/1‑61 (MQ=255)                                                                                          tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              <  2:389591‑M1/150‑90 (MQ=255)                                                                                             gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                           >  1:91534‑M1/1‑58 (MQ=255)                                                                                             gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                           >  2:903271‑M1/1‑58 (MQ=255)                                                                                             gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaa                                                           >  2:548791‑M1/1‑58 (MQ=255)                                                                                                 ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                       <  2:174788‑M1/150‑97 (MQ=255)                                                                                                  tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                                      >  2:61987‑M1/1‑53 (MQ=255)                                                                                                       aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttatgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                                 >  1:757743‑M1/1‑48 (MQ=255)                                                                                                          gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                              <  1:489469‑M1/150‑106 (MQ=255)                                                                                                               aTGCAGATCACACAAACCACTCAATTACTTAACATAAATGgataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaaggggatgttgcgcacagt                                         >  2:69274‑M1/1‑40 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  1:117778‑M1/150‑117 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  2:91534‑M1/150‑117 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  1:431088‑M1/150‑117 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   >  1:109459‑M1/1‑34 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaaa                                   >  1:297175‑M1/1‑34 (MQ=255)                                                                                                                     aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaacaggaataatcgaaatgggatgttgcgcacaggcaaaat                                   <  1:807140‑M1/150‑117 (MQ=255)                                                                                                                       cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggattaatcgaaatgggatgttgcgcacagtcaaaataa                                 >  2:46736‑M1/1‑32 (MQ=255)                                                                                                                         cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatgggatgttgcgcacagtcaaaataact                               >  1:757351‑M1/1‑30 (MQ=255)                                                                                                                             aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                           >  2:790516‑M1/1‑26 (MQ=255)                                                                                                                                 acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                       >  2:346032‑M1/1‑22 (MQ=255)                                                                                                                                       ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc                 >  1:94107‑M1/1‑16 (MQ=255)                                                                                                                                        tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca                <  1:372651‑M1/150‑136 (MQ=255)                                                                                                                                         aCTTAACATAATTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaattgaaacaggaataatcgaaatgggatgttgcgcacagtcaaattaactcaccgtaaataatcat               >  2:867043‑M1/1‑14 (MQ=255)                                                                                                                                           ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct             <  1:346032‑M1/150‑139 (MQ=255)                                                                                                                                           ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct             <  2:925086‑M1/150‑139 (MQ=255)                                                                                                                                            tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg            >  1:225223‑M1/1‑11 (MQ=255)                                                                                                                                               cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctt         >  2:192567‑M1/1‑8 (MQ=255)                                                                                                                                                cgtattGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCt        <  1:900462/145‑1 (MQ=14)                                                                                                                                                 tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata       >  1:237634‑M1/1‑6 (MQ=255)                                                                                                                                                 tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata       <  1:809252‑M1/150‑145 (MQ=255)                                                                                                                                                      gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataa  >  2:196206‑M1/1‑1 (MQ=255)                                                                                                                                                      gggTAATGACTCCAACTTATTGATAGTGATTTATGTTCAGATAATGCACGATGACTTTGTCATGCAGCTCGCCCGATTTTGAGAACGACAGCGACTTCCGACACAGCAGTGCCAGGTGCTGCCTCAGATTCTGGTTATGCCTGTCAAtta  >  1:648172/1‑149 (MQ=11)                                                                                                                                                      |                                                                                                                                                      ATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  NC_000913_3_bme_pgi/1978353‑1978651

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