BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,410	0	C	T	45.7%	216.0 / inf	541	L185L (CTC→CTT)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand): ref base C (162/132); new base T (125/122); total (287/254)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 553 total aligned reads displayed.

CCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCAC  >  NC_000913/1949319‑1949502

   gTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgc                                                                                         >  1:1023651/1‑94 (MQ=255)

    tAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgca                                                                                        >  1:3504432/1‑93 (MQ=255)

         gcgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATg                                                                                   <  1:7266999/93‑1 (MQ=255)

            ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                                                                                <  1:5578967/94‑1 (MQ=255)

             cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                                                                           <  2:4054100/82‑1 (MQ=255)

              aGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                                                                          >  2:3185657/1‑82 (MQ=255)

               gcgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                                                                                         <  2:7130646/82‑1 (MQ=255)

                 gcCGAGCAACACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGt                                                                                       >  2:5730948/1‑82 (MQ=255)

                  ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                      >  2:4330851/1‑82 (MQ=255)

                  ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                      <  2:1300180/82‑1 (MQ=255)

                   cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATcc                                                                         >  1:6447902/1‑93 (MQ=255)

                   cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATca                                                                         >  1:7123400/1‑93 (MQ=255)

                     agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCggg                                                                       >  1:2931191/1‑93 (MQ=255)

                     agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCggc                                                                       >  1:3190639/1‑93 (MQ=255)

                      gcagcaCAATAACAGTATCAACCATATCTACCCGCGCAGCGCCCGTGTGCGTAGCGTGGAAGATATCCTCAACGCGTTATGATTTACATCGGTc                                                                      >  1:7465167/1‑94 (MQ=255)

                       caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                                                                                 <  2:1578841/79‑1 (MQ=255)

                       caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                                                                                 <  2:2958544/79‑1 (MQ=255)

                        agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCtt                                                                    >  1:148550/1‑93 (MQ=255)

                          cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                                                              <  2:736912/82‑1 (MQ=255)

                            caATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                                                                            <  2:6316026/82‑1 (MQ=255)

                            caATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCgg                                                                >  1:3043603/1‑93 (MQ=255)

                            caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCga                                                                >  1:4101795/1‑93 (MQ=255)

                             aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                                                                           <  2:2150322/82‑1 (MQ=255)

                              gaaaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                              <  1:334012/92‑1 (MQ=255)

                              aTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                                                                          >  2:849915/1‑82 (MQ=255)

                              aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCac                                                              >  1:2484802/1‑93 (MQ=255)

                               tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCGCAACGCGTTCTGATTTATCTCGATCTACCGaaac                                                             >  1:2062793/1‑93 (MQ=255)

                                 aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt                                                                       >  2:3692755/1‑82 (MQ=255)

                                 aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgc                                                           >  1:627960/1‑93 (MQ=255)

                                  gaGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGt                                                          <  1:4534600/93‑1 (MQ=255)

                                  cAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc                                                                      >  2:398491/1‑82 (MQ=255)

                                   aGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCt                                                                     >  2:6088348/1‑82 (MQ=255)

                                    gTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTa                                                                    <  2:2621359/82‑1 (MQ=255)

                                     tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAc                                                                   <  2:6746134/82‑1 (MQ=255)

                                      aTCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGct                                                      >  1:1416868/1‑93 (MQ=255)

                                        taatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGc                                                                <  2:3550388/78‑1 (MQ=255)

                                        cAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                                                                <  2:5797381/82‑1 (MQ=255)

                                         aatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATcc                                                   <  1:2183048/91‑1 (MQ=255)

                                          ttcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCt                                                  <  1:998464/92‑1 (MQ=255)

                                          atgatatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                              >  2:510507/4‑82 (MQ=255)

                                          aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcca                                                  >  1:5097563/1‑93 (MQ=255)

                                           ccATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAAt                                                             <  2:4477428/82‑1 (MQ=255)

                                            cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTag                                                >  1:1855297/1‑93 (MQ=255)

                                            cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATg                                                            >  2:1457290/1‑82 (MQ=255)

                                             ctatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                               <  1:6923486/93‑1 (MQ=255)

                                             atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaag                                               >  1:566201/1‑93 (MQ=255)

                                             atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgg                                                           >  2:4666803/1‑82 (MQ=255)

                                               atatACCCGCGCATCGCCCGTGAGCGTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTc                                                         >  2:1769811/4‑82 (MQ=255)

                                               atCTTCCCGCGCATCGACAGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc                                                         <  2:2268788/82‑1 (MQ=255)

                                                tCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCg                                                        >  2:190274/1‑82 (MQ=255)

                                                   aCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAt                                                     >  2:3335160/1‑82 (MQ=255)

                                                       gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTa                                                 <  2:6898995/82‑1 (MQ=255)

                                                       gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGt                                     <  1:4018942/94‑1 (MQ=255)

                                                       ccgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTGATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                     <  1:4256227/93‑1 (MQ=255)

                                                         gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTg                                   <  1:1183403/94‑1 (MQ=255)

                                                         ccATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTg                                   <  1:6674730/93‑1 (MQ=255)

                                                           aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATg                                             <  2:4801157/82‑1 (MQ=255)

                                                              gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggtag                              >  1:3891646/1‑91 (MQ=255)

                                                              gCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGt                                          <  2:2720459/82‑1 (MQ=255)

                                                              gCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGt                                          <  2:6989426/82‑1 (MQ=255)

                                                               tccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                             <  1:5502195/93‑1 (MQ=255)

                                                               cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                         >  2:115582/1‑82 (MQ=255)

                                                               cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                         >  2:3093522/1‑82 (MQ=255)

                                                               cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTAta                             >  1:6584508/1‑93 (MQ=255)

                                                                ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                        >  2:5403472/1‑82 (MQ=255)

                                                                ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATcg                            >  1:3616184/1‑93 (MQ=255)

                                                                ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTT‑CATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCat                           >  1:6116633/1‑93 (MQ=255)

                                                                ccGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGc                                        >  2:747761/1‑82 (MQ=255)

                                                                 cGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCg                                       <  2:5033567/82‑1 (MQ=255)

                                                                   tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                         >  1:4712343/1‑94 (MQ=255)

                                                                   tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCa                         >  1:1626387/1‑94 (MQ=255)

                                                                      gCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                  >  2:2039177/1‑82 (MQ=255)

                                                                        gTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCtc                    >  1:5799560/1‑93 (MQ=255)

                                                                         tAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTg                   >  1:5321542/1‑94 (MQ=255)

                                                                           gCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGat                 >  1:3726520/1‑93 (MQ=255)

                                                                           acGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGaa                 <  1:25739/93‑1 (MQ=255)

                                                                            cGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCa                            <  2:4939477/82‑1 (MQ=255)

                                                                             ttGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAaga               <  1:846807/93‑1 (MQ=255)

                                                                             gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGAAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga               >  1:1588790/1‑94 (MQ=255)

                                                                             gTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCAc                           <  2:2368846/82‑1 (MQ=255)

                                                                              cggAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagag              <  1:6360723/93‑1 (MQ=255)

                                                                               tgAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt             <  1:542856/93‑1 (MQ=255)

                                                                               ggAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCa                         >  2:4570123/1‑82 (MQ=255)

                                                                                gAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                        >  2:6001852/1‑82 (MQ=255)

                                                                                aaaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa            <  1:5194046/93‑1 (MQ=255)

                                                                                 gagagATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTAt           <  1:6797652/93‑1 (MQ=255)

                                                                                  agagATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcc                      <  2:7065877/82‑1 (MQ=255)

                                                                                  agagATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTAtc          >  1:4836209/1‑93 (MQ=255)

                                                                                   gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                     <  2:1812887/82‑1 (MQ=255)

                                                                                   gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                     <  2:6524952/82‑1 (MQ=255)

                                                                                   gagaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCt                     >  2:4232893/1‑82 (MQ=255)

                                                                                   aagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc         <  1:4338070/93‑1 (MQ=255)

                                                                                   aagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc         <  1:3663196/93‑1 (MQ=255)

                                                                                     caTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccc       <  1:5925558/93‑1 (MQ=255)

                                                                                      aTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGa                  <  2:2543691/82‑1 (MQ=255)

                                                                                      aTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGa                  >  2:5238100/1‑82 (MQ=255)

                                                                                      aTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGccca      >  1:2781840/1‑93 (MQ=255)

                                                                                       tCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGaa                 <  2:286777/82‑1 (MQ=255)

                                                                                        ccTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGcc    >  1:2293425/1‑94 (MQ=255)

                                                                                           cAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt             <  2:6833574/82‑1 (MQ=255)

CCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCAC  >  NC_000913/1949319‑1949502

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