BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,156,591	0	T	G	22.9%	389.8 / 136.1	210	intergenic (+73/‑32)	yqhD/dkgA	aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand): ref base T (74/88); new base G (30/18); total (104/106)
Fisher's exact test for biased strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 217 total aligned reads displayed.

GCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCT  >  NC_000913/3156493‑3156689

gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCt                                                                                                    <  1:1411420/99‑1 (MQ=255)

 ccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCCGACTTAAGTTCACAACACCTCACCGGAGCCtc                                                                                                   >  1:4199508/1‑98 (MQ=255)

      tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCg                                                                                              <  1:1126990/99‑1 (MQ=255)

       ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                             >  2:1900796/1‑99 (MQ=255)

       ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                             <  2:4470598/99‑1 (MQ=255)

       ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                             <  2:2462036/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCGGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGt                                                                                            <  1:1216417/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:1341800/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:368873/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:1727711/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:193303/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:3621371/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:2020638/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:3274238/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                            <  1:3930171/99‑1 (MQ=255)

        taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGt                                                                                            <  1:4951980/99‑1 (MQ=255)

         aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCCCCGGAGCCTGCCCCGGTg                                                                                           >  2:4026604/1‑99 (MQ=255)

         aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           >  1:3072962/1‑99 (MQ=255)

         aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           >  1:1717358/1‑99 (MQ=255)

          cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGa                                                                                          >  2:3692856/1‑99 (MQ=255)

          cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                          <  2:3501877/99‑1 (MQ=255)

          cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                          <  2:1283485/99‑1 (MQ=255)

           gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                         >  1:2164432/1‑99 (MQ=255)

           gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                         <  2:2961233/99‑1 (MQ=255)

           gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTCCGTTTCCGGGCATTTCGTCCAGCCTTAAGTCCACAACACCCCCCCGGACCCTCCCCCGGTGAg                                                                                         >  2:107968/1‑99 (MQ=255)

            aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                        >  1:2580994/1‑99 (MQ=255)

             aGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggt                                                                                       <  2:1781799/99‑3 (MQ=255)

              gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCCGACTTAAGCTCACAACACCTCACCGGAGCCTGCCCCGGTGAGTTc                                                                                      >  2:3868241/1‑99 (MQ=255)

                cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    <  2:255808/99‑1 (MQ=255)

                cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    <  1:5094910/99‑1 (MQ=255)

                cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    <  2:900240/99‑1 (MQ=255)

                  ccggcTAAGCTTTTTACGTCCCAAACTTTCGTTTTCGCGCATTTCGACCAGATCGAAGGTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCagat                                                                                  <  2:1597663/96‑1 (MQ=255)

                  cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCatat                                                                                  <  2:1209939/99‑1 (MQ=255)

                   ccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                 >  2:1542704/1‑99 (MQ=255)

                    cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTCAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaa                                                                                <  2:2387827/99‑1 (MQ=255)

                     gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                               <  2:3161760/99‑1 (MQ=255)

                      cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                              <  1:3663249/99‑1 (MQ=255)

                      cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                              >  2:2635525/1‑99 (MQ=255)

                        aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                            <  2:2206003/99‑1 (MQ=255)

                           cTTTTTACGCCTCAAACTTTAGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                         >  2:3796594/1‑99 (MQ=255)

                                tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTa                                                                    <  1:3503969/99‑1 (MQ=255)

                                  cGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATACAGGAGGAACGTATg                                                                  >  2:2814985/1‑99 (MQ=255)

                                    ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                <  1:712197/99‑1 (MQ=255)

                                     cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                               >  1:2882302/1‑99 (MQ=255)

                                      tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTGACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                              >  2:4589534/1‑99 (MQ=255)

                                      tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                              >  1:4515316/1‑99 (MQ=255)

                                       cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGGGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                             <  2:289262/99‑1 (MQ=255)

                                        aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                            <  2:1830316/99‑1 (MQ=255)

                                            tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                        >  1:665187/1‑99 (MQ=255)

                                            tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                        >  1:2662969/1‑99 (MQ=255)

                                            tttCGTTTTCGGGAATTTCGTCCAGCCTTAAGTTCACAACACCTCAACGGAGCATGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTACTCCat                                                        >  2:4734543/1‑98 (MQ=255)

                                                 ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   >  2:338310/1‑99 (MQ=255)

                                                 ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   <  1:402966/99‑1 (MQ=255)

                                                 ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   >  1:4338235/1‑99 (MQ=255)

                                                        gCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGCGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                            <  1:3834435/99‑1 (MQ=255)

                                                         cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTCTGGCTAATCCAACCGTTATTAAGCt                                           >  1:41877/1‑99 (MQ=255)

                                                          aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTa                                          <  2:198459/99‑1 (MQ=255)

                                                               gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGa                                     >  1:2274098/1‑99 (MQ=255)

                                                               gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGa                                     >  1:4770551/1‑99 (MQ=255)

                                                               gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGa                                     <  1:3749684/99‑1 (MQ=255)

                                                                tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                    >  1:3472700/1‑99 (MQ=255)

                                                                 ctaGACTTAAGTTCACACCACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                   <  1:235804/97‑1 (MQ=255)

                                                                     aCTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                               <  2:3627525/99‑1 (MQ=255)

                                                                       ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                             >  1:4587173/1‑99 (MQ=255)

                                                                       ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                             <  2:105052/99‑1 (MQ=255)

                                                                           gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                         >  2:3862804/1‑99 (MQ=255)

                                                                            ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        <  1:3445592/99‑1 (MQ=255)

                                                                            ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        >  2:2880957/1‑99 (MQ=255)

                                                                            ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        >  2:1998520/1‑99 (MQ=255)

                                                                               acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                     >  2:2087619/1‑99 (MQ=255)

                                                                               acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                     <  1:4202319/99‑1 (MQ=255)

                                                                               acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                     >  1:3049140/1‑99 (MQ=255)

                                                                               acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTAATGccc                     >  2:4817567/1‑99 (MQ=255)

                                                                                 aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                   >  2:2619882/1‑99 (MQ=255)

                                                                                  tgaacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAg                  >  2:743247/5‑99 (MQ=39)

                                                                                   cacCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                 >  2:4116699/1‑99 (MQ=255)

                                                                                      cgcACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGGAATGTCATGCCCCAGCTgg              <  1:4050922/97‑1 (MQ=255)

                                                                                      cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              <  1:2447877/99‑1 (MQ=255)

                                                                                      cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              <  2:733917/99‑1 (MQ=255)

                                                                                      cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              <  1:4100021/99‑1 (MQ=255)

                                                                                      cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              <  2:2231477/99‑1 (MQ=255)

                                                                                      cTCACCGGAGCCTGCGCCGGTGAGTTCAGATAGAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGATACAGGATGGCAATGTCATGCCCCAGCTgg              <  2:345290/99‑1 (MQ=255)

                                                                                       tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  2:3755205/1‑99 (MQ=255)

                                                                                       tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:4658710/1‑99 (MQ=255)

                                                                                       tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:1449337/1‑99 (MQ=255)

                                                                                       tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:3967320/1‑99 (MQ=255)

                                                                                       tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  2:555486/1‑99 (MQ=255)

                                                                                       tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  2:3219645/1‑99 (MQ=255)

                                                                                         aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCACCCGCTATTAAGCTACAGGATGGCAATGTCACGCCCCCGCTcggcc           >  2:1067887/1‑97 (MQ=255)

                                                                                         aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:100731/1‑99 (MQ=255)

                                                                                         aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:3122485/1‑99 (MQ=255)

                                                                                         aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:554379/1‑99 (MQ=255)

                                                                                         aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAAAGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  2:3894869/1‑99 (MQ=255)

                                                                                          ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt          >  1:4695178/1‑99 (MQ=255)

                                                                                          ccGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt          >  1:1966279/1‑99 (MQ=255)

                                                                                            ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg        <  1:2723952/99‑1 (MQ=255)

                                                                                              aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc      <  1:4488151/99‑1 (MQ=255)

                                                                                                  tgctgcGGTGAGTTCGTATAAAGGGGGGGCGGAGGGCTAAGCCAACCGTTATTAGGCTGCAGGATGGGGATGTCATGCCCCAGCTGGGACTGGGCGTcg  <  1:1468590/99‑2 (MQ=255)

                                                                                                  tGCTCCGGTGAGTTCATATAAAGGGGGAAAGTTTGGCTAATCCAATCGTTATTAAGGTACAGGATGGCAATGTCATCCCCCCGCTGGCACTGGGCGGCt  >  2:495940/1‑99 (MQ=255)

                                                                                                  tGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt  >  1:4870059/1‑99 (MQ=255)

GCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCT  >  NC_000913/3156493‑3156689

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