BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	641,948	A→G	Q213R (CAG→CGG)	ybdR →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	641,948	0	A	G	100.0%	424.7 / NA	134	Q213R (CAG→CGG)	ybdR	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (71/63); total (71/63)

Only 100 of 140 total aligned reads displayed.

AGATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCG  >  NC_000913/641858‑642037

agagCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGa                                                                                         <  1:826118/93‑1 (MQ=255)

aGATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGa                                                                                         >  1:2828380/1‑93 (MQ=255)

aGATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGa                                                                                         >  2:2621306/1‑93 (MQ=255)

aGATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGa                                                                                         <  2:2382164/93‑1 (MQ=255)

 gATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGAt                                                                                        <  1:3955768/93‑1 (MQ=255)

 gATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGAt                                                                                        <  1:2659441/93‑1 (MQ=255)

 gATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGAt                                                                                        >  1:869828/1‑93 (MQ=255)

   tCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGAttt                                                                                      >  1:1952015/1‑93 (MQ=255)

      aacaaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttg                                                                                   >  2:2357613/1‑93 (MQ=255)

      aacaaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttg                                                                                   <  2:1894640/93‑1 (MQ=255)

       acaaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgt                                                                                  >  2:2051296/1‑93 (MQ=255)

       acaaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgt                                                                                  >  1:2463637/1‑93 (MQ=255)

        caaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgtt                                                                                 >  1:1549293/1‑93 (MQ=255)

        caaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgtt                                                                                 >  2:2565755/1‑93 (MQ=255)

        caaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgtt                                                                                 >  2:2351768/1‑93 (MQ=255)

        caaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgtt                                                                                 <  2:2293145/93‑1 (MQ=255)

        caaGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgtt                                                                                 <  1:1417031/93‑1 (MQ=255)

          aGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGAGAATCGCCTGTGCACGGTTGCTCGGTGCGGGAGGGATTTttgttgt                                                                               >  1:3427537/1‑93 (MQ=255)

          aGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgt                                                                               <  1:2232986/93‑1 (MQ=255)

          aGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgt                                                                               <  2:1736147/93‑1 (MQ=255)

          aGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgt                                                                               <  1:859925/93‑1 (MQ=255)

          aGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgt                                                                               >  1:3243879/1‑93 (MQ=255)

            gTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgttg                                                                             <  1:1808958/93‑1 (MQ=255)

            gTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTttgttgttg                                                                             >  2:4024514/1‑93 (MQ=255)

             ttCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGa                                                                            <  2:3162526/93‑1 (MQ=255)

              tCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGat                                                                           >  2:1768428/1‑93 (MQ=255)

              tCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGat                                                                           >  1:1388009/1‑93 (MQ=255)

                aaGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatca                                                                         <  2:1728193/93‑1 (MQ=255)

                aaGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatca                                                                         >  2:2198170/1‑93 (MQ=255)

                 aGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcat                                                                        >  1:2560538/1‑93 (MQ=255)

                  gCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatc                                                                       <  1:3118771/93‑1 (MQ=255)

                  gCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatc                                                                       >  1:3488171/1‑93 (MQ=255)

                    gTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatcat                                                                     >  2:229358/1‑93 (MQ=255)

                    gTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatcat                                                                     <  1:806466/93‑1 (MQ=255)

                    gTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatcat                                                                     >  1:1998002/1‑93 (MQ=255)

                     ttGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGatcatcatc                                                                    >  2:3391132/1‑93 (MQ=255)

                      tGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATcc                                                                   <  1:278594/93‑1 (MQ=255)

                      tGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATcc                                                                   <  1:4227179/93‑1 (MQ=255)

                       gCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTTATTATCAGccc                                                                  <  2:4565211/93‑1 (MQ=255)

                       gCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATccc                                                                  <  2:2750591/93‑1 (MQ=255)

                         aGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTa                                                                <  2:4636930/93‑1 (MQ=255)

                         aGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTa                                                                <  2:3387833/93‑1 (MQ=255)

                          gTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAc                                                               <  1:2984741/93‑1 (MQ=255)

                          gTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAc                                                               >  1:1522677/1‑93 (MQ=255)

                          gTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAc                                                               >  1:3403976/1‑93 (MQ=255)

                           tCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAcc                                                              <  2:3436680/93‑1 (MQ=255)

                           tCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAcc                                                              >  2:15925/1‑93 (MQ=255)

                           tCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTAcc                                                              <  2:1139796/93‑1 (MQ=255)

                             tATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGc                                                            >  1:723670/1‑93 (MQ=255)

                             tATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGc                                                            >  1:3275188/1‑93 (MQ=255)

                               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCtt                                                          >  1:1476380/1‑93 (MQ=255)

                               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCtt                                                          >  1:977919/1‑93 (MQ=255)

                               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGACTGTGCACGGTTGCTCGGTGCGGCACGGATTTTTGTTGTTGATCATCATCCCTACCGCtt                                                          >  1:1619575/1‑93 (MQ=255)

                                 gTGCTGGTCCTGTGGGATTGTTGACAATCGGCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGc                                                        <  1:4328322/93‑1 (MQ=255)

                                   gCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCAt                                                      >  2:3894263/1‑93 (MQ=255)

                                     tGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCAttt                                                    <  1:650390/93‑1 (MQ=255)

                                      ggTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTc                                                   <  1:4400570/93‑1 (MQ=255)

                                        tCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTcgc                                                 >  2:318693/1‑93 (MQ=255)

                                           tgtgGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTcgccgc                                              >  1:2453675/1‑93 (MQ=255)

                                           tgtgGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTcgccgc                                              >  1:1789041/1‑93 (MQ=255)

                                              gggATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGa                                           >  2:614149/1‑93 (MQ=255)

                                                gATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCCTCCCTACCGCTTGCATTTCGCCGCCGAcc                                         <  1:281837/93‑1 (MQ=255)

                                                gATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGAcc                                         >  2:1621821/1‑93 (MQ=255)

                                                 aTTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCg                                        >  1:1211041/1‑93 (MQ=255)

                                                  ttgttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGc                                       <  1:1656791/93‑1 (MQ=255)

                                                  ttgttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGc                                       >  2:3587362/1‑93 (MQ=255)

                                                    gttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTa                                     <  1:4482076/93‑1 (MQ=255)

                                                     ttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTAc                                    >  1:2148400/1‑93 (MQ=255)

                                                     ttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTAc                                    <  1:41852/93‑1 (MQ=255)

                                                     ttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTAc                                    <  2:2337370/93‑1 (MQ=255)

                                                       gACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACgg                                  <  1:4615309/93‑1 (MQ=255)

                                                           aTCGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGcgcg                              >  1:4032607/1‑93 (MQ=255)

                                                             cGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGAt                            >  1:3210936/1‑93 (MQ=255)

                                                             cGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGAt                            <  2:4456717/93‑1 (MQ=255)

                                                             cGCCTGTGCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGAt                            <  1:4009550/93‑1 (MQ=255)

                                                                  gtgCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGa                       <  2:2075031/93‑1 (MQ=255)

                                                                   tgCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGAt                      >  1:2650184/1‑93 (MQ=255)

                                                                   tgCACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGAt                      <  2:1310114/93‑1 (MQ=255)

                                                                     cACGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTa                    <  1:328056/93‑1 (MQ=255)

                                                                      aCGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTaa                   >  2:1473412/1‑93 (MQ=255)

                                                                      aCGGTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTaa                   <  1:3025416/93‑1 (MQ=255)

                                                                        ggTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAAtt                 >  2:2448186/1‑93 (MQ=255)

                                                                        ggTTGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAAtt                 >  1:915072/1‑93 (MQ=255)

                                                                          ttGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAAtttt               <  1:3618277/93‑1 (MQ=255)

                                                                          ttGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAAtttt               <  2:4492681/93‑1 (MQ=255)

                                                                           tGCTCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTtg              >  1:1280794/1‑93 (MQ=255)

                                                                              tCGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTtgatg           >  2:452809/1‑93 (MQ=255)

                                                                               cGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTtgatga          <  1:273833/93‑1 (MQ=255)

                                                                               cGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTtgatga          <  1:946298/93‑1 (MQ=255)

                                                                               cGGTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTtgatga          <  2:1206702/93‑1 (MQ=255)

                                                                                ggTGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGaa         >  2:1582282/1‑93 (MQ=255)

                                                                                  tGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGa       <  1:2722658/93‑1 (MQ=255)

                                                                                  tGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGa       >  1:3784367/1‑93 (MQ=255)

                                                                                  tGCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGa       <  2:459760/93‑1 (MQ=255)

                                                                                   gCGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGAc      <  1:1006052/93‑1 (MQ=255)

                                                                                    cGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACa     <  2:1119184/93‑1 (MQ=255)

                                                                                    cGGAACGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACa     >  1:2397112/1‑93 (MQ=255)

                                                                                       aaCGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCg  >  1:1722072/1‑93 (MQ=255)

                                                                                       aaCGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCg  >  2:3893041/1‑93 (MQ=255)

                                                                                       aaCGGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCg  <  2:3617874/93‑1 (MQ=255)

AGATCCAACAAGGTTCAAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCCGATTAATTTTGATGAAGACAGCG  >  NC_000913/641858‑642037

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