BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,506,868	A→G	P193P (CCT→CCC)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,506,868	0	A	G	100.0%	651.8 / NA	193	P193P (CCT→CCC)	glk	glucokinase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (92/101); total (92/101)

Only 100 of 204 total aligned reads displayed.

CAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGC  >  NC_000913/2506778‑2506957

cagcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCggg                                                                                         >  2:3014845/1‑93 (MQ=255)

cagcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATCGCGCGATACAAATTCACCAGCCCggg                                                                                         >  1:5200557/1‑93 (MQ=255)

 agcGCGCGTTCTGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCgggg                                                                                        <  1:1939066/93‑1 (MQ=255)

 agcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCgggg                                                                                        >  1:936809/1‑93 (MQ=255)

  gcgcgcGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGc                                                                                       >  2:4304373/1‑93 (MQ=255)

   cgcgcgTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGcc                                                                                      >  2:3101230/1‑93 (MQ=255)

    gcgcgTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCa                                                                                     <  1:331799/93‑1 (MQ=255)

     cgcgTTCGATAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGGCCGGGGCCAg                                                                                    <  1:514894/93‑1 (MQ=255)

      gcgTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGa                                                                                   >  1:3522674/1‑93 (MQ=255)

       cgTTCGGTTATATCTTTTGGCTTGAGATTTTCTGGCATGCGGTTGTCAGCTTTCACAATTGCGCGGTACAAATTCACTCGCCCCGGGCCAGaa                                                                                  >  1:391686/1‑93 (MQ=255)

        gTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGaaa                                                                                 >  2:496050/1‑93 (MQ=255)

         ttCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAg                                                                                >  1:3801963/1‑93 (MQ=255)

           cGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCa                                                                              >  1:5623765/1‑93 (MQ=255)

           cGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCa                                                                              <  2:16898/93‑1 (MQ=255)

            ggTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAc                                                                             <  2:4460502/93‑1 (MQ=255)

            ggTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAc                                                                             <  1:3108701/93‑1 (MQ=255)

            ggTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAc                                                                             >  2:64803/1‑93 (MQ=255)

               aaTATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgcg                                                                          <  2:6059420/93‑1 (MQ=255)

                atatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACgcgc                                                                         <  2:1764113/93‑1 (MQ=255)

                 tatCTTTTGTCTTGAGATTATCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCt                                                                        >  1:6273863/1‑93 (MQ=255)

                 tatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCt                                                                        >  2:4012150/1‑93 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  2:1356918/93‑1 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  1:1841124/93‑1 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  1:965422/93‑1 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       >  2:3835654/1‑93 (MQ=255)

                    cTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccg                                                                     >  2:6350652/1‑93 (MQ=255)

                     ttttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgc                                                                    >  2:1465976/1‑93 (MQ=255)

                     ttttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgc                                                                    <  2:4393635/93‑1 (MQ=255)

                     ttttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAAGCGCTccgc                                                                    <  1:711918/93‑1 (MQ=255)

                         ggCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGaa                                                                >  2:2229568/1‑93 (MQ=255)

                          gCTTGAGATTTTCTCGCCGGCGGGTGTCAGCTTTCACAATTGCGCGATACAAATTCAACAGCCCGGGGACCGAAAGCACGCGCTCCGCCGaaa                                                               >  2:4044178/1‑93 (MQ=255)

                           cTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAAc                                                              >  1:563154/1‑93 (MQ=255)

                           cTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAAc                                                              >  2:4646054/1‑93 (MQ=255)

                            ttGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACa                                                             >  1:2501208/1‑93 (MQ=255)

                            ttGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACa                                                             <  2:5541800/93‑1 (MQ=255)

                              gagaTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGGCAGAAAGCACGCGCTCCGACGAAACATg                                                           >  2:1095378/1‑93 (MQ=255)

                               agaTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGa                                                          >  2:6146814/1‑93 (MQ=255)

                                gaTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGAc                                                         >  2:2248063/1‑93 (MQ=255)

                                gaTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGAc                                                         >  2:2746314/1‑93 (MQ=255)

                                  ttttCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCa                                                       >  1:2615354/1‑93 (MQ=255)

                                     tCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAtt                                                    <  2:5134189/93‑1 (MQ=255)

                                     tCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAtt                                                    <  2:6283282/93‑1 (MQ=255)

                                       tGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTc                                                  <  1:4656136/93‑1 (MQ=255)

                                        ggCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTcc                                                 >  1:4501043/1‑93 (MQ=255)

                                           aGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgca                                              <  2:2839026/93‑1 (MQ=255)

                                            ggCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgcac                                             >  2:2082792/1‑93 (MQ=255)

                                            ggCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgcac                                             >  2:6425736/1‑93 (MQ=255)

                                               ggTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgcacgca                                          >  1:3430816/1‑93 (MQ=255)

                                                gTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCaa                                         <  2:3772653/93‑1 (MQ=255)

                                                gTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCaa                                         >  2:2931872/1‑93 (MQ=255)

                                                gTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCaa                                         <  1:4887458/93‑1 (MQ=255)

                                                 ttGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAat                                        <  1:4302407/93‑1 (MQ=255)

                                                 ttGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAat                                        >  2:307711/1‑93 (MQ=255)

                                                 ttGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAat                                        <  1:21434/93‑1 (MQ=255)

                                                  tGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAata                                       <  2:3251290/93‑1 (MQ=255)

                                                    tCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAtt                                     <  2:1629836/93‑1 (MQ=255)

                                                    tCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAtt                                     <  2:2727329/93‑1 (MQ=255)

                                                    tCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAtt                                     <  1:238825/93‑1 (MQ=255)

                                                     cAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAttt                                    <  1:947163/93‑1 (MQ=255)

                                                      agatttCACAATTGCGCGATACAAACTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATGTCCGCACGCAATAtttg                                   >  1:3576850/4‑92 (MQ=255)

                                                       gCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCg                                  <  2:486541/93‑1 (MQ=255)

                                                       gCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCg                                  >  1:1356445/1‑93 (MQ=255)

                                                         tttCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAg                                <  2:1934185/93‑1 (MQ=255)

                                                          ttCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAgg                               <  2:5987411/93‑1 (MQ=255)

                                                          ttCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAgg                               <  2:1156287/93‑1 (MQ=255)

                                                           tCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGa                              <  1:4189300/93‑1 (MQ=255)

                                                            cacaATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGAt                             >  2:2329881/1‑93 (MQ=255)

                                                             acaATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATa                            <  2:3876495/93‑1 (MQ=255)

                                                               aaTTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAAt                          >  1:4933796/1‑93 (MQ=255)

                                                                aTTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATg                         >  1:3746231/1‑93 (MQ=255)

                                                                  tGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGc                       <  1:6106240/93‑1 (MQ=255)

                                                                   gcgcGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGcc                      <  2:5716376/93‑1 (MQ=255)

                                                                    cgcgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCct                     <  1:3691638/93‑1 (MQ=255)

                                                                     gcgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctc                    <  1:2272600/93‑1 (MQ=255)

                                                                      cgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctct                   <  1:4969427/93‑1 (MQ=255)

                                                                      cgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctct                   <  1:3933910/93‑1 (MQ=255)

                                                                      cgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctct                   <  1:2888168/93‑1 (MQ=255)

                                                                        aTACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCtcttc                 <  1:6035393/93‑1 (MQ=255)

                                                                         tACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCtcttct                >  1:5832748/1‑93 (MQ=255)

                                                                         tACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCtcttct                >  2:4965742/1‑93 (MQ=255)

                                                                          aCAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTcttctt               >  1:1498652/1‑93 (MQ=255)

                                                                          aCAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTcttctt               <  2:4116772/93‑1 (MQ=255)

                                                                          aCAAATTCACCAGCCCGGGCCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTcttctt               <  1:6341341/93‑1 (MQ=255)

                                                                            aaaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCa             >  2:2509450/1‑93 (MQ=255)

                                                                             aaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCAc            >  1:3604839/1‑93 (MQ=255)

                                                                             aaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCAc            <  1:5287422/93‑1 (MQ=255)

                                                                             aaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCAc            <  1:1787231/93‑1 (MQ=255)

                                                                             aaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCAc            >  1:100040/1‑93 (MQ=255)

                                                                                tCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAt         <  1:3352886/93‑1 (MQ=255)

                                                                                tCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAt         <  2:5814891/93‑1 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTCtt        >  1:2445957/1‑93 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAtt        >  1:4759884/1‑93 (MQ=255)

                                                                                  cccAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTc       <  1:405118/92‑1 (MQ=255)

                                                                                  aCCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTc       >  2:5091144/1‑93 (MQ=255)

                                                                                  aCCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTc       >  2:4036264/1‑93 (MQ=255)

                                                                                   ccAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCg      <  1:2847867/93‑1 (MQ=255)

                                                                                    cAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCgg     <  2:3530696/93‑1 (MQ=255)

                                                                                     aGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGgc    >  2:697916/1‑93 (MQ=255)

                                                                                     aGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGgc    <  1:1128656/93‑1 (MQ=255)

                                                                                      gCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGgcg   <  2:3621401/93‑1 (MQ=255)

CAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGC  >  NC_000913/2506778‑2506957

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: `‑`