BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,587	T→C	S116G (AGT→GGT)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,587	0	T	C	98.3%	349.6 / ‑2.8	117	S116G (AGT→GGT)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (63/52); total (64/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 118 total aligned reads displayed.

ACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAC  >  NC_000913/1616499‑1616678

aCTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGt                                                                                         >  2:1519030/1‑93 (MQ=255)

aCTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGt                                                                                         >  2:647431/1‑93 (MQ=255)

 cTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGtt                                                                                        >  1:3721943/1‑93 (MQ=255)

 cTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGtt                                                                                        >  1:35425/1‑93 (MQ=255)

   gtgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCa                                                                                      <  1:2054243/93‑1 (MQ=255)

     ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                    <  1:4634019/93‑1 (MQ=255)

      gACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCt                                                                                   <  2:2601674/93‑1 (MQ=255)

      gACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCt                                                                                   <  1:2498336/93‑1 (MQ=255)

        cGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTc                                                                                 <  2:3363654/93‑1 (MQ=255)

          ccGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCa                                                                               <  1:567541/93‑1 (MQ=255)

           cGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAg                                                                              >  2:942443/1‑93 (MQ=255)

           cGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAg                                                                              >  1:199521/1‑93 (MQ=255)

            gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                             <  1:111831/93‑1 (MQ=255)

            gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                             >  2:1562423/1‑93 (MQ=255)

             aGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGct                                                                            >  1:3151323/1‑93 (MQ=255)

              gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                           >  1:2488316/1‑93 (MQ=255)

              gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                           >  1:3184700/1‑93 (MQ=255)

               cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctct                                                                          <  1:3106884/93‑1 (MQ=255)

                taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                         >  2:975761/1‑93 (MQ=255)

                 aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                        <  2:1962416/93‑1 (MQ=255)

                 aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                        <  2:3921477/93‑1 (MQ=255)

                  ataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTcc                                                                       >  2:789450/1‑93 (MQ=255)

                   taatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCg                                                                      >  2:1924793/1‑93 (MQ=255)

                   taatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCg                                                                      >  1:3580103/1‑93 (MQ=255)

                   taatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCg                                                                      <  2:3155954/93‑1 (MQ=255)

                     atcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGAc                                                                    <  1:2175367/93‑1 (MQ=255)

                     atcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGAc                                                                    <  1:3901687/93‑1 (MQ=255)

                      tcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGGTCATCTTCCAGCTCTTCCGACt                                                                   >  2:3819833/1‑93 (MQ=255)

                           gCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                              >  1:71405/1‑93 (MQ=255)

                           gCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                              <  2:1359575/93‑1 (MQ=255)

                           gCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                              >  2:4507564/1‑93 (MQ=255)

                            cTATTGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                             <  1:3011819/93‑1 (MQ=255)

                            cTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                             >  1:1753434/1‑93 (MQ=255)

                            cTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                             <  2:1967831/93‑1 (MQ=255)

                            cTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                             >  2:1937958/1‑93 (MQ=255)

                            cTATGGTGCCCGGACCGGCAGTGCTTGGCAGTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                             <  2:4332802/93‑1 (MQ=255)

                             tATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtttt                                                            >  1:2624171/1‑93 (MQ=255)

                             tATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGGTCATCTTCCAGCTCTTCCGACTTGCtttt                                                            >  1:360489/1‑93 (MQ=255)

                               tGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGCTTCATCTTCCAGCTCTTCCGACTTGCTTTTCg                                                          <  1:4340053/93‑1 (MQ=255)

                                ggTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGc                                                         >  1:4531657/1‑93 (MQ=255)

                                  ttCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCt                                                       <  2:3316784/93‑1 (MQ=255)

                                   tCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTc                                                      >  2:3276302/1‑93 (MQ=255)

                                   tCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTc                                                      <  1:3264672/93‑1 (MQ=255)

                                    cccGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTTGCCTcc                                                     >  1:1725899/1‑93 (MQ=255)

                                      cGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCgg                                                   >  1:2951711/1‑93 (MQ=255)

                                       ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                  >  1:1490869/1‑93 (MQ=255)

                                       ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                  >  2:1237714/1‑93 (MQ=255)

                                       ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                  >  2:2140897/1‑93 (MQ=255)

                                       ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                  >  2:2979999/1‑93 (MQ=255)

                                       ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                                  <  2:4494985/93‑1 (MQ=255)

                                        gACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTg                                                 <  1:124118/93‑1 (MQ=255)

                                        gACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTg                                                 >  2:2494451/1‑93 (MQ=255)

                                        gACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTg                                                 >  1:3594964/1‑93 (MQ=255)

                                         aCCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGa                                                >  1:1568811/1‑93 (MQ=255)

                                           cGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaat                                              >  1:4156737/1‑93 (MQ=255)

                                             gCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatca                                            >  2:1252195/1‑93 (MQ=255)

                                             gCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatca                                            >  2:1319938/1‑93 (MQ=255)

                                                gtgcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaatc                                         >  1:519053/4‑93 (MQ=255)

                                                gtgcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaatc                                         >  2:2272028/4‑93 (MQ=255)

                                                gtgcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaatc                                         <  2:2237741/90‑1 (MQ=255)

                                                  gcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGc                                       >  1:580444/2‑93 (MQ=255)

                                                   cTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCt                                      >  2:466734/1‑93 (MQ=255)

                                                    ttGGCATTGCGAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                     <  2:302346/93‑1 (MQ=255)

                                                    ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                     >  1:309783/1‑93 (MQ=255)

                                                     tGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCttt                                    >  1:2160880/1‑93 (MQ=255)

                                                      ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                                   <  2:417500/93‑1 (MQ=255)

                                                      ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAAGCGCTTTc                                   >  2:203984/1‑93 (MQ=255)

                                                       gCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTTGCCTCCGGTGAATTAATCGCTTTCt                                  >  2:1726762/1‑93 (MQ=255)

                                                           tGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTg                              <  1:3822073/93‑1 (MQ=255)

                                                                 cGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaaa                        <  1:4541217/93‑1 (MQ=255)

                                                                 cGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaaa                        <  2:4290000/93‑1 (MQ=255)

                                                                  ggCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAag                       >  2:91160/1‑93 (MQ=255)

                                                                   gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                      <  2:252213/93‑1 (MQ=255)

                                                                   gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                      <  2:1151885/93‑1 (MQ=255)

                                                                     acaAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGc                    >  2:3738950/1‑93 (MQ=255)

                                                                      caATGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa                   <  1:2605471/93‑1 (MQ=255)

                                                                      caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa                   <  2:1268033/93‑1 (MQ=255)

                                                                       aaaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCAt                  >  2:4520899/1‑93 (MQ=255)

                                                                       aaaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCAt                  <  1:4063051/93‑1 (MQ=255)

                                                                         aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc                >  2:4352163/1‑93 (MQ=255)

                                                                         aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc                >  2:3066118/1‑93 (MQ=255)

                                                                         aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc                <  1:836873/93‑1 (MQ=255)

                                                                         aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc                <  1:3975370/93‑1 (MQ=255)

                                                                         aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCACCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc                <  1:3385081/93‑1 (MQ=255)

                                                                          ggCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg               >  1:4638062/1‑93 (MQ=255)

                                                                           gCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGa              >  1:4532083/1‑93 (MQ=255)

                                                                             aaTATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaa            <  2:3360920/93‑1 (MQ=255)

                                                                             aaTATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaa            >  2:1693550/1‑93 (MQ=255)

                                                                              atatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaaa           <  1:1860656/93‑1 (MQ=255)

                                                                              atatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaaa           >  1:3628508/1‑93 (MQ=255)

                                                                              atatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaaa           <  1:495729/93‑1 (MQ=255)

                                                                               tatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAAc          >  2:3496381/1‑93 (MQ=255)

                                                                                 ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg        >  1:1960121/1‑93 (MQ=255)

                                                                                 ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg        >  2:3657984/1‑93 (MQ=255)

                                                                                  tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa       <  2:2247644/93‑1 (MQ=255)

                                                                                  tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa       <  1:3499182/93‑1 (MQ=255)

                                                                                  tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa       >  1:98764/1‑93 (MQ=255)

                                                                                   ggCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAt      >  1:2401160/1‑93 (MQ=255)

                                                                                    gCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCTGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATa     <  1:600712/93‑1 (MQ=255)

                                                                                       cTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAc  <  1:3441286/93‑1 (MQ=255)

ACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAAC  >  NC_000913/1616499‑1616678

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