BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,019,825	G→A	P163P (CCG→CCA)	b2881 →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,019,825	0	G	A	99.0%	722.6 / ‑4.3	205	P163P (CCG→CCA)	b2881	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (74/129); total (74/131)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 214 total aligned reads displayed.

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

ccGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAAc                                                                                         <  1:2672722/93‑1 (MQ=255)

  gcgATGCTGGCTAGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATTATCCGCAGGCCACCATTCATATCGCTCCAACtt                                                                                       >  2:3133075/1‑93 (MQ=255)

    gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTc                                                                                     <  1:3811624/93‑1 (MQ=255)

    caTGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTc                                                                                     <  2:2134273/92‑1 (MQ=255)

     aTGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTcc                                                                                    >  1:299570/1‑93 (MQ=255)

      tGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGGCACCATTGATATCGCTCCAACTTTCCg                                                                                   >  2:250184/1‑93 (MQ=255)

      tGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCg                                                                                   >  2:940301/1‑93 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  2:1271418/93‑1 (MQ=255)

       gctggctgGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGt                                                                                  <  1:4271858/93‑1 (MQ=255)

        ctggctggCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCAACCATTGATATCGCTCCAACTTTCCGTg                                                                                 >  2:4305286/1‑93 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                <  1:3561974/93‑1 (MQ=255)

         tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTga                                                                                >  1:3922866/1‑93 (MQ=255)

          ggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgac                                                                               <  1:3586288/93‑1 (MQ=255)

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGGCACCATTGATATCGCTCCAACTTTCCGTgacgac                                                                            >  2:1750909/1‑93 (MQ=255)

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacgac                                                                            >  1:20937/1‑93 (MQ=255)

             tggcAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTgacgac                                                                            >  1:3863280/1‑93 (MQ=255)

                cagcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTa                                                                         <  1:3674208/93‑1 (MQ=255)

                  gcaATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGAACTAGa                                                                       >  2:3212055/1‑93 (MQ=255)

                   caATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGaa                                                                      <  1:2907667/93‑1 (MQ=255)

                    aaTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAg                                                                     <  1:1834036/93‑1 (MQ=255)

                    aaTACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAg                                                                     <  1:1529337/93‑1 (MQ=255)

                    aaTACTATCAGGTCATTGAACTGGCGGTTCCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAg                                                                     <  2:645994/93‑1 (MQ=255)

                     aTACTATCAGGTCATTGAACTGGCGGTTGCACGGAAAAATAATACGCAGGCCACCATTGATATCGCTCCAACTTTACGTGACGACCTAGAAgc                                                                    >  2:783571/1‑92 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   >  1:577135/1‑93 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   <  1:4154283/93‑1 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   <  1:3714397/93‑1 (MQ=255)

                      tACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTc                                                                   >  2:1922816/1‑93 (MQ=255)

                        cTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAt                                                                 >  1:3057080/1‑93 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  2:3545551/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                >  2:2289977/1‑93 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  2:398037/93‑1 (MQ=255)

                         tATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCAtt                                                                <  1:2723872/93‑1 (MQ=255)

                          aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTg                                                               >  2:136415/1‑93 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  2:4368537/93‑1 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           >  2:1329869/1‑93 (MQ=255)

                              ggTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTaa                                                           <  1:2827131/93‑1 (MQ=255)

                                 cATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCa                                                        <  1:3300667/93‑1 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       <  2:2759045/93‑1 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       <  2:1687547/93‑1 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       >  1:1034439/1‑93 (MQ=255)

                                  aTTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       >  2:4543784/1‑93 (MQ=255)

                                  aTTGAACTGGCAGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAt                                                       <  1:4681314/93‑1 (MQ=255)

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      >  1:1141770/1‑93 (MQ=255)

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      <  1:2558836/93‑1 (MQ=255)

                                   ttGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCAtt                                                      >  1:806457/1‑93 (MQ=255)

                                     gAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  2:2723995/93‑1 (MQ=255)

                                     gAACTGGCGCTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTAt                                                    <  2:2641754/93‑1 (MQ=255)

                                       aCTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATcc                                                  >  2:577624/1‑93 (MQ=255)

                                        cTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCt                                                 <  2:4093643/93‑1 (MQ=255)

                                              gTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACt                                           >  2:3369429/1‑93 (MQ=255)

                                                tGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGa                                         <  2:202461/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        <  1:639663/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        <  1:1974956/93‑1 (MQ=255)

                                                 gCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGAt                                        >  2:890443/1‑93 (MQ=255)

                                                   aCGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGc                                      <  2:3244488/93‑1 (MQ=255)

                                                      tAAAAATAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCcgc                                   >  1:3397591/1‑93 (MQ=255)

                                                        aaaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGa                                 <  2:4444349/93‑1 (MQ=255)

                                                         aaaTAATCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaa                                <  2:1984542/93‑1 (MQ=255)

                                                           ataatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGaaaa                              >  2:506857/1‑93 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:2889140/93‑1 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:1496187/93‑1 (MQ=255)

                                                             aatCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATg                            <  2:1635716/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  2:3813175/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  1:1920213/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  2:2337993/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           <  1:1167110/93‑1 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           >  1:3772556/1‑93 (MQ=255)

                                                              atCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATgg                           >  1:484077/1‑93 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGt                          <  1:443492/93‑1 (MQ=255)

                                                               tCCGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGt                          >  2:1742301/1‑93 (MQ=255)

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGc                        >  1:4204422/1‑93 (MQ=255)

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGc                        <  1:555044/93‑1 (MQ=255)

                                                                 cGCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGc                        <  1:93344/93‑1 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       <  2:4231831/93‑1 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       >  2:4607884/1‑93 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       >  2:1614803/1‑93 (MQ=255)

                                                                  gCAGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCa                       >  1:2763120/1‑93 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     <  2:1859184/93‑1 (MQ=255)

                                                                    aGGCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAgg                     <  1:1668172/93‑1 (MQ=255)

                                                                     ggCCACCATTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGc                    >  2:2974350/1‑93 (MQ=255)

                                                                     ggCCACCATCGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGc                    <  2:257045/93‑1 (MQ=255)

                                                                         cacaTTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                <  2:3506631/91‑1 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                >  1:2413164/1‑93 (MQ=255)

                                                                         accaTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGaaa                <  1:3907216/93‑1 (MQ=255)

                                                                           caTTGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAAcc              <  1:615177/93‑1 (MQ=255)

                                                                             ttGATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCt            <  2:2409720/93‑1 (MQ=255)

                                                                               gATATCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGc          <  2:3714716/93‑1 (MQ=255)

                                                                                atatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCt         >  1:3208367/1‑93 (MQ=255)

                                                                                 tatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTa        >  1:4460481/1‑93 (MQ=255)

                                                                                 tatCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTa        <  2:3669907/93‑1 (MQ=255)

                                                                                  atCGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTAt       <  2:742018/93‑1 (MQ=255)

                                                                                    cGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGt     <  1:4494989/93‑1 (MQ=255)

                                                                                    cGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGt     <  1:4394749/93‑1 (MQ=255)

                                                                                    cGCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGt     <  1:2974740/93‑1 (MQ=255)

                                                                                     gCTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGtt    >  1:2334607/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   >  2:1763598/1‑93 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   <  1:4367535/93‑1 (MQ=255)

                                                                                      cTCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTg   <  1:133063/93‑1 (MQ=255)

                                                                                       tCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGa  <  2:1448902/93‑1 (MQ=255)

                                                                                       tCCAACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGa  >  1:2199654/1‑93 (MQ=255)

CCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTTCCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGA  >  NC_000913/3019735‑3019914

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