BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	529,036	A→G	pseudogene (397/491 nt)	ybbD →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	529,036	0	A	G	96.6%	388.9 / 1.4	119	pseudogene (397/491 nt)	ybbD	pseudogene
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (58/56); total (61/58)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

Only 100 of 122 total aligned reads displayed.

GTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAA  >  NC_000913/528939‑529130

gTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTg                                                                                                 >  2:1396023/1‑100 (MQ=255)

 tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTg‑‑g                                                                                              <  1:2603666/100‑1 (MQ=255)

 tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTg‑‑g                                                                                              <  1:1695094/100‑1 (MQ=255)

 tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTg‑‑g                                                                                              <  2:2268191/100‑1 (MQ=255)

  gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGggtg‑‑gt                                                                                             >  1:2216416/1‑100 (MQ=255)

   aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGgtg‑‑gtg                                                                                            <  1:39299/100‑1 (MQ=255)

   aGGGCAAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGgtg‑‑gtg                                                                                            <  2:1136894/100‑1 (MQ=255)

    gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑gtgt                                                                                           <  1:1267229/100‑1 (MQ=255)

    gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑gtgt                                                                                           >  2:197104/1‑100 (MQ=255)

    gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATG‑‑gtgt                                                                                           >  1:1163178/1‑100 (MQ=255)

     ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTa                                                                                          >  1:2420589/1‑100 (MQ=255)

     ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTa                                                                                          <  2:918464/100‑1 (MQ=255)

      gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑gtgt                                                                                           <  2:916857/98‑1 (MQ=255)

      gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTaa                                                                                         <  2:365857/100‑1 (MQ=255)

       ccAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTaaa                                                                                        >  2:2778700/1‑100 (MQ=255)

       ccAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTaaa                                                                                        >  2:2682990/1‑100 (MQ=255)

        cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAAc                                                                                       <  1:1738530/100‑1 (MQ=255)

        cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAAc                                                                                       <  2:507576/100‑1 (MQ=255)

        cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAAc                                                                                       <  1:947761/100‑1 (MQ=255)

        cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATG‑‑GTGTAAAc                                                                                       <  2:1163167/100‑1 (MQ=255)

          gggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑t                                                                                    >  1:504102/1‑100 (MQ=255)

          gggAACAAGTGGCGAAAATCGTAGCAAAGAATGATCCAGAGACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑t                                                                                    >  2:2119443/1‑100 (MQ=255)

            gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑ttt                                                                                  >  2:753261/1‑100 (MQ=255)

               cgaGTGGCGAAAAACGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑tttttt                                                                               >  1:2055324/3‑100 (MQ=255)

               cAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATG‑‑GTGTAAACC‑tttttt                                                                               >  1:1795540/1‑100 (MQ=255)

                aaGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑ttttttt                                                                              <  1:16630/100‑1 (MQ=255)

                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑tttttttt                                                                             <  2:1472702/100‑1 (MQ=255)

                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑tttttttt                                                                             >  2:929953/1‑100 (MQ=255)

                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑tttttttt                                                                             >  1:1232682/1‑100 (MQ=255)

                  gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttg                                                                            <  1:26959/100‑1 (MQ=255)

                  gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttg                                                                            <  2:368359/100‑1 (MQ=255)

                   tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttgt                                                                           <  1:1261761/100‑1 (MQ=255)

                   tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttgt                                                                           <  2:1428741/100‑1 (MQ=255)

                    ggCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttgtt                                                                          >  2:804114/1‑100 (MQ=255)

                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTgg                                                                        >  1:1474419/1‑100 (MQ=255)

                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTgg                                                                        >  2:2907197/1‑100 (MQ=255)

                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTgg                                                                        >  1:1883067/1‑100 (MQ=255)

                       gAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGt                                                                       <  1:1227441/100‑1 (MQ=255)

                        aaaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTg                                                                      <  2:1444894/100‑1 (MQ=255)

                        aaaaTCGTATAAAAAAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGTTACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTttgtt                                                                          <  2:1483601/96‑1 (MQ=255)

                           aTCGTATAAAAGAATAACACACACACAAAAAAAGTGTGGTGTAACTATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGaaa                                                                   <  2:1008213/100‑1 (MQ=255)

                            tCGTATTAAAGAATGAAAAAAAAAAAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAAt                                                                  <  2:479831/100‑1 (MQ=255)

                            tCGTATCAAAGAATGATCCAGATAC‑AAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAAtt                                                                 <  1:2852899/100‑1 (MQ=255)

                             cGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAAtt                                                                 <  2:316475/100‑1 (MQ=255)

                              gTATCAAAAAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGTTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGaaan                                                                  <  2:1005306/98‑2 (MQ=255)

                               tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTaa                                                               <  2:633517/100‑1 (MQ=255)

                                aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAAt                                                              <  1:1849544/100‑1 (MQ=255)

                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                           <  2:2668376/100‑1 (MQ=255)

                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                           >  2:375291/1‑100 (MQ=255)

                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                           >  1:2866653/1‑100 (MQ=255)

                                     aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTAc                                                         >  2:2385104/1‑100 (MQ=255)

                                          gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                    >  2:468749/1‑100 (MQ=255)

                                          gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                    >  2:2861749/1‑100 (MQ=255)

                                          gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                    >  1:2836139/1‑100 (MQ=255)

                                          gATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                    >  2:1609597/1‑100 (MQ=255)

                                            tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATTTTACGCattat                                                  >  2:2638049/1‑100 (MQ=255)

                                              cAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                >  1:694729/1‑100 (MQ=255)

                                                gagaCAAAAAAAGTAGGGTGTAAATATCGTAAGATACCAGGGGAAGAGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAAGTAATGTTACGCATTATTTta                                              <  2:1261670/97‑1 (MQ=255)

                                                gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTTTTTTTGTTATTAATGTTACACATTATTTta                                              >  2:3000044/1‑100 (MQ=255)

                                                gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                              <  1:2292633/100‑1 (MQ=255)

                                                gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATATTACGCATTATTTta                                              >  1:2982795/1‑100 (MQ=255)

                                                  tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtata                                            >  1:1254856/1‑100 (MQ=255)

                                                    cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtata                                            >  2:1192442/1‑100 (MQ=255)

                                                    cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGAGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                          >  1:2139413/1‑100 (MQ=255)

                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                         >  2:331714/1‑100 (MQ=255)

                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                         <  2:295539/100‑1 (MQ=255)

                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTGTGTTGGTGAAATTAGTGTTACGCATTATTTTATAACa                                         >  2:2043495/1‑100 (MQ=255)

                                                       aaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaaa                                       >  2:472355/1‑100 (MQ=255)

                                                               tggtgTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAg                               <  2:1882824/100‑1 (MQ=255)

                                                               tggtgTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAg                               >  1:2097109/1‑100 (MQ=255)

                                                                ggtgTAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAga                              >  1:1882841/1‑99 (MQ=255)

                                                                 gtgtAACTATGGTAAGATACAAGGGCAAGGGTGTG‑‑GTATAAATC‑TTTTTTTTGATGATGAAATTAATGTTACGCATTATTTTATACCACATATTGGa                                <  2:1860646/97‑1 (MQ=255)

                                                                 gtgtAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                             >  1:495450/1‑100 (MQ=255)

                                                                   gtAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTgn                           <  2:694723/100‑2 (MQ=255)

                                                                    tAAATATGGTAAGATACCAGGGCAAGGGTGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                          <  2:91848/100‑1 (MQ=255)

                                                                    tAAATATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                          >  1:1437176/1‑100 (MQ=255)

                                                                      aaTATGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAtt                        >  1:1390771/1‑100 (MQ=255)

                                                                       atatGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTg                       <  1:1780276/100‑1 (MQ=255)

                                                                         atGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                     <  2:11121/100‑1 (MQ=255)

                                                                          tGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACCTTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                     >  1:1233716/1‑100 (MQ=255)

                                                                          tGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                    <  2:1899829/100‑1 (MQ=255)

                                                                          tGGTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                    <  1:1245918/100‑1 (MQ=255)

                                                                           ggTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                   >  1:436515/1‑100 (MQ=255)

                                                                            gTAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa                  <  2:197884/100‑1 (MQ=255)

                                                                            gTAAGATACCAGGGCAAGGGGGTG‑‑GTATAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa                  >  1:1852631/1‑100 (MQ=255)

                                                                             tAAGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGAt                 <  2:2749532/100‑1 (MQ=255)

                                                                              aaGATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATg                <  1:1396035/100‑1 (MQ=255)

                                                                                gATACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCt              >  1:1697110/1‑100 (MQ=255)

                                                                                  tACCAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg            >  2:1648985/1‑100 (MQ=255)

                                                                                     cAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGc         <  2:1474407/100‑1 (MQ=255)

                                                                                     cAGGGCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATACTTGTGc         <  1:1599711/100‑1 (MQ=255)

                                                                                       gggCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTTTTGTTTAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag       >  2:1106209/1‑100 (MQ=255)

                                                                                       gggCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag       <  1:753268/100‑1 (MQ=255)

                                                                                       gggCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag       >  2:101407/1‑100 (MQ=255)

                                                                                        gtcAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga      <  1:2870428/98‑1 (MQ=255)

                                                                                         gCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag     <  1:481946/100‑1 (MQ=255)

                                                                                         gCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag     >  2:1792027/1‑100 (MQ=255)

                                                                                         gCAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag     <  2:1663984/100‑1 (MQ=255)

                                                                                         gCAAGGGGATG‑‑GTGTAAACC‑TTTTTTTTGTTTGTTAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag     >  2:852689/1‑100 (MQ=255)

                                                                                          cAAGGGGGTG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGt    <  2:1018707/100‑1 (MQ=255)

GTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATG‑‑GTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAA  >  NC_000913/528939‑529130

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