BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,600,078	C→T	W135* (TGG→TAG)	yjjP ←	putative structural protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,600,078	0	C	T	100.0%	821.0 / NA	230	W135* (TGG→TAG)	yjjP	putative structural protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (133/97); total (133/97)

Only 100 of 234 total aligned reads displayed.

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

gCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTAt                                                                                         >  1:1331542/1‑93 (MQ=255)

   aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCtt                                                                                      >  2:206560/1‑93 (MQ=255)

   aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCtt                                                                                      >  1:4573764/1‑93 (MQ=255)

    agaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTg                                                                                     <  2:2402181/93‑1 (MQ=255)

    agaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTg                                                                                     <  1:1718492/93‑1 (MQ=255)

    agaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTg                                                                                     >  2:3920951/1‑93 (MQ=255)

     gaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTgg                                                                                    <  1:510966/93‑1 (MQ=255)

         gTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTaa                                                                                <  2:3741947/93‑1 (MQ=255)

           gatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACg                                                                              >  1:4099609/1‑93 (MQ=255)

           gatgaCGGCACCATCCCAGCCACCGTTATTGAG‑TTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGt                                                                             <  2:3095191/93‑1 (MQ=255)

              gaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTaa                                                                           <  2:1782674/93‑1 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          <  2:3572532/93‑1 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          >  1:4232396/1‑93 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          <  1:3040708/93‑1 (MQ=255)

                cGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATg                                                                         <  2:4419276/93‑1 (MQ=255)

                 ggCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATgg                                                                        >  2:2086415/1‑93 (MQ=255)

                 ggCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATgg                                                                        >  2:4488754/1‑93 (MQ=255)

                  gCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGc                                                                       >  1:4394371/1‑93 (MQ=255)

                   cACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCt                                                                      <  1:1373023/93‑1 (MQ=255)

                   cACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCt                                                                      <  1:343769/93‑1 (MQ=255)

                     ccATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATAGCTGa                                                                    >  2:1960570/1‑93 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTGTCTTGGGTAACGTAATGGCTGaa                                                                   >  2:2588641/1‑93 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                   <  2:1028756/93‑1 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 <  2:2055716/93‑1 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 >  1:1803970/1‑93 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 <  2:621900/93‑1 (MQ=255)

                         cccATCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                <  2:4689740/93‑1 (MQ=255)

                         cccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                <  2:3364082/93‑1 (MQ=255)

                         cccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                <  2:1618706/93‑1 (MQ=255)

                            aGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGc                                                             >  2:4443273/1‑93 (MQ=255)

                             gCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAACGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCt                                                            >  2:1277922/1‑93 (MQ=255)

                              ccACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTa                                                           >  2:3207951/1‑93 (MQ=255)

                              ccACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTa                                                           <  1:2943617/93‑1 (MQ=255)

                               cACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaa                                                          >  2:1203998/1‑93 (MQ=255)

                                aCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGGAACGTAATGGCTGAATTTGGCTaaa                                                         >  1:2011919/1‑93 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        <  1:4469630/93‑1 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        >  1:1101907/1‑93 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        >  1:3961875/1‑93 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCCAAAt                                                        <  2:3028928/93‑1 (MQ=255)

                                  cGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATc                                                       >  2:3416519/1‑93 (MQ=255)

                                   gTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGGCAACCATTAAGGCTACCAGGTAGCTTGAGTAACGAAATGGATGAATTTGGCTAAATCg                                                      >  1:2889514/1‑93 (MQ=255)

                                   gTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCg                                                      >  2:3479335/1‑93 (MQ=255)

                                    ttATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGt                                                     <  1:2774314/93‑1 (MQ=255)

                                    ttATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGt                                                     >  2:4515019/1‑93 (MQ=255)

                                         gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTTTCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                >  1:4129287/1‑93 (MQ=255)

                                         gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                <  1:462482/93‑1 (MQ=255)

                                         gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                >  1:458197/1‑93 (MQ=255)

                                         gAGTTTAAAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                >  2:2883187/1‑93 (MQ=255)

                                          aGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTctct                                               >  2:1315485/1‑93 (MQ=255)

                                           gTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTa                                              <  1:374974/93‑1 (MQ=255)

                                           gTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTa                                              <  1:2495442/93‑1 (MQ=255)

                                           gTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTa                                              >  1:1842785/1‑93 (MQ=255)

                                            tttACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTAc                                             <  2:1759339/93‑1 (MQ=255)

                                               tcAGAAACAGGCCCAAGAAAGGCCAACCATTAAGGCTCCCCGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATTGTTTCTCTACgac                                          >  2:3984490/2‑91 (MQ=255)

                                                cAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCt                                         <  1:4302478/93‑1 (MQ=255)

                                                  gAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCttt                                       >  1:3779773/1‑93 (MQ=255)

                                                   aaaCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTg                                      >  2:3953286/1‑93 (MQ=255)

                                                    aaCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGt                                     >  1:2075511/1‑93 (MQ=255)

                                                     aCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTa                                    >  2:4292289/1‑93 (MQ=255)

                                                     aCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTa                                    >  2:4735370/1‑93 (MQ=255)

                                                       aGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAt                                  >  1:4039146/1‑93 (MQ=255)

                                                       aGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAt                                  <  2:1109144/93‑1 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 >  1:2736917/1‑93 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 >  1:2646251/1‑93 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 >  1:3219475/1‑93 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 >  1:3829770/1‑93 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTCTCTCTACGTCTTTGTAATc                                 >  2:4207161/1‑93 (MQ=255)

                                                          cgcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCca                               <  2:756544/93‑1 (MQ=255)

                                                          cgcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCca                               >  1:3110425/1‑93 (MQ=255)

                                                            cAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagc                             >  2:1659885/1‑93 (MQ=255)

                                                             aaGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATGGTTTCTCTACGCCTTTGTAATCcagca                            >  1:162316/1‑93 (MQ=255)

                                                             aaGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTCAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATgcagca                            >  1:162318/1‑93 (MQ=255)

                                                              aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagcag                           >  2:2360442/1‑93 (MQ=255)

                                                              aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagcag                           <  2:1240577/93‑1 (MQ=255)

                                                               gAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAga                          <  2:592014/93‑1 (MQ=255)

                                                                 aaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatg                        <  2:3350248/93‑1 (MQ=255)

                                                                  aGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatga                       >  1:1518828/1‑93 (MQ=255)

                                                                  aGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatga                       >  1:2344714/1‑93 (MQ=255)

                                                                   ggCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatgat                      >  1:1413718/1‑93 (MQ=255)

                                                                     ccAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGc                    <  2:4004135/93‑1 (MQ=255)

                                                                      cAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCGTTTGTAATCCAGCAGATGATGCt                   >  2:3227636/1‑93 (MQ=255)

                                                                      cAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCt                   >  2:1408476/1‑93 (MQ=255)

                                                                       aaCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTc                  >  1:90117/1‑93 (MQ=255)

                                                                        aCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTcc                 >  2:4566721/1‑93 (MQ=255)

                                                                         ccATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCg                >  1:1889713/1‑93 (MQ=255)

                                                                           aTTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa              >  2:1304316/1‑93 (MQ=255)

                                                                           aTTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa              >  1:1646143/1‑93 (MQ=255)

                                                                            ttAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCaa             <  1:4261835/93‑1 (MQ=255)

                                                                            ttAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCaa             >  1:1744998/1‑93 (MQ=255)

                                                                              aaGGCTACCAGCTATCTTGGTTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGa           >  1:2408468/1‑93 (MQ=255)

                                                                              aaGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGa           >  2:1921002/1‑93 (MQ=255)

                                                                                  cTACCAGCTATCTTGGGTAATGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATca       <  2:1465339/93‑1 (MQ=255)

                                                                                   tACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcac      >  1:4023580/1‑93 (MQ=255)

                                                                                   tACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcac      >  2:519875/1‑93 (MQ=255)

                                                                                     ccAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACaa    <  2:3016544/93‑1 (MQ=255)

                                                                                      cAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAt   >  1:4638909/1‑93 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  >  1:2988575/1‑93 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  >  1:3523112/1‑93 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  >  1:81325/1‑93 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  <  1:1297939/93‑1 (MQ=255)

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

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