BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,890,225	C→T	R345R (CGG→CGA)	b1808 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,890,225	0	C	T	99.1%	403.2 / ‑4.9	116	R345R (CGG→CGA)	b1808	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (57/58); minor base A (0/1); total (57/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 120 total aligned reads displayed.

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAG  >  NC_000913/1890135‑1890317

gAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCg                                                                                            >  2:233427/1‑93 (MQ=255)

 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                           <  1:1665639/93‑1 (MQ=255)

   cggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTAGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtgt                                                                                         <  1:4084205/93‑1 (MQ=255)

    ggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGtgtg                                                                                        <  2:852377/93‑1 (MQ=255)

       gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                     >  2:947032/1‑93 (MQ=255)

        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                    >  1:2173876/1‑93 (MQ=255)

         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                   <  2:1649237/93‑1 (MQ=255)

          agagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAc                                                                                  <  1:2146308/93‑1 (MQ=255)

           gagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACa                                                                                 <  2:4327275/93‑1 (MQ=255)

            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                                <  1:2081743/93‑1 (MQ=255)

            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                                >  2:1236276/1‑93 (MQ=255)

             gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                               <  1:516693/93‑1 (MQ=255)

               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                             <  1:21508/93‑1 (MQ=255)

               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                             >  1:4029300/1‑93 (MQ=255)

               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                             <  2:2302264/93‑1 (MQ=255)

               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                             <  2:2999365/93‑1 (MQ=255)

                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                            >  1:737301/1‑93 (MQ=255)

                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                            >  2:1374384/1‑93 (MQ=255)

                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                            <  1:368577/93‑1 (MQ=255)

                 gagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTgg                                                                           >  2:2072987/1‑93 (MQ=255)

                  agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                          <  1:161756/93‑1 (MQ=255)

                  agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                          >  1:1820566/1‑93 (MQ=255)

                   gTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcg                                                                         >  1:2811179/1‑93 (MQ=255)

                   gTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcg                                                                         >  1:2518026/1‑93 (MQ=255)

                    tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                        <  2:765787/93‑1 (MQ=255)

                      aaaTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgcgc                                                                      >  2:3890443/1‑93 (MQ=255)

                      aaaTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgcgc                                                                      <  1:2782357/93‑1 (MQ=255)

                       aaTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCt                                                                     <  2:1471240/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                                 <  2:2687390/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                                 <  1:44128/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                                 <  2:3603206/93‑1 (MQ=255)

                             gcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATg                                                               <  1:3485864/93‑1 (MQ=255)

                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                             >  1:2367160/1‑93 (MQ=255)

                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                             >  2:22250/1‑93 (MQ=255)

                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                             >  2:1981635/1‑93 (MQ=255)

                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                             >  1:1341601/1‑93 (MQ=255)

                                  gTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCa                                                          <  2:1239909/93‑1 (MQ=255)

                                   tGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCAt                                                         <  1:3069948/93‑1 (MQ=255)

                                   tGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCAt                                                         >  1:2578012/1‑93 (MQ=255)

                                     cATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATcc                                                       <  1:3778352/93‑1 (MQ=255)

                                     cATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATcc                                                       >  2:2533613/1‑93 (MQ=255)

                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                     >  2:4082353/1‑93 (MQ=255)

                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                     <  1:1348652/93‑1 (MQ=255)

                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATATCTTTTAGCCGCTTCAGTTGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCac                                                     >  1:3792389/1‑92 (MQ=255)

                                           tACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAg                                                 >  1:568836/1‑93 (MQ=255)

                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                                <  1:2057092/93‑1 (MQ=255)

                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                                >  1:2640919/1‑93 (MQ=255)

                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                                >  2:2994697/1‑93 (MQ=255)

                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                               <  1:111024/93‑1 (MQ=255)

                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                               <  1:2182340/93‑1 (MQ=255)

                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                              >  2:757/1‑93 (MQ=255)

                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                             <  2:4375762/93‑1 (MQ=255)

                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                             >  2:195206/1‑93 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                           <  1:864871/93‑1 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                           >  1:689791/1‑93 (MQ=255)

                                                    cTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGc                                        <  1:2491709/93‑1 (MQ=255)

                                                     tGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcc                                       >  2:2692617/1‑93 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                      >  1:4657793/1‑93 (MQ=255)

                                                       tAGCCCGGCTGATTGATCTCTTGTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                     <  2:2463512/93‑1 (MQ=255)

                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                                 >  2:5420/1‑93 (MQ=255)

                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                                 >  2:5419/1‑93 (MQ=255)

                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                                <  1:838361/93‑1 (MQ=255)

                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                                <  1:3378801/93‑1 (MQ=255)

                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                              >  1:445793/1‑93 (MQ=255)

                                                               cTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAg                             >  2:2019366/1‑93 (MQ=255)

                                                               cTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAg                             <  1:126603/93‑1 (MQ=255)

                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                           >  1:3405587/1‑93 (MQ=255)

                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                           <  1:3108985/93‑1 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                          >  1:1076534/1‑93 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                          <  1:2776942/93‑1 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                          >  2:302671/1‑93 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                          >  2:302670/1‑93 (MQ=255)

                                                                   ttgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTc                         >  2:1115418/1‑93 (MQ=255)

                                                                   ttgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTc                         >  1:4319947/1‑93 (MQ=255)

                                                                   ttgatCTCTTTTAGCCGCTGCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTc                         <  2:2253908/93‑1 (MQ=255)

                                                                    tgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCg                        >  1:3552473/1‑93 (MQ=255)

                                                                     gatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGc                       >  1:4127838/1‑93 (MQ=255)

                                                                      atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                      >  2:1446677/1‑93 (MQ=255)

                                                                       tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                     >  2:4144359/1‑93 (MQ=255)

                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                    <  2:2494772/93‑1 (MQ=255)

                                                                         tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCACGGCGGACGGCCCCAGTGACAGTTTCGCCACg                   >  1:4139459/1‑93 (MQ=255)

                                                                         tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACg                   >  2:931221/1‑93 (MQ=255)

                                                                           ttttAGCCGCTTCAGTCGTGTGCGATACACCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc                 <  2:2506068/93‑1 (MQ=255)

                                                                           ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc                 <  1:341639/93‑1 (MQ=255)

                                                                           ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc                 >  1:4393521/1‑93 (MQ=255)

                                                                           ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc                 >  1:2273215/1‑93 (MQ=255)

                                                                             ttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCAt               <  2:4527825/93‑1 (MQ=255)

                                                                              tAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATa              <  2:2985985/93‑1 (MQ=255)

                                                                               aGCCGCTTCGGTCGTGTGCGATACAACGTGGCGCGCTCGAATGCCGCATCCAGCAGGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa             <  2:4212397/93‑1 (MQ=255)

                                                                               aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa             <  2:1085196/93‑1 (MQ=255)

                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc            <  1:4060516/93‑1 (MQ=255)

                                                                                   gCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaaa         <  2:4661048/93‑1 (MQ=255)

                                                                                    cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAg        >  2:3617144/1‑93 (MQ=255)

                                                                                    cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAg        <  2:3268861/93‑1 (MQ=255)

                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt       <  2:112212/93‑1 (MQ=255)

                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt       >  1:69340/1‑93 (MQ=255)

                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt       >  1:3349201/1‑93 (MQ=255)

                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGGGACAGTTTCGCCACGTCATAACAAAGt       >  1:3349200/1‑93 (MQ=255)

                                                                                      tCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGtt      >  2:346824/1‑93 (MQ=255)

                                                                                          ccGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAg  >  2:4328604/1‑93 (MQ=255)

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAG  >  NC_000913/1890135‑1890317

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