BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,091,676	A→G	Q53R (CAG→CGG)	yqgF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,091,676	0	A	G	97.6%	391.4 / ‑0.2	126	Q53R (CAG→CGG)	yqgF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (55/68); total (56/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

Only 100 of 137 total aligned reads displayed.

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                          <  1:3343993/93‑1 (MQ=255)

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                          >  2:46368/1‑93 (MQ=255)

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                          <  1:2370106/93‑1 (MQ=255)

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                          >  2:711809/1‑93 (MQ=255)

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                          >  2:454977/1‑93 (MQ=255)

    cGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGa                                                                                      <  1:3372092/93‑1 (MQ=255)

     ggCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGAc                                                                                     <  2:197554/93‑1 (MQ=255)

      gCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACg                                                                                    >  1:2519335/1‑93 (MQ=255)

          cGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAat                                                                                <  2:1130618/93‑1 (MQ=255)

          cGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAat                                                                                <  2:3390979/93‑1 (MQ=255)

           gCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatc                                                                               <  1:4227797/93‑1 (MQ=255)

            cTCGCGCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCTGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatca                                                                              <  1:888492/93‑1 (MQ=255)

            cTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatca                                                                              <  2:1659325/93‑1 (MQ=255)

            cTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatca                                                                              <  2:4241701/93‑1 (MQ=255)

             tCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGACGAAaacat                                                                             >  1:3352853/1‑93 (MQ=255)

             tCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGCAattat                                                                             >  1:3127591/1‑93 (MQ=255)

              cGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatcatc                                                                            <  1:3179576/93‑1 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGGCGGAaacattg                                                                           >  2:515198/1‑91 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGACGCAATCAtcg                                                                           >  1:4308143/1‑93 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcg                                                                           <  1:4575673/93‑1 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTAAGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCCGACGAAATAAtcg                                                                           >  1:3920823/1‑93 (MQ=255)

                gccTTTGTCTGCAATGAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgt                                                                          <  2:4725256/92‑1 (MQ=255)

                cccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgt                                                                          <  1:187964/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCAtcgtc                                                                         >  1:2162164/1‑93 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtc                                                                         <  2:386368/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtc                                                                         <  2:1547070/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtc                                                                         <  1:2748919/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtc                                                                         <  2:347190/93‑1 (MQ=255)

                  cTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCTGACGAAAACAtcgtcg                                                                        >  2:3601226/1‑93 (MQ=255)

                   tttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCgg                                                                       >  2:3066340/1‑93 (MQ=255)

                    ttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGt                                                                      <  1:4163953/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCGGGGCGAAATCCTCGTCgggt                                                                     >  1:3618883/1‑91 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGtt                                                                     <  1:703062/93‑1 (MQ=255)

                       ccTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATAGTCGGTTTg                                                                   >  1:2572024/1‑93 (MQ=255)

                        cTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATACTCGTCGGGTTGc                                                                  >  1:2878898/1‑93 (MQ=255)

                        cTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAAACATCGTCGGttttc                                                                  >  2:1960978/1‑91 (MQ=255)

                          gCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCg                                                                <  2:4091538/93‑1 (MQ=255)

                          gCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCg                                                                <  2:4091539/93‑1 (MQ=255)

                            aaTTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGGCGGAATCCTACTCGGTTTGCCGCt                                                              >  2:1868357/1‑93 (MQ=255)

                            aaTCAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGACGGACGAAATCATCGTCGGTTTGCCGCt                                                              <  2:493085/93‑1 (MQ=255)

                              ttAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGAAGAAATTATAGTCGGGTTGCCCCTGa                                                            >  1:1248430/1‑93 (MQ=255)

                              ttAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGa                                                            <  1:90751/93‑1 (MQ=255)

                               tAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGaa                                                           <  1:1306635/93‑1 (MQ=255)

                                aaaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAAACATCGTCGGTTTGCCGCTGAat                                                          >  1:131595/1‑93 (MQ=255)

                                aaaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAat                                                          <  1:3852775/93‑1 (MQ=255)

                                 aaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAata                                                         <  2:2200291/93‑1 (MQ=255)

                                 aaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAata                                                         <  1:3607300/93‑1 (MQ=255)

                                  aGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAata‑t                                                       <  1:3176801/93‑1 (MQ=255)

                                   gCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑Tg                                                      <  2:2301963/93‑1 (MQ=255)

                                   acacaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑Tg                                                      <  2:3853536/92‑1 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGa                                                    <  2:114815/93‑1 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGa                                                    <  2:1773997/93‑1 (MQ=255)

                                        ggACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGGTTGCCGCTGAATA‑Tgggagg                                                 >  2:890309/1‑89 (MQ=255)

                                           cGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATA‑TGGACGGGAc                                              >  2:2502797/1‑93 (MQ=255)

                                            ggTACGCCGGACTGGAACATTATCGAGCGTTTACTTAAAGAGTGGCGGCCGGACGAAAACATCGTCGGGTTGCCGCTGAATA‑TGGGCGGCAcc                                             >  2:541214/1‑93 (MQ=255)

                                            ggTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCAcc                                             <  2:4499740/93‑1 (MQ=255)

                                              tACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTTGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGa                                           <  1:4622898/93‑1 (MQ=255)

                                               aCGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGGAGGCAACgcg                                          >  1:3041552/1‑91 (MQ=255)

                                                cGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTTGATA‑TTGAGGGCACCGAGc                                         >  1:2589635/1‑93 (MQ=255)

                                                 gCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCa                                        <  2:2977892/93‑1 (MQ=255)

                                                    ggACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGACGAAATCATCGTCGGTTTGTCGCTGAATA‑TGGACGGCACCGAGCAAcc                                     >  2:337889/1‑93 (MQ=255)

                                                    ggACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAAcc                                     <  1:3347440/93‑1 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGGTTGCCGCTGGATA‑TTGACGGGACCCAGGAAccc                                    >  2:4188193/1‑92 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATA‑TGGACGGCACCAGACAACCa                                    >  2:3651193/1‑93 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCATCGGTTTGTCGCTGAATA‑TAGGCGGGCGCGAGCaaacc                                    >  1:5130/1‑90 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTGACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCa                                    <  2:1547769/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCAt                                   <  1:2985587/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCCGACGAAAACATCGTCGGGTTGTCGCTGTATA‑TGGGCGGCACCGGCGAACCAt                                   >  2:1715558/1‑93 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTTGGGTTGCCGCTGAAAATTGGAGGGC‑CCGAGCAACCCtt                                  >  2:3870655/1‑93 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCAtt                                  <  1:3219289/93‑1 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCAtt                                  <  2:3993331/93‑1 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGGCGGCACCCAGCAACCATTg                                 >  2:1549919/1‑93 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATA‑TGGGCGGCACCGAGCAACCATTg                                 >  1:4659334/1‑93 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGACGAAATCGTATGCGGTTTTCCTCTTAATA‑TGGACGGCACAGTGCAACCagttc                                >  2:2871603/1‑89 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGTCGCTGAATA‑TGGACGGCACCGCGCCACCATTGa                                >  1:1927436/1‑93 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCGTCGGTTTGCCGCTGGATA‑TGGACGGCAACCAGCAAACATTGa                                >  2:3134546/1‑93 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTTAc                               >  2:1483275/1‑93 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGAc                               <  2:2341582/93‑1 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGAc                               <  2:4628212/93‑1 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGAc                               <  2:1222939/93‑1 (MQ=255)

                                                           aaCATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACt                              <  1:4440439/93‑1 (MQ=255)

                                                             cATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCACTGACTGc                            >  1:1015805/1‑93 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCCACCACTGACTGCCa                          >  1:474763/1‑93 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCa                          >  1:2221997/1‑93 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGCGGCCCGACGAAATCATCGTCGGTTTTCCGCTTAATA‑TTGACGGCACCGAGCAACCATTTACTTCCa                          >  2:3594831/1‑93 (MQ=255)

                                                                 aTCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCaga                        <  1:1375386/93‑1 (MQ=255)

                                                                 aTCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCaga                        <  2:1292593/93‑1 (MQ=255)

                                                                  tCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCGTCGGTTTGGCGCTGAATA‑TGGACGGCACCCAGCAACCATTGGCTGCCacag                       >  2:3841212/1‑93 (MQ=255)

                                                                   cGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAgc                      <  2:4440032/93‑1 (MQ=255)

                                                                    gAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAgcg                     <  1:2303879/93‑1 (MQ=255)

                                                                     aGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAgcgc                    >  2:2368047/1‑93 (MQ=255)

                                                                          ttACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTaaa               <  2:2064568/93‑1 (MQ=255)

                                                                                aaaGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGcc         <  1:1820779/93‑1 (MQ=255)

                                                                                   gagTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAAc      <  1:2431492/93‑1 (MQ=255)

                                                                                   gagTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGAGTAAATTTGCCAAc      <  1:4150904/93‑1 (MQ=255)

                                                                                    agTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAAcc     <  2:1107044/93‑1 (MQ=255)

                                                                                     gTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCg    >  2:3870698/1‑93 (MQ=255)

                                                                                     gTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCg    <  2:1779008/93‑1 (MQ=255)

                                                                                       ggcggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  <  2:1651475/93‑1 (MQ=255)

                                                                                       ggcggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  <  1:1000622/93‑1 (MQ=255)

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATA‑TGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

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