BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	908,793	A→G	C494R (TGC→CGC)	poxB ←	pyruvate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	908,793	0	A	G	99.5%	608.8 / ‑5.1	196	C494R (TGC→CGC)	poxB	pyruvate oxidase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (89/105); total (89/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 206 total aligned reads displayed.

CCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGC  >  NC_000913/908703‑908882

ccaccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCcgcgcg                                                                                         <  1:3241282/93‑1 (MQ=255)

   ccaACACCGGACCGTCGATGGAGAAGGGGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCtt                                                                                      >  2:3455837/1‑93 (MQ=255)

   ccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCtt                                                                                      <  1:4533953/93‑1 (MQ=255)

   ccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCtt                                                                                      >  1:7503701/1‑93 (MQ=255)

   ccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCtt                                                                                      <  1:2223652/93‑1 (MQ=255)

       cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGATTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCgga                                                                                  >  2:7144116/1‑92 (MQ=255)

        cccGGACCGGCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCa                                                                                 <  2:1790428/92‑1 (MQ=255)

        acCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCa                                                                                 >  2:1983179/1‑93 (MQ=255)

        acCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCa                                                                                 >  2:4783757/1‑93 (MQ=255)

         ccGGACCGTCGATGGAGAAGTCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                                                                                <  2:4686235/93‑1 (MQ=255)

         ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                                                                                <  1:2394822/93‑1 (MQ=255)

         ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                                                                                <  1:6029522/93‑1 (MQ=255)

         ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                                                                                <  1:5284388/93‑1 (MQ=255)

         ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCaa                                                                                <  1:870922/93‑1 (MQ=255)

          cGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAAt                                                                               <  1:6183831/93‑1 (MQ=255)

          cGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAAt                                                                               <  1:4201209/93‑1 (MQ=255)

           ggACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATg                                                                              >  2:5949779/1‑93 (MQ=255)

            gACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGc                                                                             <  1:1834225/93‑1 (MQ=255)

            gACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGc                                                                             >  1:4580075/1‑93 (MQ=255)

             gccGTCGATGGAGAAGGGGCGTTGCAGCGCTTCATCAACTGCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCg                                                                            <  1:6552918/92‑1 (MQ=255)

             aCCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCg                                                                            <  2:3369475/93‑1 (MQ=255)

             aCCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCg                                                                            >  2:7781413/1‑93 (MQ=255)

              ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCgg                                                                           <  1:7658711/93‑1 (MQ=255)

              ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCgg                                                                           <  1:7828541/93‑1 (MQ=255)

                   gATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCaaa                                                                      >  1:1728603/1‑93 (MQ=255)

                      ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt                                                                   <  1:6162696/93‑1 (MQ=255)

                      ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt                                                                   <  1:6989489/93‑1 (MQ=255)

                      ggAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGtt                                                                   <  2:3310301/93‑1 (MQ=255)

                       gagaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGttt                                                                  >  2:1193018/1‑93 (MQ=255)

                        agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTtg                                                                 <  2:6993379/93‑1 (MQ=255)

                             gcgcgcTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCg                                                            <  2:5276807/93‑1 (MQ=255)

                              cgcgCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATCCCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgt                                                           <  2:4520795/93‑1 (MQ=255)

                              cgcgCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgt                                                           <  1:5883992/93‑1 (MQ=255)

                                cgctGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgtgt                                                         >  2:3881822/4‑93 (MQ=255)

                                cgctGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgtgt                                                         <  1:2738130/90‑1 (MQ=255)

                                cgctGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCgtgt                                                         <  1:5422337/90‑1 (MQ=255)

                                 gctGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTa                                                        >  1:6217666/3‑93 (MQ=255)

                                 gctGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTa                                                        <  2:1178966/91‑1 (MQ=255)

                                 gctGCAGGGCTTCATCAACTTCAGACGCTTCTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATAAGGGCAAAGTTTGTGTCGTGTa                                                        >  2:2237799/3‑93 (MQ=255)

                                   tGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGt                                                      <  2:3995201/93‑1 (MQ=255)

                                    gCAGGGCTTCATCATCTTCAGACGTTTTTTCTACTTGGATACCCGTTATGCCGCGCGCATCGGTATTGCGGGCAAAGTCTGTGTCGTGTAgat                                                     >  1:7016209/1‑91 (MQ=255)

                                    gCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGtt                                                     >  1:2906241/1‑93 (MQ=255)

                                    gCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGtt                                                     >  1:5370980/1‑93 (MQ=255)

                                     cAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTc                                                    <  2:7258378/93‑1 (MQ=255)

                                     cAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTc                                                    >  2:3195840/1‑93 (MQ=255)

                                      aGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCg                                                   >  1:6652923/1‑93 (MQ=255)

                                      aGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCg                                                   <  1:3743611/93‑1 (MQ=255)

                                      aGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCg                                                   >  2:4922054/1‑93 (MQ=255)

                                       gggCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCgg                                                  <  1:6500558/93‑1 (MQ=255)

                                       gggCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCgg                                                  >  2:2133269/1‑93 (MQ=255)

                                       gggCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCgg                                                  >  1:3767227/1‑93 (MQ=255)

                                        ggCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGt                                                 >  1:3400970/1‑93 (MQ=255)

                                           ttCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGcc                                              >  2:1340957/1‑93 (MQ=255)

                                             catcaACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgt                                            >  2:4654457/1‑93 (MQ=255)

                                             catcaACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgt                                            <  1:4240068/93‑1 (MQ=255)

                                              atcaACTTCAGACGCTTTTTCTACACGGATGCCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtc                                           <  1:479488/93‑1 (MQ=255)

                                               tcaACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtca                                          <  2:1417291/93‑1 (MQ=255)

                                                caACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcag                                         >  1:5043000/1‑93 (MQ=255)

                                                caACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcag                                         >  1:4744100/1‑93 (MQ=255)

                                                caACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcag                                         <  1:7781414/93‑1 (MQ=255)

                                                 aaCTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcagt                                        <  2:6772069/93‑1 (MQ=255)

                                                 aaCTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcagt                                        >  2:161041/1‑93 (MQ=255)

                                                  aCTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcagtc                                       >  2:286938/1‑93 (MQ=255)

                                                  aCTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcagtc                                       >  2:5626298/1‑93 (MQ=255)

                                                   cTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCTCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCgtcagtca                                      <  1:7593025/93‑1 (MQ=255)

                                                    ttCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCaa                                     >  2:3789376/1‑93 (MQ=255)

                                                     tCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCaaa                                    <  2:7358256/93‑1 (MQ=255)

                                                     tCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCaaa                                    <  2:895903/93‑1 (MQ=255)

                                                     tCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCaaa                                    <  1:7221894/93‑1 (MQ=255)

                                                     tCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCaaa                                    >  1:2339713/1‑93 (MQ=255)

                                                      cAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAAt                                   >  1:3653892/1‑93 (MQ=255)

                                                      cAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAAt                                   <  2:648954/93‑1 (MQ=255)

                                                      cAGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAAt                                   <  1:6299107/93‑1 (MQ=255)

                                                       aGACGCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATa                                  <  1:5082657/93‑1 (MQ=255)

                                                           gCTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGcca                              >  1:710068/1‑93 (MQ=255)

                                                            cTTTTTCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCAGGCAAAGTTTGTGTAGTGTAGTTCGGTGCCGTCAGTCAAATAGccac                             <  1:442746/93‑1 (MQ=255)

                                                             tttttCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGccacc                            <  1:5506137/93‑1 (MQ=255)

                                                             cttttCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGccacc                            <  1:6341936/92‑1 (MQ=255)

                                                              ttttCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCACTCAAATAGccacca                           <  1:3525871/93‑1 (MQ=255)

                                                                ttCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGc                         >  1:3024981/1‑93 (MQ=255)

                                                                ttCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGCAGTTCGGTGCCGTCAGTCAAATAGCCACCAGc                         <  2:6582545/93‑1 (MQ=255)

                                                                 tCTACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCt                        <  1:1288179/93‑1 (MQ=255)

                                                                   tACACGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCttt                      >  1:4560976/1‑93 (MQ=255)

                                                                      acGGATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAt                   >  2:3260049/1‑93 (MQ=255)

                                                                        ggATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAtct                 >  2:4568101/1‑93 (MQ=255)

                                                                        ggATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAtct                 >  1:1178751/1‑93 (MQ=255)

                                                                        ggATACCCGTAATGCCGCGCGCTTCGGCAATGCGCGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAtct                 <  1:437100/93‑1 (MQ=255)

                                                                         gATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAtctc                >  1:3177021/1‑93 (MQ=255)

                                                                         gATACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCAtctc                >  1:4099465/1‑93 (MQ=255)

                                                                           tACCCGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCa              >  1:4687075/1‑93 (MQ=255)

                                                                             cccGTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGGTTGTGTCGTGTAGGTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATc            >  2:6675188/1‑93 (MQ=255)

                                                                                gTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGcc         >  2:1842093/1‑93 (MQ=255)

                                                                                gTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGcc         >  2:3080540/1‑93 (MQ=255)

                                                                                gTAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGcc         >  2:1622822/1‑93 (MQ=255)

                                                                                 tAATGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCca        <  1:1533717/93‑1 (MQ=255)

                                                                                    tGCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACaa     >  1:1243361/1‑93 (MQ=255)

                                                                                     gccgccCGCTGCGGGAATGCGGTCAAAGTTTTGGACGTGTAGTTCGGTTCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACaaa    <  1:3956222/93‑1 (MQ=255)

                                                                                     gCCGCGCGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACaaa    <  1:2035874/93‑1 (MQ=255)

                                                                                      gcgcgcgCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAg   <  1:3647240/92‑1 (MQ=255)

                                                                                       cgcgcgCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGc  <  1:4268479/93‑1 (MQ=255)

CCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTGTAGTTCGGTGCCGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGC  >  NC_000913/908703‑908882

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