BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,356	0	T	C	49.1%	37.7 / inf	1035	S178S (AGT→AGC)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (245/279); new base C (241/264); total (492/543)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-08
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1070 total aligned reads displayed.

GAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949261‑1949451

  aTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTAt                                                                                               <  2:11713551/96‑1 (MQ=255)

     cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCaa                                                                                            <  2:4143927/96‑1 (MQ=255)

       ccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAcc                                                                                          <  2:3535/96‑1 (MQ=255)

        cgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCa                                                                                         <  1:8227805/96‑1 (MQ=255)

         gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                        <  2:9810474/96‑1 (MQ=255)

         gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATAGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                        >  2:12625713/1‑96 (MQ=255)

           ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                      <  1:1496676/96‑1 (MQ=255)

            cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                                     <  1:7602711/96‑1 (MQ=255)

             gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCt                                                                                    >  1:470450/1‑96 (MQ=255)

              cgaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                   >  1:9197809/3‑96 (MQ=255)

              cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                                   <  2:7599508/96‑1 (MQ=255)

              cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGACGCCACCGCCGAGCAACACAATAACAGAATTAATCATATCTa                                                                                   >  1:6866290/1‑96 (MQ=255)

               aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAc                                                                                  >  2:4939957/1‑96 (MQ=255)

               aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAc                                                                                  >  1:11548979/1‑96 (MQ=255)

               aaTGCCTGGGAACTTGGCTTTAATCTGGGGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACACTAACAGCATTAATCATATCTAc                                                                                  >  2:11111358/1‑96 (MQ=255)

                aTGCCTTGGAACTTGGTTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCATCCATATCTAcc                                                                                 <  1:6598663/96‑1 (MQ=255)

                 tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                <  1:11712045/96‑1 (MQ=255)

                  ggcTGGGCACTTGGCTTTAATCCGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATACCAGTATCAACCATATCTACCcg                                                                               <  1:3245918/94‑1 (MQ=255)

                  gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               >  2:8510356/1‑96 (MQ=255)

                  gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               >  2:1960699/1‑96 (MQ=255)

                  gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                               <  2:2991326/96‑1 (MQ=255)

                    cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcg                                                                             <  2:9518131/96‑1 (MQ=255)

                        gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                         >  1:6428618/1‑96 (MQ=255)

                        gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATc                                                                         >  2:4013007/1‑96 (MQ=255)

                        gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATc                                                                         >  2:1656078/1‑96 (MQ=255)

                         aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                        <  1:4936198/96‑1 (MQ=255)

                         aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                        >  2:4737372/1‑96 (MQ=255)

                            ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                     <  1:2776923/96‑1 (MQ=255)

                            ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                     <  1:990147/96‑1 (MQ=255)

                            ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                                     >  2:3194656/1‑96 (MQ=255)

                            tgggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                     >  1:10737192/3‑96 (MQ=255)

                              ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGcccat                                                                   <  1:14260281/96‑3 (MQ=255)

                              ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgt                                                                   >  1:1965035/1‑96 (MQ=255)

                                cTTTAATCTGGTGATTGCCGAAGACGCCTGTCGTGCCGCCAGCGCCGAGCAGCCCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                                 <  1:5966025/96‑1 (MQ=255)

                                  ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGCGTGc                                                               <  2:2555117/96‑1 (MQ=255)

                                  ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                               <  2:10745465/96‑1 (MQ=255)

                                   tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCg                                                              >  1:10992304/1‑96 (MQ=255)

                                      tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                           >  2:10535895/1‑96 (MQ=255)

                                      tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                           >  2:4526824/1‑96 (MQ=255)

                                       cTGATGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                          >  2:10107942/1‑96 (MQ=255)

                                         ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                        <  1:7380971/96‑1 (MQ=255)

                                         ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                        <  2:8809297/96‑1 (MQ=255)

                                          gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                       <  1:3496904/96‑1 (MQ=255)

                                          gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                       >  2:7891584/1‑96 (MQ=255)

                                            gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                     >  2:13647225/1‑96 (MQ=255)

                                             aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                    <  1:8067978/96‑1 (MQ=255)

                                               tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                                  <  2:2317402/96‑1 (MQ=255)

                                                 ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                >  2:2080962/1‑96 (MQ=255)

                                                 ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                >  1:12578533/1‑96 (MQ=255)

                                                  cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                               >  1:5442353/1‑96 (MQ=255)

                                                   gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAATAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                                              >  1:2569217/1‑96 (MQ=255)

                                                   gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                                              >  2:12989333/1‑96 (MQ=255)

                                                   gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                <  2:6578427/94‑1 (MQ=255)

                                                    aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             <  1:3910188/96‑1 (MQ=255)

                                                    aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                             <  2:5209054/96‑1 (MQ=255)

                                                     aGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                            >  2:5951911/1‑96 (MQ=255)

                                                      gACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTc                                           >  1:9527475/1‑96 (MQ=255)

                                                      gACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTc                                           <  2:5806962/96‑1 (MQ=255)

                                                      gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCtt                                           <  1:7855080/96‑2 (MQ=255)

                                                         gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                        <  2:12151990/96‑1 (MQ=255)

                                                         gCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                                        <  2:5560503/96‑5 (MQ=255)

                                                          ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTTGCGTGGAAGAGATCCTCAAcg                                       <  2:10267061/96‑1 (MQ=255)

                                                           atGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                                      <  1:226381/95‑1 (MQ=255)

                                                            tGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgcg                                     >  1:13694753/1‑96 (MQ=255)

                                                             gTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGt                                    <  1:8576499/96‑1 (MQ=255)

                                                                gTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                 <  2:11330489/96‑1 (MQ=255)

                                                                 tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATg                                <  1:14009218/96‑1 (MQ=255)

                                                                   ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                              <  1:11314241/96‑1 (MQ=255)

                                                                        aGCGCCGAGCAACACAATAACAGCATTACTCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                         <  1:12244061/96‑1 (MQ=255)

                                                                         gcgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                        <  1:5251340/96‑1 (MQ=255)

                                                                          cgcCGAGCAGCACAATAACAGTATCAACCATATCTGCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                       <  1:7128977/96‑1 (MQ=255)

                                                                          cgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                       <  1:2976510/96‑1 (MQ=255)

                                                                           gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                      >  2:13996606/1‑96 (MQ=255)

                                                                           gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGGAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                      >  1:10010070/1‑96 (MQ=255)

                                                                           gcCGAGCAGCACAATAACAGTAGCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                      <  2:6964606/96‑1 (MQ=255)

                                                                           gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                      >  1:3716099/1‑96 (MQ=255)

                                                                            ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg                     <  1:4171207/96‑1 (MQ=255)

                                                                            ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg                     >  2:6260580/1‑96 (MQ=255)

                                                                             cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                    >  1:5684518/1‑96 (MQ=255)

                                                                             cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                    >  2:376853/1‑96 (MQ=255)

                                                                              gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc                   <  2:1454308/96‑1 (MQ=255)

                                                                              gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc                   >  2:11837315/1‑96 (MQ=255)

                                                                               agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCATAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGtgt                  <  1:10578607/96‑3 (MQ=255)

                                                                                gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                 <  2:14180879/96‑2 (MQ=255)

                                                                                gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                 <  2:9287916/96‑2 (MQ=255)

                                                                                 cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTGACATCGGTCTAc                >  1:10315604/1‑96 (MQ=255)

                                                                                   gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg              <  2:13344430/96‑1 (MQ=255)

                                                                                   gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg              <  2:8149701/96‑1 (MQ=255)

                                                                                    cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGc             <  1:8374078/96‑1 (MQ=255)

                                                                                     acaATAACCGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCa            >  1:12971574/1‑96 (MQ=255)

                                                                                     acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa            <  1:12011062/96‑1 (MQ=255)

                                                                                     acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa            <  1:4528355/96‑1 (MQ=255)

                                                                                     acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa            >  1:1662547/1‑96 (MQ=255)

                                                                                      gaaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa           >  1:8822989/2‑96 (MQ=255)

                                                                                         tAACAGTATCAACCATATCTACCCGCGCGTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgg        >  2:13919874/1‑96 (MQ=255)

                                                                                         tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgg        >  2:2780431/1‑96 (MQ=255)

                                                                                          aaCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATAGGTCTACCGCAATGGt       >  1:12772221/1‑96 (MQ=255)

                                                                                          aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGt       <  1:6139697/96‑1 (MQ=255)

                                                                                           aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc      <  1:3121752/96‑1 (MQ=255)

                                                                                           aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTc      <  1:13417154/96‑1 (MQ=255)

GAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT  >  NC_000913/1949261‑1949451

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`