BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,717	1	.	C	100.0%	648.0 / NA	209	T29R (ACG→AGG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (84/125); total (84/125)

Only 100 of 242 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGGA‑C‑g                                                                                         >  1:3959755/1‑93 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gc                                                                                        >  1:2849976/1‑92 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gc                                                                                        <  2:2944162/93‑2 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAATTCAAC‑TGAAGTA‑C‑gctg                                                                                      >  2:497361/1‑90 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gctg                                                                                      <  2:2878965/93‑4 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gctgg                                                                                     <  1:2143136/93‑5 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gctgg                                                                                     <  1:2407690/93‑5 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gctgg                                                                                     <  2:310397/93‑5 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑gctgg                                                                                     <  2:2604520/93‑5 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCta                                                                                >  2:897309/1‑93 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCta                                                                                <  2:2485946/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCtata                                                                              >  2:1007948/1‑93 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCtata                                                                              <  2:3926333/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCtata                                                                              <  1:1662105/93‑1 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAC‑TCAACTTGAAGTA‑C‑GCCGGCGCTATAc                                                                             >  2:703785/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATAc                                                                             <  1:646111/93‑1 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATAc                                                                             >  1:3910680/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATAc                                                                             <  2:3727998/93‑1 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCa                                                                         >  1:1962021/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCa                                                                         <  2:4584640/93‑1 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGTAATAAA‑TCAACTTGAGGTA‑CGG‑AGGCGCTATACGCCaa                                                                        >  1:1453406/1‑93 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATAccccca                                                                        >  1:1179519/1‑91 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaa                                                                        <  1:975403/93‑1 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaa                                                                       <  2:2660442/93‑1 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaa                                                                       >  1:2319633/1‑93 (MQ=255)

                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaata                                                                     <  1:246229/93‑1 (MQ=255)

                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGCCACGTAGTAAA‑TCATCTTGAAGTA‑C‑GCTGGCGCTATACCCCaaata                                                                     <  2:2772157/93‑1 (MQ=255)

                     aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataa                                                                    <  2:1353551/93‑1 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaa                                                                   >  1:3842801/1‑93 (MQ=255)

                       aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaa                                                                  <  2:3423422/93‑1 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaat                                                                 <  2:105077/93‑1 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCaaataaaata                                                                <  1:508519/93‑1 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAA‑TCAACATGAAGTACC‑A‑TGGCGCTATACCCCaaataaaata                                                                >  2:2121312/1‑93 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAaataaaataa                                                               <  2:2647376/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAaataaaataa                                                               <  2:3060635/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAaataaaataa                                                               <  1:1609919/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAaataaaataa                                                               >  1:155581/1‑93 (MQ=255)

                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAg                                                              >  2:2162037/1‑93 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:4435081/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  2:3167337/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:1361598/93‑1 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                            >  1:234728/1‑93 (MQ=255)

                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGc                                                          >  1:457875/1‑93 (MQ=255)

                                   aacACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:4073287/91‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:4330237/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:1602386/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  2:3459161/93‑1 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCGagag                                                     >  2:4407062/1‑93 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAACGGCAagag                                                     >  2:476605/1‑93 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAACGGCAagag                                                     >  2:1658156/1‑93 (MQ=255)

                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                    <  1:2930397/93‑1 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                   >  2:1441835/1‑93 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                               >  1:2964818/1‑93 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                              <  1:2867135/93‑1 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGAGGGTaa                                              >  2:2966302/1‑93 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             <  2:633591/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  2:4064459/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:2131799/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           <  1:1831893/93‑1 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCGCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  1:3621752/1‑93 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  1:192543/1‑93 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAAATTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  1:326413/1‑93 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         <  1:590214/93‑1 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:4518136/93‑1 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  1:1325155/93‑1 (MQ=255)

                                                    gATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                                     >  1:3985329/1‑93 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  >  2:680803/1‑93 (MQ=255)

                                                        tCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc                                 <  2:1823564/93‑1 (MQ=255)

                                                        tCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc                                 <  2:1461420/93‑1 (MQ=255)

                                                            cggcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                             >  1:1100887/1‑93 (MQ=255)

                                                             ggcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            >  1:2027027/1‑93 (MQ=255)

                                                             ggcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            <  2:669666/93‑1 (MQ=255)

                                                              gcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt                           <  2:3808733/93‑1 (MQ=255)

                                                              gcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt                           <  2:1190494/93‑1 (MQ=255)

                                                              gcggGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgt                           <  1:1240952/93‑1 (MQ=255)

                                                               cggGCACGTAGTAAA‑TCAACTGGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAAGAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          <  2:2053287/93‑1 (MQ=255)

                                                                gggCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATGCGCCCTGCCGCgtttgttt                         <  2:2440334/93‑1 (MQ=255)

                                                                  gCACGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTg                       <  1:4502193/93‑1 (MQ=255)

                                                                    aCGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:3569975/1‑93 (MQ=255)

                                                                     cGTAGTAAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                    <  2:3661133/93‑1 (MQ=255)

                                                                      gtagtaAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                   >  1:873360/1‑93 (MQ=255)

                                                                      gtagtaAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                   <  2:784353/93‑1 (MQ=255)

                                                                      gtagtaAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                   >  1:1005115/1‑93 (MQ=255)

                                                                       cagtaAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  >  1:2554332/2‑93 (MQ=255)

                                                                         gtaAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                <  1:4214723/93‑1 (MQ=255)

                                                                          taAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               >  2:563509/1‑93 (MQ=255)

                                                                          taAA‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               <  1:321258/93‑1 (MQ=255)

                                                                            aa‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  1:1786348/1‑93 (MQ=255)

                                                                             a‑TCAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCggg            >  2:290320/1‑93 (MQ=255)

                                                                                cAACTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGaa          <  1:1208748/93‑1 (MQ=255)

                                                                                 aaCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         >  1:3775329/1‑93 (MQ=255)

                                                                                  aCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        <  1:4252858/93‑1 (MQ=255)

                                                                                  aCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        >  2:2643328/1‑93 (MQ=255)

                                                                                  aCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        <  2:3764532/93‑1 (MQ=255)

                                                                                  aCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        <  1:1427013/93‑1 (MQ=255)

                                                                                  aCTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        >  1:1296985/1‑93 (MQ=255)

                                                                                   cTTGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       >  2:922389/1‑93 (MQ=255)

                                                                                     tGAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTc     <  1:3351176/93‑1 (MQ=255)

                                                                                      gAAGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    <  1:1582039/93‑1 (MQ=255)

                                                                                       aaGTA‑C‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt   <  2:849884/93‑1 (MQ=255)

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAA‑TCAACTTGAAGTA‑C‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

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