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–1979270	1979278	9–776	96 [0]	[4] 97	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)	85 (1.070)	61/182	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

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601                                                                                                    |                                                                                                    cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  1:530197‑M1/100‑1 (MQ=255) cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  2:1176602‑M1/100‑1 (MQ=255)  tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                    <  1:106275‑M1/100‑2 (MQ=255)  tGTACTCGTCACTACACGCACAGACAACGGAGGGGAGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  2:719489‑M1/99‑1 (MQ=255)   gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                   <  2:2241034‑M1/100‑3 (MQ=255)      cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                >  1:115087‑M1/1‑95 (MQ=255)       tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                               <  2:2227082‑M1/100‑7 (MQ=255)        gatctctacacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              <  1:1587080‑M1/95‑8 (MQ=255)        cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  2:1651879‑M1/1‑93 (MQ=255)        cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  2:1653683‑M1/1‑93 (MQ=255)        cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  1:1676048‑M1/1‑93 (MQ=255)         gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  2:2063432‑M1/100‑9 (MQ=255)         gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  1:1687903‑M1/100‑9 (MQ=255)          tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                            >  1:2305330‑M1/1‑91 (MQ=255)          tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                            >  1:885760‑M1/1‑91 (MQ=255)           cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           <  1:2235908‑M1/100‑11 (MQ=255)           cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           <  2:1675762‑M1/100‑11 (MQ=255)            aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                          >  1:1158655‑M1/1‑89 (MQ=255)            aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                          >  2:1140846‑M1/1‑89 (MQ=255)             cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtagatc                                                                                         >  2:1587076‑M1/1‑88 (MQ=255)               acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACGAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                       >  2:1123682‑M1/1‑86 (MQ=255)                cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                      >  2:1190118‑M1/1‑85 (MQ=255)                  cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    >  1:871100‑M1/1‑83 (MQ=255)                  cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    >  1:800493‑M1/1‑83 (MQ=255)                  cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  1:1884356‑M1/100‑18 (MQ=255)                    cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                  <  2:1158652‑M1/100‑20 (MQ=255)                     acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                 >  1:620409‑M1/1‑80 (MQ=255)                       aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                               <  2:797259‑M1/100‑23 (MQ=255)                        tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:739931‑M1/1‑77 (MQ=255)                        tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  2:2251361‑M1/1‑77 (MQ=255)                        tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  2:1628145‑M1/1‑77 (MQ=255)                        tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              <  1:325718‑M1/100‑24 (MQ=255)                         aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             >  1:2166459‑M1/1‑76 (MQ=255)                         aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:2251372‑M1/100‑25 (MQ=255)                         aCAACGGAGGGGGGCTGCGACTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             >  2:1407832‑M1/1‑76 (MQ=255)                           aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                           <  1:688178‑M1/100‑27 (MQ=255)                            aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                          >  1:1074441‑M1/1‑73 (MQ=255)                             cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                         >  1:1761135‑M1/1‑72 (MQ=255)                               gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                       >  1:688848‑M1/1‑70 (MQ=255)                                 ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  2:1989151‑M1/100‑33 (MQ=255)                                 atctggctgcGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  1:851117‑M1/96‑33 (MQ=255)                                  gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                    <  2:921581‑M1/100‑34 (MQ=255)                                     ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcagat                                                                 >  2:851115‑M1/1‑64 (MQ=255)                                      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                <  2:721625‑M1/100‑38 (MQ=255)                                       cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               >  1:1022776‑M1/1‑62 (MQ=255)                                        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              <  1:2236055‑M1/100‑40 (MQ=255)                                        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              >  1:911687‑M1/1‑61 (MQ=255)                                        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              <  2:871098‑M1/100‑40 (MQ=255)                                         gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  1:738542‑M1/1‑60 (MQ=255)                                         gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  2:479862‑M1/1‑60 (MQ=255)                                         gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  2:920548‑M1/1‑60 (MQ=255)                                           gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                           <  1:553576‑M1/100‑43 (MQ=255)                                               ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       <  1:1651884‑M1/100‑47 (MQ=255)                                               ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       <  1:1653688‑M1/100‑47 (MQ=255)                                               ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       <  1:1803075‑M1/100‑47 (MQ=255)                                                 tccgatctcgtgatgcaGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                     <  1:1231259‑M1/92‑49 (MQ=255)                                                 cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                     <  2:911685‑M1/100‑49 (MQ=255)                                                 cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                     >  1:158355‑M1/1‑52 (MQ=255)                                                  aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                    >  2:1019465‑M1/1‑51 (MQ=255)                                                  aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                    <  1:1315794‑M1/100‑50 (MQ=255)                                                     aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                 <  2:2256060‑M1/100‑53 (MQ=255)                                                      tccgatcttgcagatcaCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                <  2:688570‑M1/92‑54 (MQ=255)                                                      atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                <  2:739929‑M1/100‑54 (MQ=255)                                                       tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               <  2:1432413‑M1/100‑55 (MQ=255)                                                         cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                             <  1:73787‑M1/100‑57 (MQ=255)                                                         cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttagatcgga                                             >  2:1231256‑M1/1‑44 (MQ=255)                                                            gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                          >  2:1937948‑M1/1‑41 (MQ=255)                                                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgagatcgga                                        >  1:688572‑M1/1‑39 (MQ=255)                                                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                       <  1:1140848‑M1/100‑63 (MQ=255)                                                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                       >  2:726983‑M1/1‑38 (MQ=255)                                                                cAGATAACACAAAACAATCAATTACTTAACATAGATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                      <  2:2166449‑M1/100‑64 (MQ=255)                                                                   aTCATACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                   <  2:1811511‑M1/100‑67 (MQ=255)                                                                     cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                 >  2:1461517‑M1/1‑32 (MQ=255)                                                                      acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                >  1:1472162‑M1/1‑31 (MQ=255)                                                                       cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               >  1:511457‑M1/1‑30 (MQ=255)                                                                         caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatgtgttaaaaagtttttcgctacgctagcaagcaaaaa                             <  1:1628150‑M1/100‑73 (MQ=255)                                                                          aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                            <  2:100644‑M1/100‑74 (MQ=255)                                                                           aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                            <  1:1123684‑M1/99‑74 (MQ=255)                                                                              cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                        >  1:178950‑M1/1‑23 (MQ=255)                                                                                cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                      <  1:2162527‑M1/100‑80 (MQ=255)                                                                                    aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgagacaggaa                  >  2:1779357‑M1/1‑17 (MQ=255)                                                                                     ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                 >  2:417133‑M1/1‑16 (MQ=255)                                                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa               >  2:1432691‑M1/1‑14 (MQ=255)                                                                                        cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagcttagctagcaagcaaaaatgaaaaaggaataa               >  1:2076901‑M1/1‑13 (MQ=255)                                                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc             <  2:1472158‑M1/100‑89 (MQ=255)                                                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgaaagaaaaaaaaaatgaaacaggaataatc             <  2:885758‑M1/100‑89 (MQ=255)                                                                                          tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            >  1:1351433‑M1/1‑11 (MQ=255)                                                                                           aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga           >  2:1121806‑M1/1‑10 (MQ=255)                                                                                             cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  2:2108468‑M1/1‑8 (MQ=255)                                                                                             cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  1:1046306‑M1/1‑8 (MQ=255)                                                                                               tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       <  2:178951‑M1/100‑95 (MQ=255)                                                                                                naatGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg      >  2:1100652‑M1/2‑5 (MQ=255)                                                                                                aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg      >  1:399596‑M1/1‑5 (MQ=255)                                                                                                 aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagatagaaagcaaaaatgaaacaggaataatcgaaatggg     <  2:115088‑M1/100‑97 (MQ=255)                                                                                                 aaTGtataagtcatacttttgttttgggtgtatgtcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  2:1520933‑M1/1‑4 (MQ=255)                                                                                                    gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  <  2:158356‑M1/100‑100 (MQ=255)                                                                                                    |                                                                                                    CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

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