BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,178,668	A→G	intergenic (+165/‑155)	rplL → / → rpoB	50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,178,668	0	A	G	100.0%	849.5 / NA	264	intergenic (+165/‑155)	rplL/rpoB	50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (121/143); total (121/143)

Only 100 of 267 total aligned reads displayed.

TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAAT  >  NC_000913/4178578‑4178757

 ttAAGTCGCCATTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATTCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                        <  2:5289049/92‑1 (MQ=255)

 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                        <  1:1781314/93‑1 (MQ=255)

  tAAGCCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTc                                                                                       <  2:3817391/93‑1 (MQ=255)

   aaGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTcc                                                                                      <  2:2983013/93‑1 (MQ=255)

     ggCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGt                                                                                    <  2:1897871/93‑1 (MQ=255)

     ggCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGt                                                                                    >  2:2914232/1‑93 (MQ=255)

      gcgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGtt                                                                                   >  1:3303197/1‑93 (MQ=255)

       cgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTc                                                                                  <  1:602614/93‑1 (MQ=255)

       cgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTc                                                                                  <  2:147734/93‑1 (MQ=255)

        gcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCt                                                                                 <  1:4237194/93‑1 (MQ=255)

         ccAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtg                                                                                <  1:5166777/93‑1 (MQ=255)

          cAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgt                                                                               >  1:4035387/1‑93 (MQ=255)

          cAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgt                                                                               <  2:610865/93‑1 (MQ=255)

            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             >  1:4295727/1‑93 (MQ=255)

            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             <  2:2263547/93‑1 (MQ=255)

            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             >  2:1420708/1‑93 (MQ=255)

              aGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaa                                                                           >  1:145981/1‑93 (MQ=255)

               gCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaaa                                                                          <  1:5607013/93‑1 (MQ=255)

                cGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAAt                                                                         <  1:1030683/93‑1 (MQ=255)

                cGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAAt                                                                         <  1:2843810/93‑1 (MQ=255)

                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAAta                                                                        >  1:3398914/1‑92 (MQ=255)

                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                        >  2:3528686/1‑93 (MQ=255)

                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCg                                                                       >  2:1566207/1‑93 (MQ=255)

                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCg                                                                       >  2:5613829/1‑93 (MQ=255)

                   ttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGc                                                                      >  1:2015114/1‑93 (MQ=255)

                     cacaCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCaa                                                                    >  1:2598755/1‑93 (MQ=255)

                     cacaCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCaa                                                                    <  1:6098888/93‑1 (MQ=255)

                      acacTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAAt                                                                   >  1:1603853/1‑93 (MQ=255)

                      acacTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAAt                                                                   >  1:1268122/1‑93 (MQ=255)

                        acTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGa                                                                 >  2:205471/1‑93 (MQ=255)

                         cTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGaa                                                                >  2:2100240/1‑93 (MQ=255)

                          tGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGCCAGGGCGTCCGTTCTGTGTAAATCGCAATGaaa                                                               <  1:2707913/93‑1 (MQ=255)

                           gTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGaaaa                                                              >  2:4633016/1‑92 (MQ=255)

                             ttGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATgg                                                            <  1:5093092/93‑1 (MQ=255)

                               gACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGtt                                                          >  1:5368024/1‑93 (MQ=255)

                                actactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGttt                                                         >  1:1894540/1‑93 (MQ=255)

                                actactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGttt                                                         >  2:5726476/1‑93 (MQ=255)

                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                       >  2:3583478/1‑93 (MQ=255)

                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                       >  2:4035552/1‑93 (MQ=255)

                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGCGTAAATCGCAATGAAATGGTTTaa                                                       <  2:5305520/93‑1 (MQ=255)

                                   actGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAg                                                      >  2:737562/1‑93 (MQ=255)

                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                     <  2:4288311/93‑1 (MQ=255)

                                     tgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCg                                                    >  2:5382936/1‑93 (MQ=255)

                                     tgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCg                                                    >  1:1365360/1‑93 (MQ=255)

                                     tgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCg                                                    >  2:1698690/1‑93 (MQ=255)

                                     tgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGGTTAAGCg                                                    >  2:457280/1‑93 (MQ=255)

                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                 >  1:197439/1‑93 (MQ=255)

                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                 >  1:3124510/1‑93 (MQ=255)

                                         gtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGAt                                                <  2:5507073/93‑1 (MQ=255)

                                          tgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATa                                               <  2:5867479/93‑1 (MQ=255)

                                           gCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAg                                              <  2:3135962/93‑1 (MQ=255)

                                             cTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCa                                            <  1:978379/93‑1 (MQ=255)

                                              tttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                                           >  2:1734367/1‑93 (MQ=255)

                                              tttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                                           >  1:5398739/1‑93 (MQ=255)

                                               ttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAAc                                          <  2:1316313/93‑1 (MQ=255)

                                               ttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAAc                                          <  1:6150226/93‑1 (MQ=255)

                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                         <  2:4441066/93‑1 (MQ=255)

                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                         >  1:4446581/1‑93 (MQ=255)

                                                tCAATGATTGTTTGTATCGACGACTTAACATACTGCGACAGGGCGTCCGTGCTGTGTAAAGCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                         <  2:2154562/93‑1 (MQ=255)

                                                 cAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAg                                        <  2:5945502/93‑1 (MQ=255)

                                                  aaTGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAgg                                       >  1:3177760/1‑93 (MQ=255)

                                                  aaTGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAgg                                       <  1:5981191/93‑1 (MQ=255)

                                                   aTGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGc                                      <  1:2131227/93‑1 (MQ=255)

                                                    ttcTTGTTGCTATCGACGACTTAATATACTGCTACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCa                                     <  2:1977965/91‑1 (MQ=255)

                                                      cTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAtt                                   <  2:2675430/93‑1 (MQ=255)

                                                      cTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAtt                                   >  2:4753814/1‑93 (MQ=255)

                                                      cTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAtt                                   >  2:3460203/1‑93 (MQ=255)

                                                          tttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCgg                               <  1:3000332/93‑1 (MQ=255)

                                                            tCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaa                             <  1:2216045/93‑1 (MQ=255)

                                                            tCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGGAAATCGCAACGAAATGGTTGAAGCGTGAGAGCAACAGGCATTGCGGaa                             >  1:6304806/1‑93 (MQ=255)

                                                             cTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaaa                            <  2:565699/93‑1 (MQ=255)

                                                               aTCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgt                          <  2:88612/93‑1 (MQ=255)

                                                                tCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgtg                         >  1:4579882/1‑93 (MQ=255)

                                                                    cgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTcc                     <  2:3250645/93‑1 (MQ=255)

                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTATAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                   <  1:3739891/93‑1 (MQ=255)

                                                                       cTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtt                  <  2:6290785/93‑1 (MQ=255)

                                                                       cTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtt                  <  1:2931707/93‑1 (MQ=255)

                                                                       cTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtt                  <  1:3666429/93‑1 (MQ=255)

                                                                        ttAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAttt                 <  2:4050380/93‑1 (MQ=255)

                                                                        ttAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAttt                 >  1:4965205/1‑93 (MQ=255)

                                                                         tAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtttt                >  1:172843/1‑93 (MQ=255)

                                                                          aaTATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTc               <  2:245005/93‑1 (MQ=255)

                                                                          aaTATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTc               <  2:43907/93‑1 (MQ=255)

                                                                           atatACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTcc              <  1:5788420/93‑1 (MQ=255)

                                                                            tataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCg             >  1:5815672/1‑93 (MQ=255)

                                                                            tataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCg             >  2:4816294/1‑93 (MQ=255)

                                                                             ataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCgg            <  2:2513370/93‑1 (MQ=255)

                                                                             ataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCgg            <  2:3389310/93‑1 (MQ=255)

                                                                             ataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCgg            >  1:4137257/1‑93 (MQ=255)

                                                                              taCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGt           >  1:1048603/1‑93 (MQ=255)

                                                                              taCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGt           <  2:4193876/93‑1 (MQ=255)

                                                                               aCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTc          >  1:937093/1‑93 (MQ=255)

                                                                               aCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTc          >  2:106084/1‑93 (MQ=255)

                                                                                 tGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaa        <  1:2441802/93‑1 (MQ=255)

                                                                                   cGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaaca      <  2:2028752/93‑1 (MQ=255)

                                                                                    gACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaacaa     >  1:5935243/1‑93 (MQ=255)

                                                                                     aCAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaa    >  1:4689995/1‑93 (MQ=255)

                                                                                     aCAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaa    >  1:2410919/1‑93 (MQ=255)

                                                                                      cAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaaa   <  2:5065009/93‑1 (MQ=255)

                                                                                      cAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaaa   <  2:4988430/93‑1 (MQ=255)

TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAAT  >  NC_000913/4178578‑4178757

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