BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,404	+G	coding (487/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,402	1	.	G	100.0%	458.1 / NA	131	L162R (CTG→CGG)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (50/81); total (50/81)

Only 100 of 131 total aligned reads displayed.

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AA‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008346‑4008458

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:5033459/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:446582/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:8159967/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:2237438/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 >  1:337657/1‑93 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:1392991/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:4297949/93‑1 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:6922214/93‑2 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  2:2845594/93‑2 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                >  1:3614701/1‑92 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:5418012/93‑2 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:6065576/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:5872788/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               >  2:4163336/1‑91 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:4769163/1‑91 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:6064375/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:750050/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:2317198/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:8637297/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:414086/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:5805307/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTt‑‑ttgg                           >  2:6658038/1‑92 (MQ=255)

    gAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:8637298/93‑5 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  1:8162548/1‑93 (MQ=255)

           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:5721445/93‑1 (MQ=255)

           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:4840526/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:4509590/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:8788330/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     >  2:4959575/1‑93 (MQ=255)

            ggcgggcgCCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:421304/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:1244670/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:2455348/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:2829196/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:2455349/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:7587141/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAAGCTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  1:8011368/93‑1 (MQ=255)

             gcgggcggCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:1244672/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   >  1:1004117/1‑93 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:8090470/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:2568450/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  1:4501072/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:1912031/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:497721/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:3664216/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:8639006/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:2848457/93‑1 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:5213019/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:5543527/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  1:6847877/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  2:3824349/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  2:3764502/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  1:6357237/93‑1 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                >  1:6339990/1‑93 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               >  2:3745720/1‑93 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:852079/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:6245668/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:1871314/93‑1 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:1556823/1‑93 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:6899187/1‑93 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:4371439/1‑93 (MQ=255)

                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:1235066/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:3429348/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  1:8597145/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  1:2683414/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:511017/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  1:403474/1‑93 (MQ=255)

                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa         >  2:7031658/5‑93 (MQ=255)

                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:6546955/4‑93 (MQ=255)

                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:144119/4‑93 (MQ=255)

                         tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  1:4387540/4‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  2:3401072/3‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  1:8442106/3‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  1:3273534/3‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:5036154/91‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  2:412049/91‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1868207/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:3764672/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:5353936/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  1:2240486/92‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:5855352/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:5855353/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:8046208/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1302241/2‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCGGAAATg     >  1:2685785/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:8732540/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     <  1:1458219/93‑1 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:6726358/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  1:7691908/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  2:3973376/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2171973/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  2:4465606/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:5856790/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:3801018/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:26644/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:7888353/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:7935961/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:2102219/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:4728874/93‑1 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:7788741/93‑1 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAA‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:2683727/93‑1 (MQ=255)

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AA‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008346‑4008458

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