BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,176	C→G	G182R (GGA→CGA)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,176	0	C	G	99.4%	557.7 / NA	174	G182R (GGA→CGA)	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base G (102/71); minor base A (1/0); total (103/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 184 total aligned reads displayed.

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑CG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGC  >  NC_000913/4013092‑4013257

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTa                                                                            >  1:43334/1‑93 (MQ=255)

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTa                                                                            >  1:2058771/1‑93 (MQ=255)

ttAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTa                                                                            >  2:4622970/1‑93 (MQ=255)

 tAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTaa                                                                           <  2:3498309/93‑1 (MQ=255)

  agcagcCCGCAGCGCCCGGCGAGTAGTTTCAGCGCTCTCCTGCGGACGGCTGTTACCATGCCCTCCATATAAAGCGAGAATC‑GGCGCGTTAAg                                                                          <  2:2836981/93‑1 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGa                                                                         >  2:2088006/1‑93 (MQ=255)

   gcagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGa                                                                         >  2:696686/1‑93 (MQ=255)

    cagcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGAt                                                                        <  2:1933453/93‑1 (MQ=255)

     agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATc                                                                       <  2:1573919/93‑1 (MQ=255)

     agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATc                                                                       <  2:1213816/93‑1 (MQ=255)

     agcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATc                                                                       >  1:811078/1‑93 (MQ=255)

      gcCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCg                                                                      >  1:4265367/1‑93 (MQ=255)

        ccGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGt                                                                    >  1:2985391/1‑93 (MQ=255)

         cGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGtt                                                                   <  2:4575231/93‑1 (MQ=255)

         cGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGtt                                                                   <  1:948376/93‑1 (MQ=255)

          gCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTg                                                                  >  1:4453177/1‑93 (MQ=255)

           cAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGc                                                                 <  1:143447/93‑1 (MQ=255)

            aGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCg                                                                <  1:164701/93‑1 (MQ=255)

               gcCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGata                                                             >  1:646924/1‑93 (MQ=255)

                tcTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGatat                                                            >  1:4520237/2‑93 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATc                                                           <  1:2331732/93‑1 (MQ=255)

                 cTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGGTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATc                                                           >  1:3604159/1‑93 (MQ=255)

                  tGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCa                                                          >  2:1776251/1‑93 (MQ=255)

                     cgcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACa                                                       <  2:2309136/93‑1 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAg                                                      >  1:2754955/1‑93 (MQ=255)

                      gcATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAg                                                      >  2:2669978/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATCAGGCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGctt                                                   >  2:164760/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                   >  2:851336/1‑93 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                   <  1:201932/93‑1 (MQ=255)

                         tGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                   <  1:2414748/93‑1 (MQ=255)

                          ggTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttg                                                  >  1:4578140/1‑93 (MQ=255)

                          ggTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttg                                                  >  1:3908041/1‑93 (MQ=255)

                           gTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgt                                                 <  1:554750/93‑1 (MQ=255)

                           gTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgt                                                 <  1:1843972/93‑1 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                >  1:1094689/1‑93 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                <  2:4622424/93‑1 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                <  2:4259141/93‑1 (MQ=255)

                            tttCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                >  1:1356452/1‑93 (MQ=255)

                            tttCAACGCTCTCCTGAGTAATTTTGTTATCCTGACCACCATATAGGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGgttgtt                                                <  2:3461312/93‑1 (MQ=255)

                              tCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                              <  1:3991853/93‑1 (MQ=255)

                              tCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTc                                              <  2:3181629/93‑1 (MQ=255)

                                aaCGCTCTCGTGCGGAATGCTGTAATCCTGACCACCATGTAAGCCGAGAATC‑GGCGCGTAAAGATCGGTTGCGATATCAACAGGTTGTTTCgg                                            <  2:233940/93‑1 (MQ=255)

                                aaCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGGTGTTTCgg                                            >  1:862411/1‑93 (MQ=255)

                                  cGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTg                                          <  2:3027202/93‑1 (MQ=255)

                                    ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                        >  1:387091/1‑93 (MQ=255)

                                    ctctCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGaa                                        >  2:1692960/1‑93 (MQ=255)

                                     tctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                                       >  2:3519333/1‑93 (MQ=255)

                                      ctcCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                      <  2:3088112/93‑1 (MQ=255)

                                        ccTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCa                                    <  1:240560/93‑1 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAGGATCGGTTGCGATATAAACAGGCTGCTTCGGTGAAGTCAGc                                  >  1:199938/1‑93 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                  <  1:524095/93‑1 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                  >  2:2007055/1‑93 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                  <  1:988448/93‑1 (MQ=255)

                                          tGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTAGCTTCGGGGAATTCaca                                  >  2:2007054/1‑91 (MQ=255)

                                           gCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                 >  1:416932/1‑93 (MQ=255)

                                           gCGGAATGCTGTTATCCTGACCACCATATAAGACGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                 >  2:2508905/1‑93 (MQ=255)

                                             ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                               >  1:122320/1‑93 (MQ=255)

                                             ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                               >  1:3575286/1‑93 (MQ=255)

                                             ggAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                               <  1:4733095/93‑1 (MQ=255)

                                              gAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACt                              <  1:1000710/93‑1 (MQ=255)

                                              gAATGCTGTTAGCCTGACCACCATATGAGCCGAGAATC‑GGCTCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACt                              <  2:3293931/93‑1 (MQ=255)

                                               aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                             <  1:2776187/93‑1 (MQ=255)

                                               aaTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACtt                             <  1:3432981/93‑1 (MQ=255)

                                                aTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTg                            >  2:587417/1‑93 (MQ=255)

                                                aTGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTg                            <  2:3521220/93‑1 (MQ=255)

                                                  gCTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                          <  1:747794/93‑1 (MQ=255)

                                                   cTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                         >  1:2263732/1‑93 (MQ=255)

                                                   cTGTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                         >  2:2982143/1‑93 (MQ=255)

                                                     gTTATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGcc                       >  1:3962863/1‑93 (MQ=255)

                                                       tATCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg                     <  1:1733856/93‑1 (MQ=255)

                                                        aTCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGt                    >  2:4492664/1‑93 (MQ=255)

                                                        aTCCTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGt                    >  2:3730226/1‑93 (MQ=255)

                                                         tCCTGACCACCATATAAGGCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTc                   >  1:4412281/1‑93 (MQ=255)

                                                           cTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAg                 >  2:2157858/1‑93 (MQ=255)

                                                           cTGACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAg                 >  2:2233822/1‑93 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGtt               <  1:741438/93‑1 (MQ=255)

                                                             gACCACCATATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACCTGTCGCCCGTCAGtt               >  1:3354854/1‑93 (MQ=255)

                                                               ccaccaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTg             >  1:104069/1‑93 (MQ=255)

                                                                caccaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGc            >  1:4647018/1‑93 (MQ=255)

                                                                caccaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGc            >  1:3224799/1‑93 (MQ=255)

                                                                 accaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGcc           >  1:4548702/1‑93 (MQ=255)

                                                                 accaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGcc           >  2:639055/1‑93 (MQ=255)

                                                                  ccaTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCg          >  2:626305/1‑93 (MQ=255)

                                                                   caTATAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGt         <  2:364768/93‑1 (MQ=255)

                                                                    atatAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa        >  1:4264561/1‑93 (MQ=255)

                                                                    atatAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa        >  1:2489535/1‑93 (MQ=255)

                                                                    atatAAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa        >  1:2064611/1‑93 (MQ=255)

                                                                     tataAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAc       >  2:2430897/1‑93 (MQ=255)

                                                                     tataAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAc       >  2:3451792/1‑93 (MQ=255)

                                                                     tataAGCCGAGAATA‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAc       <  2:539064/93‑1 (MQ=255)

                                                                      ataAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc      >  1:1231011/1‑93 (MQ=255)

                                                                      ataAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTAcc      >  2:4101841/1‑93 (MQ=255)

                                                                       taAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa     <  1:381877/93‑1 (MQ=255)

                                                                       taAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa     <  1:1979709/93‑1 (MQ=255)

                                                                       taAGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCa     <  1:3299456/93‑1 (MQ=255)

                                                                        aaGCCGAGAATCGAC‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCAc    >  1:3136086/1‑93 (MQ=255)

                                                                        aaGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCAc    >  1:1540175/1‑93 (MQ=255)

                                                                         aGCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACg   >  2:3999826/1‑93 (MQ=255)

                                                                          gCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGc  <  2:3793892/93‑1 (MQ=255)

                                                                          gCCGAGAATC‑GGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGc  >  2:3579705/1‑93 (MQ=255)

TTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATC‑CG‑GCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGC  >  NC_000913/4013092‑4013257

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