BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,113,690	C→T	G371G (GGG→GGA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,690	0	C	T	96.5%	492.2 / 1.1	144	G371G (GGG→GGA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (4/1); new base T (74/65); total (78/66)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 148 total aligned reads displayed.

ACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCAT  >  NC_000913/4113600‑4113781

 ccAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCa                                                                                          >  2:1713756/1‑93 (MQ=255)

 ccAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCa                                                                                          >  2:4144899/1‑93 (MQ=255)

   agagTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACg                                                                                        >  2:4383630/1‑93 (MQ=255)

    gcgTCGGCCTCCCGCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGa                                                                                       <  2:227426/91‑1 (MQ=255)

      gTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGt                                                                                     >  1:900865/1‑93 (MQ=255)

        cGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCa                                                                                   <  2:1262721/93‑1 (MQ=255)

        cGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCa                                                                                   <  2:1325172/93‑1 (MQ=255)

         ggCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAg                                                                                  <  1:57585/93‑1 (MQ=255)

          gCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGa                                                                                 <  1:1345767/93‑1 (MQ=255)

          gCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGa                                                                                 <  2:4596469/93‑1 (MQ=255)

           ccTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGAc                                                                                <  2:3934480/93‑1 (MQ=255)

           ccTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGAc                                                                                <  1:3322651/93‑1 (MQ=255)

           ccTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGAc                                                                                >  1:2884085/1‑93 (MQ=255)

            cTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGAcc                                                                               <  2:53071/93‑1 (MQ=255)

             tGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCg                                                                              >  1:2061137/1‑93 (MQ=255)

              gCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGa                                                                             >  2:2676557/1‑93 (MQ=255)

                aTCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGaaa                                                                           >  2:2549123/1‑93 (MQ=255)

                aTCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGaaa                                                                           <  1:2306494/93‑1 (MQ=255)

                aTCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGaaa                                                                           <  1:2720924/93‑1 (MQ=255)

                aTCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGaaa                                                                           <  2:4100284/93‑1 (MQ=255)

                 tCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGaaaa                                                                          <  2:4712754/93‑1 (MQ=255)

                  cGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAAt                                                                         <  1:2027430/93‑1 (MQ=255)

                   gCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATc                                                                        <  2:2684753/93‑1 (MQ=255)

                     ttCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGc                                                                      >  2:4158346/1‑93 (MQ=255)

                          gCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCcgc                                                                 <  1:911736/93‑1 (MQ=255)

                          gCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCcgc                                                                 >  1:3369444/1‑93 (MQ=255)

                            ccGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCcgcgc                                                               <  2:1596957/92‑1 (MQ=255)

                            aCGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCcgcgc                                                               >  2:1912659/1‑93 (MQ=255)

                             cGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCcgcgcg                                                              >  1:1563131/1‑93 (MQ=255)

                              gacACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGCCAGACCGAAACTCGCCCCgcgcgc                                                             >  2:1699033/3‑93 (MQ=255)

                              gTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCgcgcgc                                                             >  1:22953/1‑93 (MQ=255)

                              gTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGAGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCgcgcgc                                                             >  1:1066154/1‑93 (MQ=255)

                               tCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCa                                                            >  1:1721213/1‑93 (MQ=255)

                               tCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCa                                                            >  1:2689114/1‑93 (MQ=255)

                                cACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATAGCCCCGCGCGCAt                                                           >  2:4312653/1‑93 (MQ=255)

                                 aCGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATa                                                          >  2:924734/1‑93 (MQ=255)

                                 aCGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATa                                                          <  2:3500147/93‑1 (MQ=255)

                                 aCGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATa                                                          >  1:2997749/1‑93 (MQ=255)

                                  cgcgTCTGCTAAGCAATAAACTCCAGCGTCGCACGTTTCATGTGGTTAGCGTTCACTCCACAAGTCGGACCGAAAATCACCACGCGCGCATAc                                                         <  1:4209832/93‑1 (MQ=255)

                                   gcgTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACg                                                        <  1:1790726/93‑1 (MQ=255)

                                    cgTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACgg                                                       >  1:2488469/1‑93 (MQ=255)

                                     gTCTGATAAGCAATAGCCGCCAGCGTCGCGCGTATAATTTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACggg                                                      <  1:3987335/93‑1 (MQ=255)

                                     gTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACggg                                                      >  1:4422748/1‑93 (MQ=255)

                                      tCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGt                                                     >  2:2808115/1‑93 (MQ=255)

                                       cTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTc                                                    <  2:2842559/93‑1 (MQ=255)

                                        tGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTcc                                                   >  2:2768752/1‑93 (MQ=255)

                                        tGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTcc                                                   >  2:4522035/1‑93 (MQ=255)

                                         gATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTccc                                                  <  2:3035046/93‑1 (MQ=255)

                                         gATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTccc                                                  >  2:1077302/1‑93 (MQ=255)

                                         gATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTccc                                                  <  1:776457/93‑1 (MQ=255)

                                          aTAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCa                                                 >  1:3741882/1‑93 (MQ=255)

                                          aTAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTAAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCa                                                 <  2:4127749/93‑1 (MQ=255)

                                            aaGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGt                                               >  1:3789193/1‑93 (MQ=255)

                                            aaGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGt                                               >  2:1646991/1‑93 (MQ=255)

                                             aGCAATAGACTCCAGGGCGCCCCGAATAATGTGGTTAGCGTTAACTCCCGGGGTCGGCCCGAAATTCGCCCCGCCCGCATACGGGCCCCAGtc                                              >  2:3758711/1‑92 (MQ=255)

                                             aGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGtc                                              >  1:3376209/1‑92 (MQ=255)

                                              gCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTAc                                             >  1:72909/1‑93 (MQ=255)

                                              gCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTAc                                             >  1:4337401/1‑93 (MQ=255)

                                              gCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTAc                                             <  2:3170595/93‑1 (MQ=255)

                                                aaTAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACgg                                           >  1:2036979/1‑93 (MQ=255)

                                                aaTAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACgg                                           <  2:709897/93‑1 (MQ=255)

                                                 aTAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGgc                                          >  1:3020985/1‑93 (MQ=255)

                                                  tAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGgcg                                         <  2:1878423/93‑1 (MQ=255)

                                                  tAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGgcg                                         <  1:4642160/93‑1 (MQ=255)

                                                        ccAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCa                                   >  1:3111052/1‑93 (MQ=255)

                                                        ccAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCa                                   <  1:3861723/93‑1 (MQ=255)

                                                         cAGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAg                                  >  1:1283062/1‑93 (MQ=255)

                                                          aGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGc                                 <  2:2348723/93‑1 (MQ=255)

                                                          aGCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGc                                 >  1:3173824/1‑93 (MQ=255)

                                                           gCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGcc                                <  2:2289615/93‑1 (MQ=255)

                                                           gCGTCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGcc                                >  2:850046/1‑93 (MQ=255)

                                                              tCGCGCGTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCgg                             <  1:3435695/93‑1 (MQ=255)

                                                                 cgcgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTaa                          >  1:4708272/1‑93 (MQ=255)

                                                                  gcgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTaaa                         <  1:4105696/93‑1 (MQ=255)

                                                                  gcgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTaaa                         <  2:3136027/93‑1 (MQ=255)

                                                                  gcgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTaaa                         >  2:1537400/1‑93 (MQ=255)

                                                                  gcgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTaaa                         >  2:3579682/1‑93 (MQ=255)

                                                                   cgTATAATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAAt                        <  2:3047648/93‑1 (MQ=255)

                                                                     tataATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGc                      >  2:1100472/1‑93 (MQ=255)

                                                                     tataATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGc                      >  1:2412712/1‑93 (MQ=255)

                                                                     tataATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGc                      >  2:3121021/1‑93 (MQ=255)

                                                                      ataATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGcc                     <  2:1787833/93‑1 (MQ=255)

                                                                      ataATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGcc                     >  1:1946271/1‑93 (MQ=255)

                                                                       taATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCg                    <  2:4714044/93‑1 (MQ=255)

                                                                       taATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCg                    <  2:3588642/93‑1 (MQ=255)

                                                                       taATGTGGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCg                    >  1:1733745/1‑93 (MQ=255)

                                                                       taATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCg                    <  1:4501707/93‑1 (MQ=255)

                                                                          tgtgGTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGaa                 >  1:4723227/1‑93 (MQ=255)

                                                                            tgGTTAGCGTTCACCCCACGAGGCAGACCGAAAATCGCCCCGCGCGCATACGGGTCACAGTACGGCGCACCGAGCCCGGTAAGTGCGGGAccc               >  1:1587643/1‑93 (MQ=255)

                                                                              gTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACcac             >  1:2904597/1‑93 (MQ=255)

                                                                              gTTAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGGCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACcac             >  2:514920/1‑93 (MQ=255)

                                                                                tAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACAt           <  1:1816058/93‑1 (MQ=255)

                                                                                tAGCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACAt           <  2:1899348/93‑1 (MQ=255)

                                                                                  gCGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATac         <  1:598888/93‑1 (MQ=255)

                                                                                   cGTTCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATaca        <  2:2843051/93‑1 (MQ=255)

                                                                                    gTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATacac       >  1:1799362/1‑93 (MQ=255)

                                                                                     ttCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATacaca      >  2:3942619/1‑93 (MQ=255)

                                                                                     ttCACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATacaca      >  2:744775/1‑93 (MQ=255)

                                                                                       cACTCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACAcc    <  1:4617556/93‑1 (MQ=255)

                                                                                         ccccACGAGTAAGACCGAAAAGCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCCCCCAGCCCGGTAAATGCCGGAACCACATACACACCAt  <  1:2166661/93‑1 (MQ=255)

ACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCAT  >  NC_000913/4113600‑4113781

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