BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,091,351	C→T	G168D (GGC→GAC)	rhaD ←	rhamnulose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,091,351	0	C	T	100.0%	927.1 / NA	248	G168D (GGC→GAC)	rhaD	rhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (117/131); total (117/131)

Only 100 of 258 total aligned reads displayed.

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

  ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTc                                                                                       >  1:6140202/1‑93 (MQ=255)

    tCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATAGCACCAGACGCTCGGTGCTGTCTTccc                                                                                     >  2:4245019/1‑93 (MQ=255)

     cGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTcccc                                                                                    >  2:1376801/1‑92 (MQ=255)

        cTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAg                                                                                 >  2:4792874/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGt                                                                                >  1:5528784/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGt                                                                                <  1:2002638/93‑1 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:1295005/1‑93 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               >  2:6153345/1‑93 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               >  1:65722/1‑93 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               >  2:245981/1‑93 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               <  1:771836/93‑1 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGtt                                                                               >  2:1263240/1‑93 (MQ=255)

           gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTg                                                                              <  2:1051905/93‑1 (MQ=255)

           gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTg                                                                              <  1:3569799/93‑1 (MQ=255)

           gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTg                                                                              <  2:3430202/93‑1 (MQ=255)

             cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGc                                                                            >  2:5796837/1‑93 (MQ=255)

              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           >  1:4170818/1‑93 (MQ=255)

               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          <  1:1009797/93‑1 (MQ=255)

               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          >  1:6103658/1‑93 (MQ=255)

               aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                          >  2:581754/1‑93 (MQ=255)

                tttCGTCCGTGCCGGGCACCATCGACGGCAAAATGCCAACGCCATCGGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                                         >  1:5371928/1‑93 (MQ=255)

                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  2:6535004/93‑1 (MQ=255)

                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  2:6122289/93‑1 (MQ=255)

                   cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                                      <  1:4524624/93‑1 (MQ=255)

                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     >  2:280624/1‑93 (MQ=255)

                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     >  1:6654462/1‑93 (MQ=255)

                     tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    >  1:6031879/1‑93 (MQ=255)

                     tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    <  1:5410792/93‑1 (MQ=255)

                      ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                   <  1:4868695/93‑1 (MQ=255)

                      ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                   >  1:4483417/1‑93 (MQ=255)

                         tGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                                                >  2:3236652/1‑93 (MQ=255)

                            cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             >  2:5400569/1‑93 (MQ=255)

                             gggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCggg                                                            >  1:3133139/1‑92 (MQ=255)

                             gggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgt                                                            <  1:3013019/93‑1 (MQ=255)

                              ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                                           >  2:2830553/1‑93 (MQ=255)

                              ggCACCATCCACGACAAAATGCCAACGCCATACGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTAGGCGAGTGAAGACCGCGgtg                                                           <  1:6723633/93‑1 (MQ=255)

                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  1:2103913/93‑1 (MQ=255)

                                cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                                         <  2:2428869/93‑1 (MQ=255)

                                cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                                         <  1:3847640/93‑1 (MQ=255)

                                 aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                        <  2:4212129/93‑1 (MQ=255)

                                  ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       >  2:3940963/1‑93 (MQ=255)

                                   cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      <  1:5312515/93‑1 (MQ=255)

                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     <  2:5509698/93‑1 (MQ=255)

                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  2:81264/1‑93 (MQ=255)

                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  2:3171425/1‑93 (MQ=255)

                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  1:5849552/1‑93 (MQ=255)

                                       cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                                  <  2:255658/93‑1 (MQ=255)

                                          ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               >  1:3061228/1‑93 (MQ=255)

                                          ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               >  1:1017106/1‑93 (MQ=255)

                                          ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               <  2:693492/93‑1 (MQ=255)

                                          ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               <  1:5493087/93‑1 (MQ=255)

                                           gCAAAATGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              <  1:1977680/93‑1 (MQ=255)

                                             aaaaTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            <  1:5758331/93‑1 (MQ=255)

                                             aaaaTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                            >  2:3371324/1‑93 (MQ=255)

                                               aaTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  1:1479010/93‑1 (MQ=255)

                                               aaTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          >  2:1920685/1‑93 (MQ=255)

                                                aTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         >  1:4371503/1‑93 (MQ=255)

                                                aTGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         <  1:172323/93‑1 (MQ=255)

                                                 tGCCAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                                        <  1:4948950/93‑1 (MQ=255)

                                                   ccAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                      >  1:192782/1‑93 (MQ=255)

                                                     aaCGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    >  2:6483052/1‑93 (MQ=255)

                                                      aCGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   >  1:4570053/1‑93 (MQ=255)

                                                       cGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                                  <  2:5602051/93‑1 (MQ=255)

                                                          cATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGa                               >  2:1184838/1‑93 (MQ=255)

                                                            tCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  1:1903911/93‑1 (MQ=255)

                                                             ccGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  2:261166/93‑1 (MQ=255)

                                                             ccGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:2493282/93‑1 (MQ=255)

                                                             ccGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                            <  1:6427403/93‑1 (MQ=255)

                                                              cgtgAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:6719687/90‑1 (MQ=255)

                                                              cGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  1:1641808/93‑1 (MQ=255)

                                                              cGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           >  2:2297910/1‑93 (MQ=255)

                                                              cGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:4443701/93‑1 (MQ=255)

                                                               gggAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg                          <  2:1781227/93‑1 (MQ=255)

                                                                ggAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         <  1:3906473/93‑1 (MQ=255)

                                                                 gAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                        <  1:699878/93‑1 (MQ=255)

                                                                    tACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGggtggt                     >  1:2361172/1‑93 (MQ=255)

                                                                    tACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     <  1:5176586/93‑1 (MQ=255)

                                                                     accaccAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                    >  2:4045822/1‑93 (MQ=255)

                                                                        accaGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 >  1:4059661/1‑93 (MQ=255)

                                                                        accaGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 >  1:338648/1‑93 (MQ=255)

                                                                        accaGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 >  1:296956/1‑93 (MQ=255)

                                                                         ccaGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                >  2:3010537/1‑93 (MQ=255)

                                                                         ccaGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  2:5207954/93‑1 (MQ=255)

                                                                          caGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               >  1:5020479/1‑93 (MQ=255)

                                                                           aGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtgg              >  1:1821956/1‑93 (MQ=255)

                                                                             acgcTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCa            <  2:5076531/90‑1 (MQ=255)

                                                                              cgcTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAg           <  1:6345207/91‑1 (MQ=255)

                                                                              cgcTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGATGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAg           >  2:3557919/3‑93 (MQ=255)

                                                                               gcTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          >  2:4125935/2‑93 (MQ=255)

                                                                               gcTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          <  2:3776674/92‑1 (MQ=255)

                                                                                 tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGc        <  1:4108657/93‑1 (MQ=255)

                                                                                   ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCAt      <  2:5573688/93‑1 (MQ=255)

                                                                                   ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCAt      >  1:4425574/1‑93 (MQ=255)

                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  1:3197449/1‑93 (MQ=255)

                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  2:6087904/1‑93 (MQ=255)

                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     >  1:6201328/1‑93 (MQ=255)

                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAC‑TAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa    <  1:1375646/93‑1 (MQ=255)

                                                                                     tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa    >  2:2126168/1‑93 (MQ=255)

                                                                                      tctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  1:1833217/92‑1 (MQ=255)

                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   >  2:3089014/1‑93 (MQ=255)

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

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