BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,806	(T)_8→9	coding (27/633 nt)	ycfQ ←	repressor for bhsA(ycfR)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,798	1	.	T	95.2%	339.3 / 9.1	105	coding (35/633 nt)	ycfQ	repressor for bhsA(ycfR)
Reads supporting (aligned to +/- strand): ref base . (4/1); new base T (51/48); total (55/50)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 121 total aligned reads displayed.

AGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCG  >  NC_000913/1168699‑1168890

 ggTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑t                                                                                              >  2:2313666/1‑100 (MQ=255)

   tCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt                                                                                            >  2:1031145/1‑100 (MQ=255)

    ccGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑t                                                                                          >  2:1704742/1‑100 (MQ=255)

    ccGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑t                                                                                          >  2:2315268/1‑100 (MQ=255)

       cAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑tttt                                                                                       <  1:299377/100‑1 (MQ=255)

        aaGAGAAGTCGCTTCATAACCGTGTTTCCAGAACAATGTCATCGCCTTATCAAGCGCGGCATCCCTGTCGAACCCTTTTGGGCGGCCACGGC‑‑ttt‑ttttt                                                                                      <  1:2430297/100‑1 (MQ=255)

        aaGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑ttttt                                                                                      <  2:829352/100‑1 (MQ=255)

        aaGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑ttttt                                                                                      <  1:619342/100‑1 (MQ=255)

          gagaAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑tttttta                                                                                    <  1:2682316/100‑3 (MQ=255)

          gagaAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑tttttta                                                                                    <  1:11764/100‑3 (MQ=255)

          gagaAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑tttttt                                                                                     <  1:2078594/99‑2 (MQ=255)

            gaAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑ttt‑ttttttaca                                                                                  <  2:2982648/100‑5 (MQ=255)

              aGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTAcaca                                                                                 >  2:328255/1‑100 (MQ=255)

               gTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACAt                                                                                >  2:1945729/1‑100 (MQ=255)

               gTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACAt                                                                                >  1:1794554/1‑100 (MQ=255)

                  gCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACAtt                                                                               >  2:2628280/1‑98 (MQ=255)

                  gCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATtgt                                                                             >  2:2201310/1‑100 (MQ=255)

                  gCTTCAGAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATtg                                                                              <  2:2087064/99‑1 (MQ=255)

                   cTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATtgtg                                                                            >  1:1293092/1‑100 (MQ=255)

                     tCATAACCGTGTTGCCAGAACAATTTCATGGCTTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGtt                                                                          <  1:1990714/100‑1 (MQ=255)

                     tCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGtt                                                                          <  2:1293079/100‑1 (MQ=255)

                     tCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGtt                                                                          <  1:991098/100‑1 (MQ=255)

                     tCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGtt                                                                          <  1:2467716/100‑1 (MQ=255)

                      cATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTg                                                                         >  2:1791333/1‑100 (MQ=255)

                       aTAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTg                                                                         >  2:2380183/1‑99 (MQ=255)

                       aTAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGa                                                                        >  2:146704/1‑100 (MQ=255)

                       aTAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGa                                                                        >  2:146498/1‑100 (MQ=255)

                        tAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAg                                                                       <  1:2294984/100‑1 (MQ=255)

                        tAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAg                                                                       <  1:2270832/100‑1 (MQ=255)

                        tAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGCGTTGa                                                                        >  2:2864183/1‑99 (MQ=255)

                        naaCCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAg                                                                       >  1:1795114/2‑100 (MQ=255)

                          aCCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTc                                                                     <  1:357808/100‑1 (MQ=255)

                           ccGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTc                                                                     <  1:2926428/99‑1 (MQ=255)

                             gtgtTGCCAGAACAGTTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAg                                                                       <  2:132150/95‑1 (MQ=255)

                             gtgtTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTATACATTGTGTTGAGTCAGt                                                                  >  2:2518265/1‑100 (MQ=255)

                             gtgtTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGt                                                                  >  1:1730599/1‑100 (MQ=255)

                                ttGCCAGAACAAGTTCATGGCCTGATCAAGCGCGGCATCTCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGc                                                               <  2:2633626/100‑1 (MQ=255)

                                  gCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGcc                                                              >  2:291057/1‑99 (MQ=255)

                                  gCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGcc                                                              <  1:639236/99‑1 (MQ=255)

                                  gCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCa                                                             >  1:506634/1‑100 (MQ=255)

                                   ccAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCat                                                            <  2:1550006/100‑1 (MQ=255)

                                   ccAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGGCAGTTGCCat                                                            >  2:2090626/1‑100 (MQ=255)

                                    cAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCatg                                                           >  2:1481893/1‑100 (MQ=255)

                                    cAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCatg                                                           >  1:362131/1‑100 (MQ=255)

                                    cAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCatg                                                           >  1:303322/1‑100 (MQ=255)

                                    cAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCatg                                                           <  1:2920106/100‑1 (MQ=255)

                                    cAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑TTT‑TTTTTAAACATTGTGTTGAGTCAGTTGCCatga                                                          >  2:1710655/1‑100 (MQ=255)

                                      gAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAgatcggaagag                                                         >  1:2651745/1‑90 (MQ=255)

                                        aCAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCcatcat                                                         >  2:1969348/1‑98 (MQ=255)

                                         cAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGcc                                                      >  1:2812028/1‑100 (MQ=255)

                                          aaTTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCg                                                     >  1:1403166/1‑100 (MQ=255)

                                           agttCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGt                                                    <  2:553154/98‑1 (MQ=255)

                                           aTTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGt                                                    <  1:1676324/100‑1 (MQ=255)

                                            nttCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGtt                                                   >  1:2670930/2‑100 (MQ=255)

                                                   gCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTg                                            <  2:1617999/100‑1 (MQ=255)

                                                   gCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTg                                            <  2:1423601/100‑1 (MQ=255)

                                                   gCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTg                                            >  1:2156326/1‑100 (MQ=255)

                                                    ccTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTgg                                           >  1:341037/1‑100 (MQ=255)

                                                    ccTTATCAACCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTgg                                           <  2:79753/100‑1 (MQ=255)

                                                       tATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTg                                            >  1:2234960/1‑97 (MQ=255)

                                                        aTCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGAc                                       >  1:1059072/1‑100 (MQ=255)

                                                        aTCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGAc                                       <  2:2233537/100‑1 (MQ=255)

                                                        aTCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGAc                                       >  1:1936729/1‑100 (MQ=255)

                                                          cAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACtg                                     >  1:1038425/1‑100 (MQ=255)

                                                          cAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACtg                                     >  1:1522212/1‑100 (MQ=255)

                                                            aGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACtgtg                                   >  1:1085476/1‑100 (MQ=255)

                                                             gcgcGGAATCCCTGTCGAACACTTTTGGGCGGCCCCGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACtgt                                    <  1:1925338/98‑1 (MQ=255)

                                                              cgcgGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGtgaa                                 <  1:146499/100‑1 (MQ=255)

                                                               gcgGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGtgaat                                <  2:2254704/100‑1 (MQ=255)

                                                                cgGGATCCATGTCGTATACTTTTCGGCGGCCACGGCA‑TTT‑TTTTTACACATTGTGTTGAGTCAGTAGCCATGATGCCGTTGTACCTGGGGAc                                       <  2:554888/92‑1 (MQ=255)

                                                                cgGCTTCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGtgaa                                 <  1:146705/98‑1 (MQ=255)

                                                                     tctttGTCGAATATTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAg                           <  2:1364527/95‑1 (MQ=255)

                                                                     tCCCTGTCGAACATTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAg                           <  2:913572/99‑1 (MQ=255)

                                                                     tCCCTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGCTGCCGTTGTACCTGGTGACTGTGAATGAAAgg                          >  2:1318301/1‑100 (MQ=255)

                                                                     acctcgTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAg                           <  2:2430275/94‑1 (MQ=255)

                                                                      cccTATCTAACACTTTTGGGCGGCCACGGC‑‑TTTATTTTTACTCATTGTGTTGAGTCAAATGCCATGAACCCTTTCTACCTGATTTCTGTTAATGAAAGGt                         >  2:1817632/1‑100 (MQ=255)

                                                                       ccTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtt                        >  1:391739/1‑100 (MQ=255)

                                                                       ccTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtt                        >  1:2599324/1‑100 (MQ=255)

                                                                       ccTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtt                        <  2:1920063/100‑1 (MQ=255)

                                                                       ccTGTCGAACACTTTTGGGCGGCCACGGC‑‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtta                       >  1:949618/1‑100 (MQ=255)

                                                                        cTGTCGAACCCTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtta                       <  1:1791348/100‑1 (MQ=255)

                                                                        cTGTCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGtta                       >  1:347529/1‑100 (MQ=255)

                                                                         tGTCGAACACGTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGttat                      <  1:2315288/100‑1 (MQ=255)

                                                                           tCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGttatta                    <  2:391734/100‑1 (MQ=255)

                                                                           tCGAACACTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGttatta                    >  1:878427/1‑100 (MQ=255)

                                                                               acacTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATaaa                <  1:23382/100‑1 (MQ=255)

                                                                                  cTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaat             >  1:2703460/1‑100 (MQ=255)

                                                                                  cTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaat             >  2:704792/1‑100 (MQ=255)

                                                                                  cTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaat             >  2:2964912/1‑100 (MQ=255)

                                                                                  cTTTTGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaat             <  1:675444/100‑1 (MQ=255)

                                                                                  cTTTTGGGCGGCCACGGC‑‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaata            >  2:130643/1‑100 (MQ=255)

                                                                                   ttttGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaata            <  2:506628/100‑1 (MQ=255)

                                                                                     ttGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAaataat          <  2:1702850/100‑1 (MQ=255)

                                                                                      tGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATc         <  1:2090644/100‑1 (MQ=255)

                                                                                      tGGGCGGCCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATc         <  2:1038415/100‑1 (MQ=255)

                                                                                        ggcggcCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATTAAAGGTTATTATAAAAATAATCAc       >  2:754570/1‑100 (MQ=255)

                                                                                         gcggcCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCAcc      <  2:451224/100‑1 (MQ=255)

                                                                                            gcCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTcc   >  1:1546301/1‑100 (MQ=255)

                                                                                            gcCACGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTcc   >  1:563240/1‑100 (MQ=255)

                                                                                               aCGGCT‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCg  >  2:588845/1‑98 (MQ=255)

AGGTCCGCAAGAGAAGTCGCTTCATAACCGTGTTGCCAGAACAATTTCATGGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGC‑‑TTT‑TTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCG  >  NC_000913/1168699‑1168890

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