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%	776.7 / NA	221	G168D (GGC→GAC)	rhaD	rhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (103/118); total (103/118)

Only 100 of 239 total aligned reads displayed.

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

   gTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                      >  1:1290768/1‑93 (MQ=255)

   gTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                      >  2:7243807/1‑93 (MQ=255)

   gTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                      >  1:2843144/1‑93 (MQ=255)

    tCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                                     >  1:8715692/1‑93 (MQ=255)

       ccTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                                  >  2:3452046/1‑93 (MQ=255)

       ccTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                                  >  2:4591834/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                <  1:5343916/93‑1 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  2:749959/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  1:5449709/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                >  2:5284189/1‑93 (MQ=255)

         tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGCGCTGTCTTCCCACAGt                                                                                >  2:8111549/1‑93 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATCCCACCAGACACACGGTGCTGTCTTCCCACAGtt                                                                               <  2:6699211/93‑1 (MQ=255)

          ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                               <  1:8344224/93‑1 (MQ=255)

           gCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTg                                                                              >  1:7259859/1‑93 (MQ=255)

            ccGATTTCGTCCGTGCCGGGCACCATCCAGGGCAAAATGCCAACGCCATCCGGGAATACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             <  1:7053862/93‑1 (MQ=255)

            ccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             <  1:6575547/93‑1 (MQ=255)

            ccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                                             >  1:1466366/1‑93 (MQ=255)

              gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                           <  1:3942583/93‑1 (MQ=255)

                  tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                                       <  2:3178505/93‑1 (MQ=255)

                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACTCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  1:5286664/93‑1 (MQ=255)

                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  2:199934/93‑1 (MQ=255)

                    gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  2:3719871/93‑1 (MQ=255)

                    gTCCGTGCCGCGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                                     <  2:3335784/93‑1 (MQ=255)

                     tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                    >  2:4691876/1‑93 (MQ=255)

                      ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                   >  2:5788606/1‑93 (MQ=255)

                          gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                               <  2:1378380/93‑1 (MQ=255)

                          gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                               <  1:5199199/93‑1 (MQ=255)

                            cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCgg                                                             >  2:6429718/1‑93 (MQ=255)

                              ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                                           >  1:7984153/1‑93 (MQ=255)

                               gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                                          <  2:7774237/93‑1 (MQ=255)

                                 aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                        >  1:4324964/1‑93 (MQ=255)

                                 aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                        <  1:2766745/93‑1 (MQ=255)

                                  ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       >  2:5257190/1‑93 (MQ=255)

                                  ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       >  2:2295610/1‑93 (MQ=255)

                                  ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       <  2:6996972/93‑1 (MQ=255)

                                  ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                       <  2:8083045/93‑1 (MQ=255)

                                   cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      >  1:614876/1‑93 (MQ=255)

                                   cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                      >  2:6179011/1‑93 (MQ=255)

                                    ctCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     <  2:5936880/92‑1 (MQ=255)

                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     >  1:5931844/1‑93 (MQ=255)

                                    aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                     >  1:806892/1‑93 (MQ=255)

                                     tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                    >  1:1104699/1‑93 (MQ=255)

                                       cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                                  <  2:7350954/93‑1 (MQ=255)

                                          ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                               >  2:1065286/1‑93 (MQ=255)

                                           gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                              >  1:6864103/1‑93 (MQ=255)

                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  1:800518/93‑1 (MQ=255)

                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          <  1:7547309/93‑1 (MQ=255)

                                               aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                          >  2:692183/1‑93 (MQ=255)

                                                aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         >  2:6274987/1‑93 (MQ=255)

                                                aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                         <  2:2568463/93‑1 (MQ=255)

                                                  cccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                       <  2:6556451/92‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACGCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:2150679/93‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCGCAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:2931825/93‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  2:2913353/93‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  2:6580851/93‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:6014050/93‑1 (MQ=255)

                                                    cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                                     <  1:701323/93‑1 (MQ=255)

                                                     aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                    <  2:1832359/93‑1 (MQ=255)

                                                      aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   <  2:503635/93‑1 (MQ=255)

                                                      aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   <  2:4853354/93‑1 (MQ=255)

                                                      aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                                   >  2:4798602/1‑93 (MQ=255)

                                                         ccgtCCGAGATTAACACCAGACACTAAGAGCTGTCTTCCCAAAGTAGGCAAGTAAAGACCGCGGTGACGTTTTAAAATACATAGGTGAGGGCg                                <  2:8059324/90‑1 (MQ=255)

                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  1:2621722/93‑1 (MQ=255)

                                                            tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                             <  1:3135893/93‑1 (MQ=255)

                                                            tCCGGGAATACCACCAGACACGCGGGGCTGTCTTCCCTCAGTTGGCGAGTGAAGACCGCGGAGTCGGGGTCAAGTACATAGGTGAGGGCGATc                             >  2:7842500/1‑93 (MQ=255)

                                                              cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:7627671/93‑1 (MQ=255)

                                                              cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:7978425/93‑1 (MQ=255)

                                                              cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                           <  2:5003180/93‑1 (MQ=255)

                                                               gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg                          >  1:7915556/1‑93 (MQ=255)

                                                               gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg                          <  1:1361296/93‑1 (MQ=255)

                                                                ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt                         >  2:1935606/1‑93 (MQ=255)

                                                                 gAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                        <  2:7030200/93‑1 (MQ=255)

                                                                  aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                       >  2:5100861/1‑93 (MQ=255)

                                                                  aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                       >  1:4897432/1‑93 (MQ=255)

                                                                  aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                       <  2:3491135/93‑1 (MQ=255)

                                                                   aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg                      >  2:4011046/1‑93 (MQ=255)

                                                                    tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTGGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggg                     >  1:2727616/1‑92 (MQ=255)

                                                                    tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                     <  1:5748649/93‑1 (MQ=255)

                                                                      ccaccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtgg                   <  2:2957929/93‑1 (MQ=255)

                                                                        cccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 <  2:4237431/92‑1 (MQ=255)

                                                                        accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 >  1:2528108/1‑93 (MQ=255)

                                                                        accaGACACTCAGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                 >  2:447430/1‑93 (MQ=255)

                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgt                <  1:4858349/93‑1 (MQ=255)

                                                                         ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgg                >  1:293184/1‑92 (MQ=255)

                                                                          caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg               >  1:3238118/1‑93 (MQ=255)

                                                                           aGACTCTCGGTGATGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtgg              <  1:3750777/93‑1 (MQ=255)

                                                                            gACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc             <  2:8468804/93‑1 (MQ=255)

                                                                               acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          >  2:8493100/1‑93 (MQ=255)

                                                                               acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt          <  1:4728676/93‑1 (MQ=255)

                                                                                cTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTg         >  1:3331142/1‑93 (MQ=255)

                                                                                  cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa       >  1:3664285/1‑93 (MQ=255)

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

                                                                                    gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATg     <  2:3300860/93‑1 (MQ=255)

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

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

                                                                                     tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa    <  1:8532944/93‑1 (MQ=255)

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

                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt   <  1:7326241/93‑1 (MQ=255)

                                                                                      gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGAATGAt   >  1:6644419/1‑93 (MQ=255)

                                                                                       ctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc  <  1:1041428/93‑1 (MQ=255)

GCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATC  >  NC_000913/4091261‑4091440

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