BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,758,509	T→C	D123G (GAC→GGC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,758,509	0	T	C	83.3%	624.9 / 121.4	287	D123G (GAC→GGC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base T (29/19); new base C (119/120); total (148/139)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 316 total aligned reads displayed.

CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGT  >  NC_000913/3758417‑3758600

  ttGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGccc                                                                                           >  1:1006978/1‑93 (MQ=255)

  ttGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGccc                                                                                           >  2:5546463/1‑93 (MQ=255)

  ttGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGccc                                                                                           <  1:2520284/93‑1 (MQ=255)

   tGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCa                                                                                          <  2:4392400/93‑1 (MQ=255)

   tGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCa                                                                                          <  1:6247295/93‑1 (MQ=255)

    gTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg                                                                                         >  1:4760001/1‑93 (MQ=255)

      tATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGAt                                                                                       >  1:7218862/1‑93 (MQ=255)

       aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATc                                                                                      <  1:112829/93‑1 (MQ=255)

       aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATc                                                                                      >  1:5109912/1‑93 (MQ=255)

        ttGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCa                                                                                     >  2:4949266/1‑93 (MQ=255)

          gACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATa                                                                                   <  2:3612881/93‑1 (MQ=255)

          gACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATa                                                                                   <  1:3314604/93‑1 (MQ=255)

             gATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTc                                                                                <  1:290239/93‑1 (MQ=255)

             gATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTc                                                                                >  1:2264340/1‑93 (MQ=255)

              aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCg                                                                               >  2:583057/1‑93 (MQ=255)

               tACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGa                                                                              <  2:1261535/93‑1 (MQ=255)

               tACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGa                                                                              <  2:3198249/93‑1 (MQ=255)

                 cAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAgg                                                                            >  1:1624921/1‑93 (MQ=255)

                   ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGt                                                                          <  1:7129509/93‑1 (MQ=255)

                   ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGt                                                                          <  1:6658372/93‑1 (MQ=255)

                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                                       >  1:130604/1‑93 (MQ=255)

                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCAc                                                                       >  1:7494075/1‑93 (MQ=255)

                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCAc                                                                       <  1:952979/93‑1 (MQ=255)

                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCAc                                                                       <  2:1520426/93‑1 (MQ=255)

                       tCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACt                                                                      >  2:4163234/1‑93 (MQ=255)

                       tCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACt                                                                      >  1:4395027/1‑93 (MQ=255)

                       tCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACt                                                                      >  2:6247330/1‑93 (MQ=255)

                        tttCCCCCCCCGTAATCCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTg                                                                     <  1:6151114/92‑1 (MQ=255)

                        cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTg                                                                     <  2:4038150/93‑1 (MQ=255)

                        cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTg                                                                     <  1:5767130/93‑1 (MQ=255)

                           ccGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCg                                                                  <  1:3664210/93‑1 (MQ=255)

                            cGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                                 >  2:1016131/1‑93 (MQ=255)

                            cGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGa                                                                 <  2:6291938/93‑1 (MQ=255)

                             gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGaa                                                                <  1:4876412/93‑1 (MQ=255)

                              cccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAAc                                                               <  1:3136454/93‑1 (MQ=255)

                              cccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAAc                                                               >  2:1350649/1‑93 (MQ=255)

                                cccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGc                                                             >  1:4553944/1‑93 (MQ=255)

                                 ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                                            >  2:3743082/1‑93 (MQ=255)

                                   gTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                                                          >  1:4474242/1‑93 (MQ=255)

                                   gTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                                                          >  2:2641662/1‑93 (MQ=255)

                                    tAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAg                                                         <  1:2410388/93‑1 (MQ=255)

                                    tAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGCCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAg                                                         <  2:5892993/93‑1 (MQ=255)

                                     aaTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                                                        >  2:7296612/1‑93 (MQ=255)

                                      aTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                                       <  1:215582/93‑1 (MQ=255)

                                      aTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                                       >  2:7404954/1‑93 (MQ=255)

                                      aTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                                       <  1:6102682/93‑1 (MQ=255)

                                        aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                                     <  1:2762789/93‑1 (MQ=255)

                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                                  >  2:2994632/1‑93 (MQ=255)

                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                                  >  2:3455929/1‑93 (MQ=255)

                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                                  <  1:4599944/93‑1 (MQ=255)

                                             acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                                >  1:427985/1‑93 (MQ=255)

                                             acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                                <  1:7367600/93‑1 (MQ=255)

                                              cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAAc                                               <  2:1600575/93‑1 (MQ=255)

                                                gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACgg                                             >  1:1866894/1‑93 (MQ=255)

                                                 cagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGc                                            >  1:125399/1‑93 (MQ=255)

                                                 cagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGc                                            <  1:1787423/93‑1 (MQ=255)

                                                  agcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                                           <  1:417569/93‑1 (MQ=255)

                                                   gcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt                                          >  2:2329104/1‑93 (MQ=255)

                                                   gcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt                                          >  1:6457663/1‑93 (MQ=255)

                                                    cTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTc                                         >  2:3834956/1‑93 (MQ=255)

                                                      gcgcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCa                                       <  1:653443/93‑1 (MQ=255)

                                                           aGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgcg                                  >  2:1741353/1‑93 (MQ=255)

                                                           aGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgcg                                  >  1:196957/1‑93 (MQ=255)

                                                           aGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgcg                                  >  2:2549331/1‑93 (MQ=255)

                                                               cccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt                              >  2:6541271/1‑93 (MQ=255)

                                                               cccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt                              >  1:3471208/1‑93 (MQ=255)

                                                               cccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt                              <  1:587568/93‑1 (MQ=255)

                                                                ccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTc                             >  1:683401/1‑93 (MQ=255)

                                                                ccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTc                             <  1:4002657/93‑1 (MQ=255)

                                                                ccGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTc                             <  1:6329162/93‑1 (MQ=255)

                                                                 cGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc                            >  1:461074/1‑93 (MQ=255)

                                                                 cGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc                            <  1:3784953/93‑1 (MQ=255)

                                                                   aTCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc                          <  2:791392/93‑1 (MQ=255)

                                                                    tCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGcc                         >  2:4296599/1‑93 (MQ=255)

                                                                     cgcgATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCt                        <  2:548166/93‑1 (MQ=255)

                                                                       cgATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGt                      <  1:1503708/93‑1 (MQ=255)

                                                                        gATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTa                     <  2:6714064/93‑1 (MQ=255)

                                                                        gATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTa                     >  2:1813047/1‑93 (MQ=255)

                                                                         aTGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCGCCGCGGCTTCCGCCTGTaa                    <  1:803794/93‑1 (MQ=255)

                                                                          tGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa                   <  1:2034173/93‑1 (MQ=255)

                                                                          tGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa                   >  2:276305/1‑93 (MQ=255)

                                                                            tacgcgTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc                 <  2:3936242/91‑1 (MQ=255)

                                                                            tCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc                 <  2:1086192/93‑1 (MQ=255)

                                                                            tCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc                 >  2:5461124/1‑93 (MQ=255)

                                                                            tCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc                 >  1:5083592/1‑93 (MQ=255)

                                                                             ccGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCt                <  2:2582085/93‑1 (MQ=255)

                                                                              cgcgTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTc               <  2:5274118/93‑1 (MQ=255)

                                                                               gcgTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCg              <  2:4570459/93‑1 (MQ=255)

                                                                               gcgTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCg              <  2:3041622/93‑1 (MQ=255)

                                                                               gcgTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCg              >  2:2257277/1‑93 (MQ=255)

                                                                                 gTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGcc            >  2:6068570/1‑93 (MQ=255)

                                                                                 gTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGcc            <  1:392384/93‑1 (MQ=255)

                                                                                  tCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGccc           >  1:4858009/1‑93 (MQ=255)

                                                                                  tCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGccc           >  1:5637221/1‑93 (MQ=255)

                                                                                      gCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAgg       >  1:1404337/1‑93 (MQ=255)

                                                                                       cgtcgcCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGt      >  1:3523605/1‑93 (MQ=255)

                                                                                        gtcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGtt     >  2:3305159/1‑93 (MQ=255)

                                                                                         tcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTg    >  1:5162836/1‑93 (MQ=255)

                                                                                         tcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTg    >  2:7059329/1‑93 (MQ=255)

                                                                                           gtcCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGt  >  2:687300/1‑93 (MQ=255)

CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGT  >  NC_000913/3758417‑3758600

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