BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,071,652	A→G	S114G (AGC→GGC)	yihV →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,071,652	0	A	G	86.6%	738.1 / 103.8	306	S114G (AGC→GGC)	yihV	putative kinase
Reads supporting (aligned to +/- strand): ref base A (17/24); new base G (147/118); total (164/142)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 316 total aligned reads displayed.

CCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACG  >  NC_000913/4071560‑4071742

cccGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCa                                                                                            >  1:1829776/1‑93 (MQ=255)

  cGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGc                                                                                          <  2:1684265/93‑1 (MQ=255)

   gTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGcc                                                                                         >  2:1545806/1‑93 (MQ=255)

    ttACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGccc                                                                                        <  2:5025779/93‑1 (MQ=255)

    ttACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGccc                                                                                        >  2:4626138/1‑93 (MQ=255)

     tACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCg                                                                                       >  1:2026922/1‑93 (MQ=255)

       cacCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGa                                                                                     <  1:6323662/93‑1 (MQ=255)

        acCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGAt                                                                                    >  2:5312690/1‑93 (MQ=255)

        acCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGAt                                                                                    <  1:2126645/93‑1 (MQ=255)

         ccAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATc                                                                                   >  1:446809/1‑93 (MQ=255)

          cAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATct                                                                                  >  2:6887645/1‑93 (MQ=255)

           aaaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctg                                                                                 >  1:2252601/1‑93 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctgc                                                                                <  2:1830641/93‑1 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctgc                                                                                <  1:4781628/93‑1 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctgc                                                                                >  1:6087723/1‑93 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctgc                                                                                <  1:125329/93‑1 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATctgc                                                                                <  1:5317080/93‑1 (MQ=255)

            aaCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATctgc                                                                                <  2:804307/93‑1 (MQ=255)

             aCGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATctgct                                                                               <  2:7309474/93‑1 (MQ=255)

               ggTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCtgctgc                                                                             >  2:5198777/1‑93 (MQ=255)

                gTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGcc                                                                            >  1:5785063/1‑93 (MQ=255)

                 tataACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCt                                                                           >  2:4053951/1‑93 (MQ=255)

                 tataACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCt                                                                           >  1:6479285/1‑93 (MQ=255)

                    aaCCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGAc                                                                        >  1:2810996/1‑93 (MQ=255)

                    aaCCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGAc                                                                        <  1:1318564/93‑1 (MQ=255)

                     aCCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACg                                                                       >  2:6588534/1‑93 (MQ=255)

                      ccAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGc                                                                      >  1:761984/1‑93 (MQ=255)

                          gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCagag                                                                  >  2:6824952/1‑93 (MQ=255)

                          gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCagag                                                                  <  2:2571775/93‑1 (MQ=255)

                          gCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCagag                                                                  <  2:2191040/93‑1 (MQ=255)

                            gAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTg                                                                >  2:3305939/1‑93 (MQ=255)

                             aaaTCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTgg                                                               <  1:72930/93‑1 (MQ=255)

                             aaaTCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTgg                                                               >  2:4001009/1‑93 (MQ=255)

                              aaTCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGt                                                              >  1:405550/1‑93 (MQ=255)

                              aaTCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGt                                                              >  2:4769621/1‑93 (MQ=255)

                               aTCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGtt                                                             >  2:4857483/1‑93 (MQ=255)

                                tCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTg                                                            <  1:6681204/93‑1 (MQ=255)

                                  ttCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTgga                                                          <  1:5175615/93‑1 (MQ=255)

                                  ttCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTgga                                                          <  1:6198341/93‑1 (MQ=255)

                                   tCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTTCCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                         <  1:7361678/93‑1 (MQ=255)

                                   tCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                         <  2:3734320/93‑1 (MQ=255)

                                   tCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggag                                                         >  2:5960058/1‑93 (MQ=255)

                                    cGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggagg                                                        <  1:684397/93‑1 (MQ=255)

                                     gCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTggagga                                                       >  2:4388877/1‑93 (MQ=255)

                                       aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGaaa                                                     <  1:2156694/93‑1 (MQ=255)

                                       aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGaaa                                                     >  1:1637201/1‑93 (MQ=255)

                                       aaTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGaaa                                                     >  2:5654249/1‑93 (MQ=255)

                                        aTCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAAt                                                    >  2:2496776/1‑93 (MQ=255)

                                         tCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAAtt                                                   >  1:3786614/1‑93 (MQ=255)

                                           cGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                 >  2:3682447/1‑93 (MQ=255)

                                           cGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                 <  2:3053355/93‑1 (MQ=255)

                                           cGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGa                                                 <  1:3567482/93‑1 (MQ=255)

                                            gCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                >  1:1015336/1‑93 (MQ=255)

                                            gCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                >  2:2299003/1‑93 (MQ=255)

                                            gCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAt                                                >  1:3490233/1‑93 (MQ=255)

                                              catcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAttt                                              <  1:2894929/93‑1 (MQ=255)

                                              catcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAttt                                              >  1:3083136/1‑93 (MQ=255)

                                              catcatGGGGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGAttt                                              >  1:1513850/1‑93 (MQ=255)

                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtc                                             >  2:4883960/1‑93 (MQ=255)

                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtc                                             >  2:5592562/1‑93 (MQ=255)

                                               atcatGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTtc                                             <  2:2691242/93‑1 (MQ=255)

                                                    ggtggATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCa                                        <  2:3378926/93‑1 (MQ=255)

                                                    ggtggATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCa                                        <  1:94137/93‑1 (MQ=255)

                                                        gATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTgg                                    <  2:4683060/93‑1 (MQ=255)

                                                        gATACCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTgg                                    <  1:5974206/93‑1 (MQ=255)

                                                           aCCAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAt                                 <  1:3390130/93‑1 (MQ=255)

                                                             cAAAGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgt                               >  2:4145028/1‑93 (MQ=255)

                                                              aaaGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                              >  2:6954393/1‑93 (MQ=255)

                                                              aaaGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                              <  2:2093191/93‑1 (MQ=255)

                                                              aaaGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgtt                              >  1:3668868/1‑93 (MQ=255)

                                                               aaGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttg                             <  1:3950551/93‑1 (MQ=255)

                                                                aGGCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttgt                            <  2:5937794/93‑1 (MQ=255)

                                                                 ggCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTtgtg                           >  1:7013801/1‑93 (MQ=255)

                                                                  gCGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGc                          <  2:7505385/93‑1 (MQ=255)

                                                                   cGAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCt                         <  2:1212389/93‑1 (MQ=255)

                                                                    gAGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTg                        >  1:5598625/1‑93 (MQ=255)

                                                                     aGCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                       <  1:459361/93‑1 (MQ=255)

                                                                      gCGGATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                      <  1:2339962/93‑1 (MQ=255)

                                                                         gATAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGa                   <  2:1037643/93‑1 (MQ=255)

                                                                           tAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATg                 >  2:6173566/1‑93 (MQ=255)

                                                                           tAATCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATg                 >  1:2381393/1‑93 (MQ=255)

                                                                           tAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATg                 >  1:485422/1‑93 (MQ=255)

                                                                              tCATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTAc              >  2:127196/1‑93 (MQ=255)

                                                                               cATTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACg             >  2:3103172/1‑93 (MQ=255)

                                                                                aTTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc            <  1:5491850/93‑1 (MQ=255)

                                                                                aTTAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGc            >  1:4497485/1‑93 (MQ=255)

                                                                                 ttAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCt           <  2:6327410/93‑1 (MQ=255)

                                                                                 ttAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCt           <  2:1468157/93‑1 (MQ=255)

                                                                                 ttAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCt           >  1:1772524/1‑93 (MQ=255)

                                                                                 ttAACTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCt           <  1:5863338/93‑1 (MQ=255)

                                                                                 ttAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCt           >  1:4241297/1‑93 (MQ=255)

                                                                                    aCTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGc        <  2:6718553/93‑1 (MQ=255)

                                                                                    aCTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGc        <  2:2776342/93‑1 (MQ=255)

                                                                                     cTACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCa       >  1:3211773/1‑93 (MQ=255)

                                                                                     cTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCa       <  1:2568197/93‑1 (MQ=255)

                                                                                      tACTCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCac      <  2:6085716/93‑1 (MQ=255)

                                                                                      tACCCCGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCac      >  2:636330/1‑93 (MQ=255)

                                                                                      tACCCAGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCac      <  2:5062769/93‑1 (MQ=255)

                                                                                        ccccGGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacga    <  2:234193/93‑1 (MQ=255)

                                                                                         cccggcccggATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCacgac   <  1:2495422/93‑1 (MQ=255)

CCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACG  >  NC_000913/4071560‑4071742

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