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	99.5%	1264.8 / ‑3.5	375	D123G (GAC→GGC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (185/188); total (186/189)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 412 total aligned reads displayed.

TTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTG  >  NC_000913/3758419‑3758598

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

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

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

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

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

        gACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATa                                                                                 >  1:2797625/1‑93 (MQ=255)

         aCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATAc                                                                                <  2:738798/93‑1 (MQ=255)

         aCGATACAGGCATCTTCCCCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCTCGATGTCCGCGGCCGGCGTCGCCCAGATCATAc                                                                                <  2:6296274/93‑1 (MQ=255)

          cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACt                                                                               <  1:478785/93‑1 (MQ=255)

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

             tACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGa                                                                            >  2:711316/1‑93 (MQ=255)

              aCAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAg                                                                           >  2:20436/1‑93 (MQ=255)

              aCAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAg                                                                           >  2:3417655/1‑93 (MQ=255)

                aGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAggg                                                                         <  2:5762573/93‑1 (MQ=255)

                aGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAggg                                                                         <  1:401144/93‑1 (MQ=255)

                 ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAgggc                                                                        >  1:1945050/1‑92 (MQ=255)

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

                 ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGt                                                                        >  1:7587912/1‑93 (MQ=255)

                    aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGTTCATACGCGAGGGTCAc                                                                     >  1:7817417/1‑93 (MQ=255)

                     tCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGTCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTAGCCCAGATCATACTCGAGGGTCACt                                                                    >  2:5119574/1‑93 (MQ=255)

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

                       ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTgg                                                                  <  1:1663472/93‑1 (MQ=255)

                        tCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCa                                                                                   >  2:1463169/1‑75 (MQ=255)

                         ccGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCg                                                                >  1:4158148/1‑93 (MQ=255)

                          cGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGa                                                               >  1:6376257/1‑93 (MQ=255)

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

                             ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACg                                                            >  1:2290433/1‑93 (MQ=255)

                             ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACg                                                            >  2:3195971/1‑93 (MQ=255)

                             ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACg                                                            <  2:2967421/93‑1 (MQ=255)

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

                                cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCAt                                                         <  2:7913893/93‑1 (MQ=255)

                                cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCAt                                                         >  2:2434179/1‑93 (MQ=255)

                                  tAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAg                                                       >  1:2054137/1‑93 (MQ=255)

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

                                    agccGGCACAGCAGCCGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                                     <  2:1403014/90‑1 (MQ=255)

                                    aTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCTAACGCATAGcc                                                     <  2:8686924/93‑1 (MQ=255)

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

                                     tACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCt                                                    >  2:5289797/1‑93 (MQ=255)

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

                                      aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                                   >  1:449177/1‑93 (MQ=255)

                                       cGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                                  >  2:6798818/1‑93 (MQ=255)

                                        ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcg                                                 >  1:3309879/1‑93 (MQ=255)

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

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

                                           acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                              >  2:5197501/1‑93 (MQ=255)

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

                                           acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                              >  1:7729717/1‑93 (MQ=255)

                                           acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                              >  2:6379347/1‑93 (MQ=255)

                                           acaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                                              >  1:5445597/1‑93 (MQ=255)

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

                                             agcagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACg                                            <  1:7879841/93‑1 (MQ=255)

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

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

                                              gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGAGTCACTAGCGCACGCATAGCCCGCGCAACgg                                           >  1:5140872/1‑93 (MQ=255)

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

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

                                                agcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                                         >  1:8639592/1‑93 (MQ=255)

                                                agcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                                         >  2:7068525/1‑93 (MQ=255)

                                                 gcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt                                        <  2:447117/93‑1 (MQ=255)

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

                                                     cgcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCAc                                    <  2:7623906/93‑1 (MQ=255)

                                                     cgcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCAc                                    >  2:2590682/1‑93 (MQ=255)

                                                      gcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCAcc                                   <  1:4562858/93‑1 (MQ=255)

                                                       ccAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcg                                  >  1:3460596/1‑93 (MQ=255)

                                                        cAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgc                                 <  1:6885920/93‑1 (MQ=255)

                                                        cAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgc                                 >  2:5376677/1‑93 (MQ=255)

                                                         cgcgcCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgcg                                <  1:5310652/92‑1 (MQ=255)

                                                          gcgcCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCgg                               <  2:5850787/93‑1 (MQ=255)

                                                            gcCCGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCt                             >  1:5617713/1‑93 (MQ=255)

                                                            agccGATCCCGCTGTCCGCGTCCGGCGTTGCCCAGAGCAGACTCGATGGTCACTTGCGATAGCATGGCCTGCGCAACGGCTTCCACCGCGGCt                             <  2:1162558/91‑1 (MQ=255)

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

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

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

                                                               cGATCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc                          >  2:821678/1‑93 (MQ=255)

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

                                                                 aTCGCGATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc                        >  1:1008225/1‑93 (MQ=255)

                                                                    gcgATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCtt                     >  1:1732564/1‑92 (MQ=255)

                                                                    gcgATGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTg                     <  2:3648578/93‑1 (MQ=255)

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

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

                                                                       aTGTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaa                  <  1:8272310/93‑1 (MQ=255)

                                                                         gTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAg                <  1:7999312/93‑1 (MQ=255)

                                                                         gTCCGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAg                >  2:4942038/1‑93 (MQ=255)

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

                                                                          tCCGCGTCCGCCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc               <  2:1235580/93‑1 (MQ=255)

                                                                          gccGCGTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGc               <  1:3926527/92‑1 (MQ=255)

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

                                                                             gcgTCCGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTcc            >  2:2285577/1‑92 (MQ=255)

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

                                                                                gccGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGccc         <  2:1613798/92‑1 (MQ=255)

                                                                                 ccGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCa        <  2:6193019/93‑1 (MQ=255)

                                                                                 ccGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCa        >  1:3029875/1‑93 (MQ=255)

                                                                                  cGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCaa       >  2:2759899/1‑93 (MQ=255)

                                                                                  cGGCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCaa       <  2:1911991/93‑1 (MQ=255)

                                                                                  cGGCGTCGCCCAGATCATACTCGAGGGTCACTGACGAACGCATAGCCTGCGCAACGGCTTCCACCGCTGCTTCCGCCTGTAAAGCTCGCCCac       >  1:6271233/1‑92 (MQ=255)

                                                                                    gCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAgg     >  2:391491/1‑93 (MQ=255)

                                                                                    gCGTCGCCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAgg     <  1:2244173/93‑1 (MQ=255)

                                                                                      gtcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGtt   <  2:1754573/93‑1 (MQ=255)

                                                                                       tcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGTGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTg  >  1:1251900/1‑93 (MQ=255)

                                                                                       tcgccCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTg  <  1:8372743/93‑1 (MQ=255)

TTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTG  >  NC_000913/3758419‑3758598

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