BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	979,388	T→C	R1280R (CGT→CGC)	mukB →	MukB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	979,388	0	T	C	99.4%	531.1 / ‑5.1	165	R1280R (CGT→CGC)	mukB	MukB
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (75/89); total (76/89)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 174 total aligned reads displayed.

GAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCGGA  >  NC_000913/979299‑979477

  gCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAac                                                                                      <  1:3094150/93‑1 (MQ=255)

    aGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgc                                                                                    >  2:4421576/1‑93 (MQ=255)

    aGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgc                                                                                    >  1:1250194/1‑93 (MQ=255)

    aGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgc                                                                                    <  2:4150013/93‑1 (MQ=255)

    aGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgc                                                                                    <  2:644971/93‑1 (MQ=255)

      aaCCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgcac                                                                                  <  1:4533150/93‑1 (MQ=255)

        ccgtatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAacgcacgc                                                                                <  1:4114311/93‑1 (MQ=255)

          gtatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCa                                                                              <  2:482021/93‑1 (MQ=255)

           tatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCAt                                                                             <  1:1420021/93‑1 (MQ=255)

           tatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCAt                                                                             <  1:2218724/93‑1 (MQ=255)

           tatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCAt                                                                             <  2:2610234/93‑1 (MQ=255)

           tatccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCAt                                                                             >  2:1913361/1‑93 (MQ=255)

            atccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATg                                                                            >  2:202413/1‑93 (MQ=255)

             tccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGc                                                                           >  1:144341/1‑93 (MQ=255)

             tccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGc                                                                           >  1:2223734/1‑93 (MQ=255)

             tccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGc                                                                           >  1:2223736/1‑93 (MQ=255)

             tccgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGc                                                                           >  2:4039503/1‑93 (MQ=255)

               cgtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTa                                                                         <  1:1414720/93‑1 (MQ=255)

                gtatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTAc                                                                        <  2:855380/93‑1 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACt                                                                       >  1:2949924/1‑93 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACt                                                                       <  1:736894/93‑1 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACt                                                                       >  1:4369467/1‑93 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACt                                                                       >  2:3511101/1‑93 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACt                                                                       <  2:3851363/93‑1 (MQ=255)

                 tatGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACCCACGCCATGCTACt                                                                       >  1:280447/1‑93 (MQ=255)

                  atGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTTCTg                                                                      <  1:3356446/93‑1 (MQ=255)

                  atGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTg                                                                      >  2:3112103/1‑93 (MQ=255)

                     cTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAt                                                                   >  1:1296174/1‑93 (MQ=255)

                     cTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAt                                                                   <  1:2238521/93‑1 (MQ=255)

                      tCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAtg                                                                  <  1:4682580/93‑1 (MQ=255)

                       cAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAtgt                                                                 >  1:114309/1‑93 (MQ=255)

                        aaCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAtgtg                                                                >  1:3991125/1‑93 (MQ=255)

                        aaCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAtgtg                                                                <  2:4435630/93‑1 (MQ=255)

                        aaCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAtgtg                                                                >  2:119455/1‑93 (MQ=255)

                        aaCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGAtgtg                                                                >  1:4294491/1‑93 (MQ=255)

                         aCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGc                                                               >  2:1443415/1‑93 (MQ=255)

                         aCCAGGGGTTGCAGAACGTATCGTTTGGTCAGGGGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGc                                                               <  1:4576508/93‑1 (MQ=255)

                           cAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGctc                                                             >  2:2096635/1‑93 (MQ=255)

                             ggggTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGctctc                                                           >  2:2297794/1‑93 (MQ=255)

                                 ttGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGaa                                                       >  2:377189/1‑93 (MQ=255)

                                    cAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAg                                                    >  1:3880402/1‑93 (MQ=255)

                                     aGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGc                                                   >  1:3567058/1‑93 (MQ=255)

                                     aGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGc                                                   >  2:1681609/1‑93 (MQ=255)

                                      gAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCa                                                  <  2:1394583/93‑1 (MQ=255)

                                      gAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCa                                                  >  2:3918306/1‑93 (MQ=255)

                                      gAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCa                                                  <  1:3405689/93‑1 (MQ=255)

                                      gAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCa                                                  >  1:3015638/1‑93 (MQ=255)

                                         cGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGa                                               <  2:4341768/93‑1 (MQ=255)

                                         cGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGa                                               <  2:829758/93‑1 (MQ=255)

                                          gTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGag                                              >  1:3117242/1‑93 (MQ=255)

                                            aTCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGagca                                            >  1:522697/1‑93 (MQ=255)

                                            aTCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGagca                                            <  2:1751295/93‑1 (MQ=255)

                                             tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGagcag                                           >  1:1862121/1‑93 (MQ=255)

                                              cGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGagcagc                                          >  2:1576574/1‑93 (MQ=255)

                                               gTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAgcagca                                         >  2:2966524/1‑93 (MQ=255)

                                                 ttGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAGAGCACGAGCAGCAta                                       >  1:3298878/1‑92 (MQ=255)

                                                 ttGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATc                                       <  2:2366796/93‑1 (MQ=255)

                                                 ttGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATc                                       >  1:3298877/1‑93 (MQ=255)

                                                  tGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCa                                      <  2:3144888/93‑1 (MQ=255)

                                                  tGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCa                                      <  1:2321189/93‑1 (MQ=255)

                                                  gggTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCa                                      <  2:188414/92‑1 (MQ=255)

                                                   ggTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAg                                     <  2:2151812/93‑1 (MQ=255)

                                                   ggTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAg                                     <  2:978880/93‑1 (MQ=255)

                                                     tCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGa                                   >  2:3223562/1‑93 (MQ=255)

                                                      cAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGAt                                  <  2:215592/93‑1 (MQ=255)

                                                      cAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGAt                                  <  2:215593/93‑1 (MQ=255)

                                                      cAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATACTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGAt                                  >  1:104920/1‑93 (MQ=255)

                                                          tGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGt                              <  1:2854515/93‑1 (MQ=255)

                                                          tGAACAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGt                              <  2:2598422/93‑1 (MQ=255)

                                                           gAACAGCGCGCGTCTCAACGTGAACGTGCGCGAAACGCACGCAATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGtt                             <  2:4003150/93‑1 (MQ=255)

                                                             aCAGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTa                           >  2:123554/1‑93 (MQ=255)

                                                               aGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAAc                         <  1:486655/93‑1 (MQ=255)

                                                               aGCGTGCGTCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAAc                         <  1:4744278/93‑1 (MQ=255)

                                                                gcgtgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACa                        >  2:719536/1‑93 (MQ=255)

                                                                gcgtgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACa                        >  2:4457512/1‑93 (MQ=255)

                                                                 cgtgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAg                       >  1:373876/1‑93 (MQ=255)

                                                                  gtgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGTATCTGTTTAACAGc                      <  1:4574670/93‑1 (MQ=255)

                                                                  gtgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGc                      <  1:71389/93‑1 (MQ=255)

                                                                   tgcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCa                     <  2:4668721/93‑1 (MQ=255)

                                                                    gcgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCaa                    <  1:2550525/93‑1 (MQ=255)

                                                                     cttctcAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAAc                   <  2:115637/91‑1 (MQ=255)

                                                                     cgtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAAc                   <  2:4016363/93‑1 (MQ=255)

                                                                      gtCTCAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGAACGAGCAGCATCAGGATATGTTTAACAGCAAcc                  >  1:2260208/1‑93 (MQ=255)

                                                                       tctcAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCg                 >  2:933281/1‑93 (MQ=255)

                                                                       tctcAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCg                 >  2:3937385/1‑93 (MQ=255)

                                                                         tcAACGTGAACGTGCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGtt               >  1:4194355/1‑93 (MQ=255)

                                                                            acgtgaacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGa            >  1:4015802/1‑93 (MQ=255)

                                                                            acgtgaacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGa            >  2:1150160/1‑93 (MQ=255)

                                                                             cgtgaacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGAc           <  1:4362481/93‑1 (MQ=255)

                                                                              gtgaacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGAcc          >  1:1027876/1‑93 (MQ=255)

                                                                               tgaacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCt         >  2:1050528/1‑93 (MQ=255)

                                                                                 aacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTtc       <  1:823866/93‑1 (MQ=255)

                                                                                 aacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTtc       <  1:2068342/93‑1 (MQ=255)

                                                                                 aacgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTtc       <  1:1518554/93‑1 (MQ=255)

                                                                                  acgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTtct      <  1:3577865/93‑1 (MQ=255)

                                                                                  acgtgCGCGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTtct      <  1:2440879/93‑1 (MQ=255)

                                                                                    gtgcgcgAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCg    <  2:1392928/93‑1 (MQ=255)

                                                                                    gtgcgcgAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCg    <  1:640118/93‑1 (MQ=255)

                                                                                    gtgcgcgAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCg    >  2:2049180/1‑93 (MQ=255)

                                                                                      gcgcgAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCGGa  <  1:3208207/93‑1 (MQ=255)

GAGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAACAGCACGAGCAGCATCAGGATCTGTTTAACAGCAACCGTTTGACCTTCTCGGA  >  NC_000913/979299‑979477

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