BRESEQ :: Evidence

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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	125 [1]	[0] 124	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)	103 (1.160)	71/162	NT	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
*	?	NC_000913	1979279 =	0 (0.000)	103 (1.160)	71/162	NT	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 125 total aligned reads displayed.

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978406‑1978601

 aCTCGTCATTACACGCACATACAACGTAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                 <  1:188953‑M1/100‑5 (MQ=255)

   tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                               <  2:222952‑M1/100‑7 (MQ=255)

    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  2:2612208‑M1/1‑93 (MQ=255)

     gTCACTACACGCACATACAACGGTGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  1:1680844‑M1/100‑9 (MQ=255)

     gTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  1:803241‑M1/100‑9 (MQ=255)

         ntaCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                         >  2:119135‑M1/2‑88 (MQ=255)

          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                        <  1:2445319‑M1/100‑14 (MQ=255)

          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                        >  1:1817627‑M1/1‑87 (MQ=255)

              cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    >  1:317082‑M1/1‑83 (MQ=255)

              cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  2:59704‑M1/100‑18 (MQ=255)

              cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  1:1489270‑M1/100‑18 (MQ=255)

                cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                  >  2:2003690‑M1/1‑81 (MQ=255)

                cacaTACAACGGAGGGGCGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCCATTACTTACCATAAATGtataagtcatacttttgtt                                                                                  <  2:252958‑M1/100‑20 (MQ=255)

                  caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                <  1:859162‑M1/100‑22 (MQ=255)

                    tACAACGGAGGGGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              <  1:1096269‑M1/100‑24 (MQ=255)

                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:482252‑M1/1‑77 (MQ=255)

                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:152105‑M1/1‑77 (MQ=255)

                     aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:1929500‑M1/100‑25 (MQ=255)

                      cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  2:234209‑M1/1‑75 (MQ=255)

                      cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            <  2:695692‑M1/100‑26 (MQ=255)

                      cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  2:527821‑M1/1‑75 (MQ=255)

                      cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  1:2369005‑M1/1‑75 (MQ=255)

                      cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             >  1:987737‑M1/1‑75 (MQ=255)

                       aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcatacttttgttttgggtg                                                                           <  2:181890‑M1/100‑27 (MQ=255)

                       aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                           >  2:1650659‑M1/1‑74 (MQ=255)

                           gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                       <  2:782008‑M1/100‑31 (MQ=255)

                             ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  2:953598‑M1/100‑33 (MQ=255)

                             ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     >  2:2547442‑M1/1‑68 (MQ=255)

                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                  >  1:857999‑M1/1‑65 (MQ=255)

                                 ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  2:1476936‑M1/1‑64 (MQ=255)

                                 ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  1:2097532‑M1/1‑64 (MQ=255)

                                 ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 <  1:1515619‑M1/100‑37 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                <  2:152104‑M1/100‑38 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                >  2:620171‑M1/1‑63 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                <  2:495475‑M1/100‑38 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 <  1:2152682‑M1/99‑37 (MQ=255)

                                   cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               >  2:345046‑M1/1‑62 (MQ=255)

                                    tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              <  2:1213976‑M1/100‑40 (MQ=255)

                                     gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgtt                                                             >  1:355168‑M1/1‑60 (MQ=255)

                                     gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             <  1:1356451‑M1/100‑41 (MQ=255)

                                     gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             <  1:1650675‑M1/100‑41 (MQ=255)

                                     gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             >  2:1518122‑M1/1‑60 (MQ=255)

                                     gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                             <  1:2556645‑M1/100‑41 (MQ=255)

                                      cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            <  1:620175‑M1/100‑42 (MQ=255)

                                       gATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                           <  1:1068206‑M1/100‑43 (MQ=255)

                                        aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                          <  1:1991663‑M1/100‑44 (MQ=255)

                                        aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                           <  2:1827831‑M1/99‑43 (MQ=255)

                                         ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                         <  1:1804907‑M1/100‑45 (MQ=255)

                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       <  1:2025642‑M1/100‑47 (MQ=255)

                                              aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                    >  2:2425254‑M1/1‑51 (MQ=255)

                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactattgttataggtatatttaaatatgttaaaaagat                                                     >  1:145678‑M1/1‑50 (MQ=255)

                                                 aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATAtataaaacatactattgatttgggtgtatttcaatctgttaaaaagtttttc                                                 >  1:116869‑M1/1‑48 (MQ=255)

                                                   tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               >  1:1582510‑M1/1‑46 (MQ=255)

                                                   tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttaggtgtatttcaatctgttaaaaagtttttcgc                                               >  2:2464814‑M1/1‑46 (MQ=255)

                                                     cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                             <  2:948183‑M1/100‑57 (MQ=255)

                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattttaatctgttaaaaagtttttagata                                             >  2:108159‑M1/1‑43 (MQ=255)

                                                       tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                           >  1:222726‑M1/1‑42 (MQ=255)

                                                         aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  2:1656667‑M1/100‑61 (MQ=255)

                                                         aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttaggtgtatttcaatctgttaaaaagtttttcgctacgct                                         >  1:1044414‑M1/1‑40 (MQ=255)

                                                          tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  2:1457933‑M1/1‑39 (MQ=255)

                                                             aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     <  1:330799‑M1/100‑65 (MQ=255)

                                                             aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     >  2:884515‑M1/1‑36 (MQ=255)

                                                             aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                      <  2:279899‑M1/99‑64 (MQ=255)

                                                                tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                  >  2:591713‑M1/1‑33 (MQ=255)

                                                                 cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                 >  1:2500627‑M1/1‑32 (MQ=255)

                                                                  acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttaggtgtatttcaatatgttaaaaagtttttagctacgctagcaagcaa                                >  1:1434495‑M1/1‑31 (MQ=255)

                                                                   cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               >  1:187290‑M1/1‑30 (MQ=255)

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             >  2:2008864‑M1/1‑28 (MQ=255)

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                             >  2:1252852‑M1/1‑28 (MQ=255)

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctaataaaaagtttttcgctacgctagcaagcaaaaa                             >  1:221653‑M1/1‑28 (MQ=255)

                                                                          cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                        >  2:902796‑M1/1‑23 (MQ=255)

                                                                          cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                        <  2:1316508‑M1/100‑78 (MQ=255)

                                                                            cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  2:607138‑M1/99‑79 (MQ=255)

                                                                              cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttttgctacgctagcaagcaaaaatgaaacag                     >  2:2214556‑M1/1‑19 (MQ=255)

                                                                              cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  2:516141‑M1/100‑82 (MQ=255)

                                                                              cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  2:176541‑M1/100‑82 (MQ=255)

                                                                              cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  2:1008017‑M1/100‑82 (MQ=255)

                                                                              cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                      >  2:1013916‑M1/1‑19 (MQ=255)

                                                                               aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                   >  1:476971‑M1/1‑18 (MQ=255)

                                                                               aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                   >  1:518779‑M1/1‑18 (MQ=255)

                                                                                aTTACTTAATATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                  >  2:767266‑M1/1‑17 (MQ=255)

                                                                                   aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa               >  1:1086176‑M1/1‑14 (MQ=255)

                                                                                    nttaACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttattcgctacgatagcaagcaaaaatgaaacaggaataat              >  2:1839843‑M1/2‑13 (MQ=255)

                                                                                    cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat              >  2:1505375‑M1/1‑13 (MQ=255)

                                                                                       aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttcaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga           >  2:342341‑M1/1‑10 (MQ=255)

                                                                                        aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacaataaaaaaaaaaaatgaaacaggaataatcgaa          <  2:1511709‑M1/100‑92 (MQ=255)

                                                                                          cataatgtataagtcatacttttggtttgggtgtatttcaatttgttaaaaagttttttgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  2:673561‑M1/4‑7 (MQ=255)

                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       >  2:2010062‑M1/1‑6 (MQ=255)

                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       <  1:2203131‑M1/100‑95 (MQ=255)

                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       >  2:263729‑M1/1‑6 (MQ=255)

                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacactagcaagcaaaaatgaaacaggaataatcgaaat        >  1:797542‑M1/1‑6 (MQ=255)

                                                                                            aaaTGtataagtcatactttcgttttgagtgtatttcaatctgataaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgcaa         >  2:1813658‑M1/1‑5 (MQ=255)

                                                                                             tgtGGGTTATGATTTCAACTTATTGATAGTGTTTTAGGTTCGTATAATGCCCGACGACTGTGTCATGCAGCTCCACCGAGTGTGAGAACGACa            <  2:2000172/91‑1 (MQ=11)

                                                                                             aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  1:1179843‑M1/1‑4 (MQ=255)

                                                                                             aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     >  2:147270‑M1/1‑4 (MQ=255)

                                                                                             aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctaacaagcaaaaatgaaacaggaataatcgaaat        >  1:2214309‑M1/1‑4 (MQ=255)

                                                                                              aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga    >  2:1154288‑M1/1‑3 (MQ=255)

                                                                                               tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat   <  2:375623‑M1/100‑99 (MQ=255)

                                                                                                gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  >  1:2319650‑M1/1‑1 (MQ=255)

                                                                                                gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  >  1:1714035‑M1/1‑1 (MQ=255)

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978406‑1978601

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