BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,525	(C)_6→7	coding (402/630 nt)	yqiE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,519	1	.	C	97.7%	701.5 / 7.7	219	coding (408/630 nt)	yqiE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (110/104); total (114/105)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 262 total aligned reads displayed.

CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACA  >  NC_000913/3175427‑3175606

 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑c                                                                                        <  1:5500559/93‑1 (MQ=255)

  aaTATCTTCGTTTTCATCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cc                                                                                       <  1:2986250/93‑1 (MQ=255)

  aaTATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCTTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cc                                                                                       <  2:6155217/93‑1 (MQ=255)

   atatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑ccc                                                                                      >  2:815037/1‑93 (MQ=255)

    tatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccc                                                                                     >  1:887740/1‑93 (MQ=255)

    tatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccc                                                                                     >  1:3864094/1‑93 (MQ=255)

      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccc                                                                                   >  1:5227984/1‑93 (MQ=255)

      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccc                                                                                   >  2:318320/1‑93 (MQ=255)

       cttcttTTCCACCACCCAGACCGTAAATACCGCTTGCGGTCGTGGCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑ccccccc                                                                                  <  1:454834/93‑2 (MQ=255)

        ttCGTTTTCATCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccccg                                                                                 <  2:1035824/93‑3 (MQ=255)

        ttCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccccg                                                                                 <  1:6441722/93‑3 (MQ=255)

           gTTTTCATCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑CCCCCCcgggc                                                                              <  1:1576614/93‑3 (MQ=255)

            ttttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCt                                                                              >  2:4387017/1‑93 (MQ=255)

            ttttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCt                                                                              <  2:2799198/93‑1 (MQ=255)

             tttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCtt                                                                             >  1:5563939/1‑93 (MQ=255)

              ttCATCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTg                                                                            <  1:3550215/93‑1 (MQ=255)

               tcatcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGc                                                                           >  2:2325999/1‑93 (MQ=255)

               tcatcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGc                                                                           >  1:2658635/1‑93 (MQ=255)

                  tcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAg                                                                        >  2:2396651/1‑93 (MQ=255)

                  acgcccAGACAGTAAATACCGCTGGCGGTCGTGGCGTCCACTTCCCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAg                                                                        <  1:2504129/89‑1 (MQ=255)

                    aGCCAGACGGTGAATCCGGCTTGCGGTCGTGGCGTCCACTCCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGa                                                                      <  2:2707712/93‑1 (MQ=255)

                    aGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGa                                                                      >  2:1428671/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGaa                                                                     >  2:1550378/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGaa                                                                     >  1:216528/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGaa                                                                     >  2:5645033/1‑93 (MQ=255)

                       cAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAAc                                                                   >  2:2932262/1‑93 (MQ=255)

                        aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACt                                                                  <  1:109071/93‑1 (MQ=255)

                         gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACtt                                                                 >  2:3236348/1‑93 (MQ=255)

                         gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACtt                                                                 >  1:720093/1‑93 (MQ=255)

                           ccGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTaa                                                               <  1:4067974/93‑1 (MQ=255)

                           ccGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTaa                                                               <  1:2091625/93‑1 (MQ=255)

                            cGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAac                                                              >  1:4141965/1‑93 (MQ=255)

                            cGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAac                                                              <  2:4802095/93‑1 (MQ=255)

                               gAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACAcc                                                           >  1:817133/1‑93 (MQ=255)

                               gAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACAcc                                                           >  1:5387174/1‑93 (MQ=255)

                                 aTACCGCTTGCGGTCGTGGCGTCTACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCgg                                                         >  1:5934316/1‑93 (MQ=255)

                                 aTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCgg                                                         >  1:6233247/1‑93 (MQ=255)

                                 aTACCGCTGGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCgg                                                         <  2:6307204/93‑1 (MQ=255)

                                 aTACCGCTGGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCgg                                                         <  1:417676/93‑1 (MQ=255)

                                  tcccGCTGGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCCCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                        <  2:4693264/91‑1 (MQ=255)

                                  tACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                        >  1:3704340/1‑93 (MQ=255)

                                   aCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGtt                                                       <  1:1851107/93‑1 (MQ=255)

                                    ccGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGttg                                                      >  2:6068012/1‑92 (MQ=255)

                                    ccGCTTGCGGTCGTGGCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGttt                                                      <  2:3534671/93‑1 (MQ=255)

                                     cGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTg                                                     <  1:1081426/93‑1 (MQ=255)

                                      gCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                    <  2:1663573/93‑1 (MQ=255)

                                      gCTTGCGGTCATGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                    <  2:6564294/93‑1 (MQ=255)

                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                   <  1:3800750/93‑1 (MQ=255)

                                        ttGCGGTCGTGGCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTc                                                  <  1:3490712/93‑1 (MQ=255)

                                          gCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTgggcg                                                 >  2:4744726/1‑90 (MQ=255)

                                             gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                             <  2:5011651/93‑1 (MQ=255)

                                              tCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                            <  1:4604905/93‑1 (MQ=255)

                                                   gCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                       >  2:885712/1‑93 (MQ=255)

                                                   gCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                       >  1:6346222/1‑93 (MQ=255)

                                                       ccccTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                   <  1:6630605/90‑1 (MQ=255)

                                                       ccACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                   <  2:6703886/93‑1 (MQ=255)

                                                        cACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTc                                  >  2:3856705/1‑93 (MQ=255)

                                                          cTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                >  2:2582660/1‑93 (MQ=255)

                                                          cTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                >  1:1188107/1‑93 (MQ=255)

                                                           ttCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                               >  2:6471208/1‑93 (MQ=255)

                                                           ttCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                               >  1:1737942/1‑93 (MQ=255)

                                                           ttCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                               >  2:1236023/1‑93 (MQ=255)

                                                           ttCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                               >  1:5019721/1‑93 (MQ=255)

                                                            tcccccACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                              <  2:525041/90‑1 (MQ=255)

                                                            tcccccACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                              <  1:4437822/90‑1 (MQ=255)

                                                            tCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                              >  2:5897730/1‑93 (MQ=255)

                                                             cGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGcc                             >  1:5712990/1‑93 (MQ=255)

                                                              gCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCct                            >  2:2619308/1‑93 (MQ=255)

                                                               cccACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctc                           <  1:613075/93‑1 (MQ=255)

                                                                ccaccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                          <  1:5769209/93‑1 (MQ=255)

                                                                ccaccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                          >  2:2486232/1‑93 (MQ=255)

                                                                ccaccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                          <  1:1352594/93‑1 (MQ=255)

                                                                ccaccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                          <  1:2335116/93‑1 (MQ=255)

                                                                 caccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                         <  1:3289399/93‑1 (MQ=255)

                                                                 caccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                         >  2:6630260/1‑93 (MQ=255)

                                                                  accaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTc                        >  2:515513/1‑93 (MQ=255)

                                                                    caTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACGCCGGTTTGGTCCGTTTGACTATCAGTCCCACCTCTTCaa                      <  1:4249124/93‑1 (MQ=255)

                                                                    caTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCaa                      >  2:4142988/1‑93 (MQ=255)

                                                                     aTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAAt                     >  1:2969020/1‑93 (MQ=255)

                                                                     aTAATTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAAt                     >  1:4028465/1‑93 (MQ=255)

                                                                       aaTTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                   <  2:921809/93‑1 (MQ=255)

                                                                        aTTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGc                  >  1:3173178/1‑93 (MQ=255)

                                                                        aTTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGc                  >  1:3968134/1‑93 (MQ=255)

                                                                        aTTGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCACCTCTTCAATCGc                  <  1:6601447/93‑1 (MQ=255)

                                                                         ttGACGACCGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCt                 <  2:3725129/93‑1 (MQ=255)

                                                                          tGACGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCtt                <  2:5731503/93‑1 (MQ=255)

                                                                            aCGAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcg              <  2:5431434/93‑1 (MQ=255)

                                                                              gAACGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgcg            >  2:2157530/1‑93 (MQ=255)

                                                                               aaCGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           <  1:5958915/93‑1 (MQ=255)

                                                                               aaCGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           >  2:1112905/1‑93 (MQ=255)

                                                                                aCGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGAc          <  1:2702845/93‑1 (MQ=255)

                                                                                aCGCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTACCAGTCCCGCCTCTTCAATCGCTTCGCGAc          <  2:4924095/93‑1 (MQ=255)

                                                                                  gCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACgg        <  2:4276258/93‑1 (MQ=255)

                                                                                  gCTCACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACgg        <  2:755320/93‑1 (MQ=255)

                                                                                   cTCCCTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACggg       <  1:3606126/93‑1 (MQ=255)

                                                                                    tCACTGGTCC‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGc      <  2:5489000/93‑1 (MQ=255)

                                                                                     cACTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCa     >  2:1851881/1‑93 (MQ=255)

                                                                                       cTGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc   >  2:5289470/1‑93 (MQ=255)

                                                                                        tGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTACCCCCGGTTTGGTCCGTTTGACTATCATTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACa  <  1:3899637/93‑1 (MQ=255)

                                                                                        tGGT‑GC‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACa  >  1:4837544/1‑93 (MQ=255)

CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACA  >  NC_000913/3175427‑3175606

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