BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,925	T→C	W305R (TGG→CGG)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,925	0	T	C	100.0%	1359.2 / NA	324	W305R (TGG→CGG)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (162/162); total (162/162)

Only 100 of 324 total aligned reads displayed.

AGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGT  >  NC_000913/3103827‑3104022

aGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCgg                                                                                                 <  2:5065114/101‑1 (MQ=255)

aGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCgg                                                                                                 <  2:5064877/101‑1 (MQ=255)

 gTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGc                                                                                                <  2:5065029/101‑1 (MQ=255)

 gTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGc                                                                                                <  1:5065017/101‑1 (MQ=255)

  tggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCt                                                                                               <  2:5065255/101‑1 (MQ=255)

   ggctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCtt                                                                                              <  2:5064900/101‑1 (MQ=255)

    gctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTc                                                                                             <  1:5065035/101‑1 (MQ=255)

     ctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTcc                                                                                            <  1:5064977/101‑1 (MQ=255)

         cgcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCgtg                                                                                        <  1:5064859/101‑1 (MQ=255)

         cgcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCgtg                                                                                        >  1:5065516/1‑101 (MQ=255)

         cgcAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCgtg                                                                                        >  2:5065514/1‑101 (MQ=255)

           cAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTgtc                                                                                      >  2:5065519/1‑101 (MQ=255)

             aCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTgtcgt                                                                                    <  2:5065010/101‑1 (MQ=255)

               ggCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCa                                                                                  >  2:5065524/1‑101 (MQ=255)

                  aGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTc                                                                               <  2:5064847/101‑1 (MQ=255)

                    aTTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCAc                                                                             <  2:5065138/101‑1 (MQ=255)

                    aTTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCAc                                                                             >  1:5065527/1‑101 (MQ=255)

                     ttGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCAcc                                                                            >  2:5065528/1‑101 (MQ=255)

                      tgctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCg                                                                           >  1:5065531/1‑101 (MQ=255)

                         tgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCggat                                                                        <  1:5065208/101‑3 (MQ=255)

                           ccGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGc                                                                      >  1:5065537/1‑101 (MQ=255)

                            cGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCa                                                                     <  1:5064686/101‑1 (MQ=255)

                             gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                    >  2:5065542/1‑101 (MQ=255)

                             gATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                    >  1:5065541/1‑101 (MQ=255)

                              aTAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATg                                                                   <  2:5065383/101‑1 (MQ=255)

                               tAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                  >  1:5065546/1‑101 (MQ=255)

                                  ccTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                               >  2:5065551/1‑101 (MQ=255)

                                   cTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                              >  1:5065553/1‑101 (MQ=255)

                                    tGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGaa                                                             >  2:5065556/1‑101 (MQ=255)

                                      aCGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                           <  2:5064712/101‑1 (MQ=255)

                                       cGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGc                                                          >  2:5065559/1‑101 (MQ=255)

                                       cGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGc                                                          >  2:5065563/1‑101 (MQ=255)

                                        gCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCa                                                         <  2:5065209/101‑1 (MQ=255)

                                           aCTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATg                                                      >  1:5065569/1‑101 (MQ=255)

                                            cTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgc                                                     <  1:5064798/101‑1 (MQ=255)

                                            cTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgc                                                     <  2:5065289/101‑1 (MQ=255)

                                            cTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgc                                                     <  1:5065064/101‑1 (MQ=255)

                                             tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgcg                                                    >  1:5065574/1‑101 (MQ=255)

                                             tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgcg                                                    <  1:5064894/101‑1 (MQ=255)

                                             tGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgcg                                                    >  2:5065573/1‑101 (MQ=255)

                                               aCCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTgtcgtc                                                                                   <  2:5065179/68‑1 (MQ=255)

                                               aCCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGct                                                  <  2:5065187/101‑1 (MQ=255)

                                                ccGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctc                                                 >  1:5065579/1‑101 (MQ=255)

                                                 cgcgTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGctct                                                <  1:5064817/101‑1 (MQ=255)

                                                   cgTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTgg                                              <  2:5064939/101‑1 (MQ=255)

                                                    gTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGt                                             >  2:5065580/1‑101 (MQ=255)

                                                    gTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGt                                             <  1:5065153/101‑1 (MQ=255)

                                                     tttCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGta                                            <  1:5064778/101‑1 (MQ=255)

                                                       tCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGct                                                  >  2:5065585/1‑93 (MQ=255)

                                                        cGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATaa                                         <  1:5065310/101‑1 (MQ=255)

                                                        cGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATaa                                         <  2:5065001/101‑1 (MQ=255)

                                                         ggCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAAc                                        >  2:5065590/1‑101 (MQ=255)

                                                         ggCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAAc                                        <  1:5065100/101‑1 (MQ=255)

                                                         ggCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAAc                                        >  1:5065587/1‑101 (MQ=255)

                                                          gCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACt                                       <  2:5065040/101‑1 (MQ=255)

                                                           cATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                      <  1:5065196/101‑1 (MQ=255)

                                                           cATACCTTCAGCCATTTCCACTTAGATATTGTGCCCATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                      <  1:5064962/101‑1 (MQ=255)

                                                            aTACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTa                                     <  1:5064971/101‑1 (MQ=255)

                                                            aTACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTa                                     >  2:5065594/1‑101 (MQ=255)

                                                             tACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                    <  1:5065023/101‑1 (MQ=255)

                                                             tACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAg                                    >  1:5065597/1‑101 (MQ=255)

                                                              aCCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgc                                   >  2:5065602/1‑101 (MQ=255)

                                                              aCCTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgc                                   <  1:5065476/101‑1 (MQ=255)

                                                                cTTCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcgc                                 <  1:5064925/101‑1 (MQ=255)

                                                                 ttCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCa                                >  1:5065605/1‑101 (MQ=255)

                                                                  tCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                               >  2:5065610/1‑101 (MQ=255)

                                                                  tCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                               >  2:5065606/1‑101 (MQ=255)

                                                                  tCAGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCaa                               <  1:5065252/101‑1 (MQ=255)

                                                                    aGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                             <  2:5065117/101‑1 (MQ=255)

                                                                    aGCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                             <  1:5065230/101‑1 (MQ=255)

                                                                     gCCATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccg                            >  1:5065614/1‑101 (MQ=255)

                                                                      ccATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                           <  2:5064928/69‑1 (MQ=255)

                                                                      ccATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgc                           >  2:5065618/1‑101 (MQ=255)

                                                                       cATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgcc                          >  2:5065622/1‑101 (MQ=255)

                                                                       cATTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgcc                          >  1:5065623/1‑101 (MQ=255)

                                                                        aTTTCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgcct                         >  1:5065626/1‑100 (MQ=255)

                                                                          gtCCACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTATCGCAACCGCCGTc                       >  1:5065629/2‑101 (MQ=255)

                                                                             cACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAg                     <  1:5065565/100‑1 (MQ=255)

                                                                             cACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGt                    >  2:5065634/1‑101 (MQ=255)

                                                                             cACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGt                    >  2:5065637/1‑101 (MQ=255)

                                                                             cACTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGt                    <  2:5064959/101‑1 (MQ=255)

                                                                              aCTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                   <  1:5064870/101‑1 (MQ=255)

                                                                               cTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTg                  >  2:5065641/1‑101 (MQ=255)

                                                                               cTTAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTg                  <  2:5065019/101‑1 (MQ=255)

                                                                                ttAGATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                   >  1:5065642/1‑99 (MQ=255)

                                                                                   gATATTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCt              >  2:5065648/1‑101 (MQ=255)

                                                                                    atatTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCgg                                                                          >  1:5065652/1‑40 (MQ=255)

                                                                                    atatTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACGGCCGTCAGTTGGCCTa             <  1:5064986/101‑1 (MQ=255)

                                                                                    atatTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTa             <  2:5065493/101‑1 (MQ=255)

                                                                                    atatTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTa             >  1:5065649/1‑101 (MQ=255)

                                                                                     tatTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAg            <  1:5064996/101‑1 (MQ=255)

                                                                                      atTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGc           >  2:5065657/1‑101 (MQ=255)

                                                                                      atTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGc           >  1:5065658/1‑101 (MQ=255)

                                                                                      atTGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGc           <  1:5065132/101‑1 (MQ=255)

                                                                                       ttGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg          >  1:5065662/1‑101 (MQ=255)

                                                                                       ttGTGCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCg          <  1:5064949/101‑1 (MQ=255)

                                                                                          tgCCTATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCt       >  2:5065665/1‑101 (MQ=255)

                                                                                              tATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATg                                                                   <  2:5064968/37‑1 (MQ=255)

                                                                                              tATGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg   <  2:5065098/101‑1 (MQ=255)

                                                                                               aTGCGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt  >  1:5065671/1‑101 (MQ=255)

AGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGT  >  NC_000913/3103827‑3104022

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