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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979268	1979278	11–776	76 [0]	[1] 78	insB1–insA	insB1, insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	71 (1.000)	53/266	0.2	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644                                                                                                                                           |                                                                                                                                               ctgtacaCTTATCTAAGATGTTTCCTAAACCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  1:782856‑M1/140‑6 (MQ=255) cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTACGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  2:369199‑M1/143‑6 (MQ=255)  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtaaaag                                                                                                                                          >  2:778033‑M1/1‑137 (MQ=255)     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  1:1269706‑M1/1‑134 (MQ=255)     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTT‑AACAAGGCCAGGGTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  1:606547‑M1/1‑134 (MQ=255)       acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                     >  1:1118452‑M1/1‑132 (MQ=255)               aaGATTTTTCCTAAATCGACGCAACTGTGCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                             >  1:874072‑M1/1‑124 (MQ=255)                  aTTTTTTCTATATCGACGTACATTTACTCGTTACTACACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATG‑CGTGATGCAGATCACACAATACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                          <  2:319001‑M1/143‑23 (MQ=255)                   tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCACATACAAATTATTGTTTtataagtcatacttttgtcttgg                                                                                                                         >  1:239511‑M1/1‑120 (MQ=255)                   ttgctccgaAATCGACGTGACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAAGTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:924814‑M1/139‑24 (MQ=255)                         cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacctttgtgttgggtgtat                                                                                                                   >  1:609930‑M1/1‑114 (MQ=255)                         cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGACTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtat                                                                                                                   >  2:112204‑M1/1‑114 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGCTTTTCAATACTG‑AGTGATGCAGACCACACAAAACACTCAATTACTTAACATACCTGtatacgtcatactattgctctgggtgtatt                                                                                                                  >  2:302010‑M1/1‑113 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                  >  2:116027‑M1/1‑113 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttagttttggatgtatt                                                                                                                  >  2:22034‑M1/1‑113 (MQ=255)                           aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                 <  2:61579‑M1/143‑32 (MQ=255)                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCGCACAAAACACCCAATTACTTAACATAAATGtataagtcatacgtttgttttgggtgtatttca                                                                                                               >  2:193200‑M1/1‑110 (MQ=255)                               cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATAtataagtcatacttttgtattgggtgtatatcaat                                                                                                             >  2:410136‑M1/1‑108 (MQ=255)                               cGACGCAACGGTACTCGTCTCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             <  2:1272329‑M1/143‑36 (MQ=255)                                gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                            >  2:726242‑M1/1‑107 (MQ=255)                                 aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATCCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcaaacttttgtttttgttgtatataaatct                                                                                                           >  2:33321‑M1/1‑106 (MQ=255)                                       cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     <  1:6077‑M1/143‑44 (MQ=255)                                        tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                    >  1:346078‑M1/1‑99 (MQ=255)                                         gTACTCGTCTCTACACGCACGTACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACAAAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:1261094‑M1/143‑46 (MQ=255)                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttttttgggtgtatttccatctgttaaaaagt                                                                                                  >  1:759304‑M1/1‑97 (MQ=255)                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   >  1:1210928‑M1/1‑97 (MQ=255)                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:1210928‑M1/142‑46 (MQ=255)                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctcttaaaaag                                                                                                  >  2:500654‑M1/1‑97 (MQ=255)                                             tCGTCAGTACACGCACTTACAACGGAGGTGCGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACTCTCAATTACTTAACTTAAACGtataagtcatacttttgttttgggtgtatttcaatctgtgaaaaagttt                                                                                               <  2:371520‑M1/143‑50 (MQ=255)                                             tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgaatttcaatctgttaaaaagttt                                                                                               >  1:448446‑M1/1‑94 (MQ=255)                                               gTCACTACACGCACCTGCAACGGAGGCGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                             <  2:520701‑M1/143‑52 (MQ=255)                                                   cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  1:477862‑M1/1‑88 (MQ=255)                                                         gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   >  2:114463‑M1/1‑82 (MQ=255)                                                           acaTACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttgttcgcttcgctagc                                                                                 >  2:852493‑M1/1‑80 (MQ=255)                                                              tACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                              <  2:247179‑M1/143‑67 (MQ=255)                                                                cAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                            >  1:263929‑M1/1‑75 (MQ=255)                                                                    ggAGGTGCGCTTCGATTTTCAATAATG‑CGTGATGCATATCACATAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  2:520324‑M1/143‑73 (MQ=255)                                                                    ggAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:1047981‑M1/143‑73 (MQ=255)                                                                    ggAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:25332‑M1/143‑73 (MQ=255)                                                                    ggAGGGGGGCTGCACTTTTCAACAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  2:621368‑M1/143‑73 (MQ=255)                                                                     gAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       <  1:138081‑M1/143‑74 (MQ=255)                                                                         ggggCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   <  1:1293057‑M1/143‑78 (MQ=255)                                                                           ggCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                                 <  1:1215019‑M1/143‑80 (MQ=255)                                                                            gCTGCGATTTTCAATAATG‑CGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                <  1:980149‑M1/143‑81 (MQ=255)                                                                             cTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctacgctagcaagccaaaatgaaatagg                                                               >  2:263353‑M1/1‑62 (MQ=255)                                                                                    tttCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaccgtttatgtgcacgcgagcaagcaaaaatgaaacaggaataatc                                                        >  2:271345‑M1/1‑55 (MQ=255)                                                                                      tCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagcttctgtcgacgcgagcaagcaaaaatgaaacaggaataatcga                                                      >  2:272223‑M1/1‑53 (MQ=255)                                                                                          taatG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggactcatcgaaatg                                                  >  1:326789‑M1/1‑49 (MQ=255)                                                                                          taagG‑CGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                  <  1:39227‑M1/143‑95 (MQ=255)                                                                                           tcacgcgtccagcatatcacccaCAACACTCACCTGCTTAACCTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                  <  1:42094‑M1/132‑95 (MQ=255)                                                                                           aatG‑CGTCATGCAGATCACGTAGAACACTCAAGTACTTAACGTAAATGtataagtcatacttttgttttgggtgtgtttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                 <  2:661478‑M1/143‑96 (MQ=255)                                                                                                    aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                         >  1:976131‑M1/1‑40 (MQ=255)                                                                                                      gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaatcatgaaacaggaatcatcgaaatgggatactgcgc                                       >  2:1020145‑M1/1‑38 (MQ=255)                                                                                                      gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  2:882346‑M1/1‑38 (MQ=255)                                                                                                        aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  2:502281‑M1/1‑36 (MQ=255)                                                                                                         gATTACACAAAAGACTCAATTACTTAACATAAATGtataagtcatacttttgttttggggggatttcaatctgttaaaaagtttttcgctacgcgagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:554637‑M1/143‑109 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:706081‑M1/143‑109 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  1:778033‑M1/143‑109 (MQ=255)                                                                                                           tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                  <  1:344713‑M1/143‑111 (MQ=255)                                                                                                            cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtc                                 <  1:922488‑M1/143‑112 (MQ=255)                                                                                                               acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              >  1:1201659‑M1/1‑29 (MQ=255)                                                                                                                 aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaaagtcaaaat                            >  1:976328‑M1/1‑27 (MQ=255)                                                                                                                  aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                           >  1:726642‑M1/1‑26 (MQ=255)                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                   >  2:862537‑M1/1‑18 (MQ=255)                                                                                                                            ttACGTAACATAAATGtataagtcatacttttgttttgggtgtagttcaatctgttaaaaagtgtttcgctacgctagcaagcaaaaatgaaccaagaatgatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                 <  2:1076361‑M1/143‑128 (MQ=255)                                                                                                                             tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaacagttcttcgctacgctagcatgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                >  1:971041‑M1/1‑15 (MQ=255)                                                                                                                             tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtccaaataactcaccgtaa                >  1:763756‑M1/1‑15 (MQ=255)                                                                                                                               cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatgggatgttgcgcacagtcaaaataactcaccataaat              >  1:590472‑M1/1‑13 (MQ=255)                                                                                                                               cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat              >  1:1215701‑M1/1‑13 (MQ=255)                                                                                                                                ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       >  2:176850‑M1/1‑12 (MQ=255)                                                                                                                                ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       <  1:176850‑M1/85‑74 (MQ=255)                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaacaaacaacaagtgagaagcgaaatgggatgttgcgcacagtcaaaataactcaccgtaactaatcat        >  1:609107‑M1/1‑7 (MQ=255)                                                                                                                                        aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     >  2:849203‑M1/1‑4 (MQ=255)                                                                                                                                           gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctatcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:1266816‑M1/1‑1 (MQ=255)                                                                                                                                           gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccctaaacaagcatccgctg  >  1:155765‑M1/1‑1 (MQ=255)                                                                                                                                           gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatcgaaatgggatgttgcggacagacaaaataactcgcggtaagtaatcatatggta  >  1:742535‑M1/1‑1 (MQ=255)                                                                                                                                           |                                                                                                                                               CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG‑CGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644

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