BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,805	T→C	W222R (TGG→CGG)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,805	0	T	C	100.0%	903.9 / NA	276	W222R (TGG→CGG)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (148/128); total (148/128)

Only 100 of 293 total aligned reads displayed.

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

tATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACtggcgg                                                                                         >  2:4564784/1‑93 (MQ=255)

tATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACtggcgg                                                                                         >  2:3701504/1‑93 (MQ=255)

 aTCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCggg                                                                                        >  2:5985073/1‑93 (MQ=255)

 aTCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCggg                                                                                        >  2:6033442/1‑93 (MQ=255)

  tCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGc                                                                                       <  1:6055621/93‑1 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                  <  2:358902/93‑1 (MQ=255)

       tGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGctacta                                                                                  >  2:1585286/1‑93 (MQ=255)

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCa                                                                              <  2:4406122/93‑1 (MQ=255)

           ggCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCa                                                                              <  1:3253631/93‑1 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                             <  2:3986987/93‑1 (MQ=255)

            gCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAg                                                                             >  2:5545145/1‑93 (MQ=255)

             cGAAAGTCGCAAAACGCAGAGGCCTGCCTTATCCCCACGTCAATCAGCAAGCCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAgg                                                                            <  1:166857/93‑1 (MQ=255)

              gAAAGTCGCAAAGCGCAAAGACCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATCCCGACTGGCGGGCTACTATGCAGGc                                                                           <  1:967246/93‑1 (MQ=255)

              gAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGc                                                                           >  1:3115956/1‑93 (MQ=255)

               aaaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCt                                                                          >  1:3787229/1‑93 (MQ=255)

               aaaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCt                                                                          <  1:698719/93‑1 (MQ=255)

                aaGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTg                                                                         >  1:2411481/1‑93 (MQ=255)

                 aGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGGTACTATGCAGGCTgt                                                                        >  1:1516277/1‑92 (MQ=255)

                 aGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGc                                                                        >  1:2959008/1‑93 (MQ=255)

                   tCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAg                                                                      <  2:763331/93‑1 (MQ=255)

                    cGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAgg                                                                     <  2:1238744/93‑1 (MQ=255)

                     gCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAggg                                                                    >  1:6302273/1‑93 (MQ=255)

                       aaaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCGACTATGCAGGCTGCAGGGTa                                                                  >  1:3834136/1‑93 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTaa                                                                 >  2:182326/1‑93 (MQ=255)

                        aaaCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTaa                                                                 <  1:6459292/93‑1 (MQ=255)

                          aCGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGt                                                               <  1:889859/93‑1 (MQ=255)

                           cGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTg                                                              <  1:3904790/93‑1 (MQ=255)

                            gCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGc                                                             >  1:138058/1‑93 (MQ=255)

                             cAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGcc                                                            <  1:3690228/93‑1 (MQ=255)

                             cAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGcc                                                            <  2:58880/93‑1 (MQ=255)

                              aaaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCCACTATGCAGGCTGCCGGGGAAGGGccc                                                           >  2:4195038/1‑92 (MQ=255)

                               aaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCAt                                                          >  1:1083027/1‑93 (MQ=255)

                               aaGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCAt                                                          <  2:4243104/93‑1 (MQ=255)

                                aGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATg                                                         >  2:1765363/1‑93 (MQ=255)

                                  gcctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTTTGCAGGCTGCAGGGTAAGTGCCATGaa                                                       >  1:610193/1‑93 (MQ=255)

                                   cctgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGaaa                                                      <  1:3053204/93‑1 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGGAGGCTGCAGGGTAAGTGCCATGAAAgg                                                    >  2:5170395/1‑93 (MQ=255)

                                     tgcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAgg                                                    >  1:1259800/1‑93 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGTGTAAGTGCCATGAAAGGt                                                   <  1:3178464/93‑1 (MQ=255)

                                      gcctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGACGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGt                                                   >  2:6158392/1‑93 (MQ=255)

                                       cctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTc                                                  <  1:3989920/93‑1 (MQ=255)

                                        ctTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCg                                                 >  1:4316696/1‑93 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGc                                                >  2:4890443/1‑93 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGc                                                <  2:888919/93‑1 (MQ=255)

                                         ttATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGc                                                >  2:4995770/1‑93 (MQ=255)

                                          tATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGcc                                               >  1:2605758/1‑93 (MQ=255)

                                            tCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTg                                             >  2:3420533/1‑93 (MQ=255)

                                             ccccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGTAAGGTAAGTGCCATGAAAGGTCGCCTGt                                            >  2:6424812/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGtt                                           >  1:4408200/1‑93 (MQ=255)

                                              cccACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGtt                                           <  1:4168655/93‑1 (MQ=255)

                                                 ccGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  1:5943365/92‑1 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        >  1:1353268/1‑93 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGa                                        <  1:1850267/93‑1 (MQ=255)

                                                 aCGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGAGGGCTACTATGCAGGCCGCAGGGGAAATGCCATGAAAGGTCGCCTGTgaga                                        >  2:6592724/1‑93 (MQ=255)

                                                  cGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGAt                                       >  2:635355/1‑93 (MQ=255)

                                                    tCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGc                                     <  2:1300036/93‑1 (MQ=255)

                                                       aTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGt                                  >  1:5095341/1‑93 (MQ=255)

                                                       aTCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGt                                  >  1:5976734/1‑93 (MQ=255)

                                                        tCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTc                                 <  2:3898432/93‑1 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTccc                               >  2:2860951/1‑93 (MQ=255)

                                                          aGCAAGGCGAGATCGAAGCTGATGCAGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTccc                               <  2:6481377/93‑1 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                              >  2:2064113/1‑93 (MQ=255)

                                                           gCAAGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCa                              >  1:4598562/1‑93 (MQ=255)

                                                             aaGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACt                            >  1:4762044/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGTCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  2:2294492/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  2:6699885/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  1:197013/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           <  1:3651476/93‑1 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  1:4385120/1‑93 (MQ=255)

                                                              aGGCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTc                           >  2:2560303/1‑93 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                          >  1:1204437/1‑93 (MQ=255)

                                                               ggCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCa                          <  1:5420458/93‑1 (MQ=255)

                                                                gCGAGATCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAg                         >  1:252661/1‑93 (MQ=255)

                                                                 ggagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGt                        <  2:2211578/92‑1 (MQ=255)

                                                                  gagaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGtt                       >  1:5176954/1‑93 (MQ=255)

                                                                   agaTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTc                      <  2:3311438/93‑1 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    <  1:2195379/93‑1 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    >  1:2758766/1‑93 (MQ=255)

                                                                     aTCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTccc                    <  1:1426794/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   <  1:1301269/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   <  1:2112583/93‑1 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   >  2:1181239/1‑93 (MQ=255)

                                                                      tCGAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCt                   <  2:561370/93‑1 (MQ=255)

                                                                        gAAGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGGAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTaa                 >  1:3446433/1‑93 (MQ=255)

                                                                          aGCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACg               <  2:831198/93‑1 (MQ=255)

                                                                           gCTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACgg              <  1:6587909/93‑1 (MQ=255)

                                                                            cTGATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACggg             >  1:4109416/1‑93 (MQ=255)

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCg           >  2:5064239/1‑93 (MQ=255)

                                                                              gATGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGCAAGGGCCATGAAAGGGCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCg           >  1:3582860/1‑93 (MQ=255)

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGt          >  2:556743/1‑93 (MQ=255)

                                                                               aTGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGt          <  1:819463/93‑1 (MQ=255)

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGtt         <  1:5568222/93‑1 (MQ=255)

                                                                                tGCCGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGtt         >  2:4287429/1‑93 (MQ=255)

                                                                                  ccGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTgg       <  1:4860421/93‑1 (MQ=255)

                                                                                   cGACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgc      <  2:338174/93‑1 (MQ=255)

                                                                                    gACTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcg     >  2:2002910/1‑93 (MQ=255)

                                                                                     aCTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGgcgc    <  2:2398332/93‑1 (MQ=255)

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca   >  1:6534238/1‑93 (MQ=255)

                                                                                      cTGGCGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgca   <  2:1138810/93‑1 (MQ=255)

                                                                                       tggcggGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCgcag  >  1:6199717/1‑93 (MQ=255)

TATCCGGTGCTGGCGAAAGTCGCAAAACGCAAAGGCCTGCCTTATCCCCACGTCAATCAGCAAGGCGAGATCGAAGCTGATGCCGACTGGTGGGCTACTATGCAGGCTGCAGGGTAAGTGCCATGAAAGGTCGCCTGTTAGATGCTGTCCCACTCAGTTCCCTAACGGGCGTTGGCGCAG  >  NC_000913/3822715‑3822894

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