BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	361,958	A→G	V149A (GTT→GCT)	lacY ←	galactoside permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	361,958	0	A	G	80.7%	414.1 / 102.9	207	V149A (GTT→GCT)	lacY	galactoside permease
Reads supporting (aligned to +/- strand): ref base A (21/19); new base G (65/102); total (86/121)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 219 total aligned reads displayed.

GATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTA  >  NC_000913/361867‑362050

gATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAac                                                                                             <  2:7234632/93‑1 (MQ=255)

 aTGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAgca                                                                                            <  1:679287/93‑1 (MQ=255)

  tGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAgcac                                                                                           >  1:4955673/1‑93 (MQ=255)

  tGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAgcac                                                                                           >  1:3878414/1‑93 (MQ=255)

   gAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAAcaca                                                                                          <  2:6933561/93‑1 (MQ=255)

     gTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGc                                                                                        <  2:2771342/93‑1 (MQ=255)

      tGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGcc                                                                                       <  2:1627574/93‑1 (MQ=255)

       gCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCa                                                                                      >  1:920572/1‑93 (MQ=255)

         acaGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCaaa                                                                                    <  1:2959135/93‑1 (MQ=255)

         acaGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCaaa                                                                                    >  2:2009219/1‑93 (MQ=255)

            gCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACAt                                                                                 <  2:6243487/93‑1 (MQ=255)

             ccAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATc                                                                                <  2:6183903/93‑1 (MQ=255)

             ccAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATc                                                                                <  1:2639049/93‑1 (MQ=255)

              cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATcc                                                                               <  2:5599421/93‑1 (MQ=255)

              cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATcc                                                                               <  2:3479226/93‑1 (MQ=255)

              cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATcc                                                                               <  2:3168237/93‑1 (MQ=255)

               agagCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcg                                                                              >  1:2872946/1‑93 (MQ=255)

                gagCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgc                                                                             >  1:4600786/1‑93 (MQ=255)

                 agCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcg                                                                            <  1:651548/93‑1 (MQ=255)

                  gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgc                                                                           <  2:1242907/93‑1 (MQ=255)

                  gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgc                                                                           <  1:7209960/93‑1 (MQ=255)

                  gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgc                                                                           >  1:3100478/1‑93 (MQ=255)

                  gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCcgcgc                                                                           <  1:511021/93‑1 (MQ=255)

                     cagccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgcgcg                                                                        >  2:26946/1‑93 (MQ=255)

                     cagccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgcgcg                                                                        >  1:6675094/1‑93 (MQ=255)

                     cagccagAAAACAAACTGATTATTGAGGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCcgcgcgcg                                                                        <  1:3344528/93‑1 (MQ=255)

                        ccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGAcc                                                                     >  2:6412413/1‑93 (MQ=255)

                        ccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGAcc                                                                     >  1:178996/1‑93 (MQ=255)

                         cagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGTCCAGCCAGCACAGCCAAACATCCGCGCGCGACCa                                                                    <  2:514128/93‑1 (MQ=255)

                          agAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCaa                                                                   >  2:735141/1‑93 (MQ=255)

                           gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCaaa                                                                  >  2:2407004/1‑93 (MQ=255)

                           gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCaaa                                                                  >  1:7400298/1‑93 (MQ=255)

                           gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCaaa                                                                  <  1:5867253/93‑1 (MQ=255)

                           gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCaaa                                                                  >  1:4252614/1‑93 (MQ=255)

                           gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCaaa                                                                  <  1:7092700/93‑1 (MQ=255)

                              aacaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTc                                                               >  1:4642867/1‑93 (MQ=255)

                               acaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCg                                                              >  1:2369327/1‑93 (MQ=255)

                               acaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCg                                                              >  2:185299/1‑93 (MQ=255)

                               acaaacTGATTATTGATGG‑GAACCTGATGCCGACAATCGAGGCACACAGCGCCCAGGCAGCACAGCCAAACATCCGCGAGCGGCCAAATTCGa                                                             >  1:631115/1‑93 (MQ=255)

                                 aaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGaa                                                            >  2:39590/1‑93 (MQ=255)

                                 aaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGaa                                                            >  2:3317121/1‑93 (MQ=255)

                                  aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                                           >  2:2947935/1‑93 (MQ=255)

                                  aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                                           >  2:6505238/1‑93 (MQ=255)

                                    cTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAAtt                                                         >  1:2507332/1‑93 (MQ=255)

                                     tGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTa                                                        <  1:1939540/93‑1 (MQ=255)

                                       attattGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACt                                                      <  2:439697/93‑1 (MQ=255)

                                        ttattGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTg                                                     <  2:3761496/93‑1 (MQ=255)

                                        ttattGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTg                                                     >  1:577661/1‑93 (MQ=255)

                                         tattGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGc                                                    <  1:5699540/93‑1 (MQ=255)

                                           ttGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGa                                                  <  1:1798730/93‑1 (MQ=255)

                                           ttGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGa                                                  <  2:7265644/93‑1 (MQ=255)

                                           ttGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGa                                                  >  1:2784560/1‑93 (MQ=255)

                                            tGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGAc                                                 <  2:2818450/93‑1 (MQ=255)

                                             gATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACg                                                <  2:5583570/93‑1 (MQ=255)

                                             gATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACg                                                <  1:3637772/93‑1 (MQ=255)

                                             gATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACg                                                >  2:5352398/1‑93 (MQ=255)

                                              aTGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACgg                                               >  2:50302/1‑93 (MQ=255)

                                              aTGGGGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACgg                                               <  1:3841667/93‑1 (MQ=255)

                                               tGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGc                                              >  1:5616991/1‑93 (MQ=255)

                                               tGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGc                                              <  1:1265404/93‑1 (MQ=255)

                                               tGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGc                                              >  1:2027629/1‑93 (MQ=255)

                                                  tGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGa                                           >  1:4343489/1‑93 (MQ=255)

                                                    aaCATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACt                                         <  1:6283899/93‑1 (MQ=255)

                                                      cATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACttt                                       <  2:2082269/93‑1 (MQ=255)

                                                      cATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACttt                                       <  2:4371010/93‑1 (MQ=255)

                                                       atgatgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCTAAATTACTGCGACGGCTGACTTtc                                      >  2:4789821/1‑93 (MQ=255)

                                                        tgatgCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtct                                     >  2:5412810/1‑93 (MQ=255)

                                                        tgatgCCGACAAACGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtct                                     >  1:1063411/1‑93 (MQ=255)

                                                         gatgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtctc                                    >  2:7158028/1‑93 (MQ=255)

                                                           tgCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCaa                                  <  1:219772/93‑1 (MQ=255)

                                                            gCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                 <  2:5924651/93‑1 (MQ=255)

                                                            gCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                 >  1:3366441/1‑93 (MQ=255)

                                                            gCCGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                 <  1:862222/93‑1 (MQ=255)

                                                             ccGACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATa                                >  1:5572325/1‑93 (MQ=255)

                                                               gACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa                              <  1:4389356/93‑1 (MQ=255)

                                                               gACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa                              <  1:1470236/93‑1 (MQ=255)

                                                               gACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa                              <  1:4527906/93‑1 (MQ=255)

                                                               gACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa                              <  2:3353493/93‑1 (MQ=255)

                                                               gACAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa                              <  2:5529012/93‑1 (MQ=255)

                                                                aCAATCGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAAt                             >  1:5994660/1‑93 (MQ=255)

                                                                aCAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAAt                             >  1:4278890/1‑93 (MQ=255)

                                                                     cGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc                        <  2:2522329/93‑1 (MQ=255)

                                                                     cGAGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc                        <  1:5471529/93‑1 (MQ=255)

                                                                       aGGCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTa                      >  2:6755695/1‑93 (MQ=255)

                                                                       aGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTa                      >  1:1155495/1‑93 (MQ=255)

                                                                        ggCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc                     <  2:4176482/93‑1 (MQ=255)

                                                                        ggCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc                     <  1:5977152/93‑1 (MQ=255)

                                                                         gCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAct                    <  1:6070515/93‑1 (MQ=255)

                                                                         gCACACAGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAct                    >  2:2634115/1‑93 (MQ=255)

                                                                            cacaGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActgct                 <  2:3867810/93‑1 (MQ=255)

                                                                             acaGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActgctg                >  2:1092208/1‑93 (MQ=255)

                                                                              caGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTgg               <  2:1979392/93‑1 (MQ=255)

                                                                              caGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTgg               >  2:7309570/1‑93 (MQ=255)

                                                                              caGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTgg               <  2:1376514/93‑1 (MQ=255)

                                                                               aGCGCCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGgc              >  1:6948086/1‑93 (MQ=255)

                                                                                  gcCCAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCa           >  1:3508802/1‑93 (MQ=255)

                                                                                   cccAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCAc          >  1:652601/1‑93 (MQ=255)

                                                                                   cccAGCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCAc          <  2:1508692/93‑1 (MQ=255)

                                                                                    ccagccagCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCAcc         <  1:4997031/93‑1 (MQ=255)

                                                                                       gCCAGCACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGc      <  2:6701907/93‑1 (MQ=255)

GATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTA  >  NC_000913/361867‑362050

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