BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,227,714	A→G	L3P (CTG→CCG)	pepE ←	peptidase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,227,714	0	A	G	88.0%	729.1 / 82.2	266	L3P (CTG→CCG)	pepE	peptidase E
Reads supporting (aligned to +/- strand): ref base A (16/16); new base G (113/121); total (129/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 267 total aligned reads displayed.

GAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCG  >  NC_000913/4227624‑4227803

gAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGt                                                                                         >  2:5729540/1‑93 (MQ=255)

 aGCGGCGACCCTGTAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGtt                                                                                        <  2:5183044/93‑1 (MQ=255)

  gcggcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTc                                                                                       >  1:4880/1‑93 (MQ=255)

  gcggcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTc                                                                                       <  1:2431455/93‑1 (MQ=255)

    ggcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCa                                                                                     <  2:1672821/93‑1 (MQ=255)

     gcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCAt                                                                                    <  2:510282/93‑1 (MQ=255)

     gcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCAt                                                                                    <  2:5637744/93‑1 (MQ=255)

     gcgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGCCTTGACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCAt                                                                                    <  1:594680/93‑1 (MQ=255)

      cgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCAtt                                                                                   <  1:271027/93‑1 (MQ=255)

      cgACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCAGTTCCAtt                                                                                   >  1:843087/1‑93 (MQ=255)

        aCCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTAc                                                                                 <  2:3664022/93‑1 (MQ=255)

        aCCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTAc                                                                                 <  1:1328333/93‑1 (MQ=255)

        aCCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTAc                                                                                 >  1:315034/1‑93 (MQ=255)

         cccTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACt                                                                                >  2:5354648/1‑93 (MQ=255)

          ccTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACtt                                                                               >  2:719452/1‑93 (MQ=255)

           cTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTtc                                                                              <  1:70221/93‑1 (MQ=255)

           cTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTtc                                                                              <  1:913436/93‑1 (MQ=255)

           cTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCAGTTCCATTACTtc                                                                              >  2:3013302/1‑93 (MQ=255)

            tGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTtct                                                                             >  1:5867/1‑93 (MQ=255)

             gCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTtctc                                                                            >  1:1117202/1‑93 (MQ=255)

             gCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTtctc                                                                            >  2:5904977/1‑93 (MQ=255)

              cAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTcc                                                                           >  1:5690297/1‑93 (MQ=255)

               aaCTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCg                                                                          <  1:4615324/93‑1 (MQ=255)

               aaCTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCg                                                                          >  2:1892761/1‑93 (MQ=255)

                  tGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGtt                                                                       <  2:3155203/93‑1 (MQ=255)

                  tGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGtt                                                                       >  2:1958835/1‑93 (MQ=255)

                      cAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTa                                                                   >  2:2024575/1‑93 (MQ=255)

                           aTTAGCGGCAGTGCATGTTCCAGCCAGTCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGaaa                                                              <  2:554834/93‑1 (MQ=255)

                           aTTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGaaa                                                              <  2:4029419/93‑1 (MQ=255)

                            ttAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAAt                                                             >  1:5567965/1‑93 (MQ=255)

                            ttAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAAt                                                             <  2:2942442/93‑1 (MQ=255)

                             tAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATc                                                            <  1:4523455/93‑1 (MQ=255)

                              aGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCa                                                           <  2:65603/93‑1 (MQ=255)

                              aGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCa                                                           >  2:1148623/1‑93 (MQ=255)

                               gCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCAc                                                          <  2:585838/93‑1 (MQ=255)

                               gCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCAc                                                          <  1:2233061/93‑1 (MQ=255)

                               gCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCAc                                                          <  1:855480/93‑1 (MQ=255)

                                  gCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACggg                                                       >  2:3840413/1‑93 (MQ=255)

                                   cAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGa                                                      <  1:3592978/93‑1 (MQ=255)

                                    agggCATGTTCCAGCCAGGCTTTACCGGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAg                                                     <  2:1022343/90‑1 (MQ=255)

                                    aGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAg                                                     <  2:4570457/93‑1 (MQ=255)

                                     gTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAgg                                                    >  1:6435856/1‑93 (MQ=255)

                                     gTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAgg                                                    >  2:4729674/1‑93 (MQ=255)

                                      tGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                   >  1:4844670/1‑93 (MQ=255)

                                       gCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                  <  1:6292262/93‑1 (MQ=255)

                                        cATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCAt                                                 <  2:6467212/93‑1 (MQ=255)

                                        cATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCAt                                                 <  1:4775360/93‑1 (MQ=255)

                                         aTGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATg                                                >  2:3318098/1‑93 (MQ=255)

                                         aTGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATg                                                >  2:4455633/1‑93 (MQ=255)

                                           gTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCt                                              >  2:3716644/1‑93 (MQ=255)

                                             tCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTAc                                            >  1:5246095/1‑93 (MQ=255)

                                             tCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTAc                                            >  1:2557547/1‑93 (MQ=255)

                                               cagccagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCa                                          >  1:6183132/1‑93 (MQ=255)

                                                 gccagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACt                                        >  1:784250/1‑93 (MQ=255)

                                                 gccagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACt                                        >  2:6032707/1‑93 (MQ=255)

                                                 gccagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACt                                        >  1:1048071/1‑93 (MQ=255)

                                                  ccagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTg                                       >  2:1312678/1‑93 (MQ=255)

                                                  ccagGCTTTACCCGGCAGCGTCGAGTGACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTg                                       <  1:3383519/93‑1 (MQ=255)

                                                   cagGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGc                                      <  1:616833/93‑1 (MQ=255)

                                                     ggCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGa                                    >  2:2253131/1‑93 (MQ=255)

                                                     ggCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGGTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGa                                    >  1:3970642/1‑93 (MQ=255)

                                                      gCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGaa                                   <  2:4951556/93‑1 (MQ=255)

                                                        tttACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAAtt                                 >  2:4079027/1‑93 (MQ=255)

                                                         ttACCCGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAAttt                                <  1:1925641/93‑1 (MQ=255)

                                                            cccGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAc                             <  1:953482/93‑1 (MQ=255)

                                                            cccGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAc                             <  2:378696/93‑1 (MQ=255)

                                                             ccGGCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAcc                            >  1:5145543/1‑93 (MQ=255)

                                                                gCAGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                         <  1:1727666/93‑1 (MQ=255)

                                                                  aGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTa                       >  2:5285706/1‑93 (MQ=255)

                                                                  aGCGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTa                       >  2:4219002/1‑93 (MQ=255)

                                                                    cGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAGAACTTACt                     <  1:6539018/93‑1 (MQ=255)

                                                                    cGTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACt                     <  1:4689228/93‑1 (MQ=255)

                                                                    cGTCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACt                     <  2:2248091/93‑1 (MQ=255)

                                                                     gTCGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTg                    >  2:786485/1‑93 (MQ=255)

                                                                     gTCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTg                    >  2:6304893/1‑93 (MQ=255)

                                                                      tCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGa                   >  1:2202252/1‑93 (MQ=255)

                                                                       cGAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGAc                  >  2:2080870/1‑93 (MQ=255)

                                                                        gAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACa                 >  2:5940677/1‑93 (MQ=255)

                                                                        gAGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACa                 <  1:4236211/93‑1 (MQ=255)

                                                                         aGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACAt                >  1:1652869/1‑93 (MQ=255)

                                                                         aGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACAt                <  1:296253/93‑1 (MQ=255)

                                                                         aGTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACAt                >  2:3489934/1‑93 (MQ=255)

                                                                         aGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACAt                <  1:280877/93‑1 (MQ=255)

                                                                          gTTACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATg               >  2:2513085/1‑93 (MQ=255)

                                                                            tACTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGc             >  2:653665/1‑93 (MQ=255)

                                                                             aCTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGcc            >  1:2521132/1‑92 (MQ=255)

                                                                              cTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTc           <  1:1799411/93‑1 (MQ=255)

                                                                              cTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTc           <  1:729025/93‑1 (MQ=255)

                                                                              cTCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTc           >  2:3946447/1‑93 (MQ=255)

                                                                               tCAATAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCa          >  2:1636628/1‑93 (MQ=255)

                                                                               tCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCa          >  2:2789405/1‑93 (MQ=255)

                                                                                  aTAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGa       >  2:2663420/1‑93 (MQ=255)

                                                                                  aTAAAAGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGa       >  1:5322520/1‑93 (MQ=255)

                                                                                     aaaGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGt    >  2:3260861/1‑93 (MQ=255)

                                                                                     aaaGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGt    >  1:1481665/1‑93 (MQ=255)

                                                                                     aaaGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGt    <  2:2353352/93‑1 (MQ=255)

                                                                                      aaGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTc   >  1:2352625/1‑93 (MQ=255)

                                                                                      aaGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTc   >  2:3052858/1‑93 (MQ=255)

                                                                                       aGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg  <  2:6362230/93‑1 (MQ=255)

                                                                                       aGCGGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg  <  1:1180024/93‑1 (MQ=255)

GAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGGCTTTACCCGGCAGCGTCGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCG  >  NC_000913/4227624‑4227803

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: `‑`