BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,470,983	T→C	Y95C (TAC→TGC)	yfdO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,470,983	0	T	C	99.3%	450.8 / ‑5.0	145	Y95C (TAC→TGC)	yfdO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (94/50); total (94/51)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 149 total aligned reads displayed.

TGCTCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGTACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGAA  >  NC_000913/2470895‑2471071

  cTCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGc                                                                                    <  2:3194231/93‑1 (MQ=255)

  cTCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGc                                                                                    >  1:4085742/1‑93 (MQ=255)

   tCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCt                                                                                   <  1:837430/93‑1 (MQ=255)

   tCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCt                                                                                   >  1:298906/1‑93 (MQ=255)

   tCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCt                                                                                   >  1:157993/1‑93 (MQ=255)

     aCCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTga                                                                                 >  1:4044718/1‑93 (MQ=255)

      cccccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgac                                                                                >  2:3548604/1‑93 (MQ=255)

       ccccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgccg                                                                               >  1:3065858/1‑93 (MQ=255)

       ccccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacg                                                                               <  2:1179177/93‑1 (MQ=255)

        cccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgg                                                                              >  1:941242/1‑93 (MQ=255)

        cccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgg                                                                              >  2:3175517/1‑93 (MQ=255)

         ccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgga                                                                             <  1:938986/93‑1 (MQ=255)

         ccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgga                                                                             <  1:843198/93‑1 (MQ=255)

         ccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgga                                                                             <  2:1988275/93‑1 (MQ=255)

         ccTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacgga                                                                             <  1:2524255/93‑1 (MQ=255)

          cTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTggcggac                                                                            >  1:4106010/1‑93 (MQ=255)

          cTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacggac                                                                            >  1:1781945/1‑93 (MQ=255)

          cTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacggac                                                                            >  1:3878663/1‑93 (MQ=255)

          cTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacggac                                                                            <  1:2113499/93‑1 (MQ=255)

           tGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTgacggacg                                                                           >  1:484435/1‑93 (MQ=255)

              aTCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCg                                                                        >  1:3274777/1‑93 (MQ=255)

              aTCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCg                                                                        >  2:1022159/1‑93 (MQ=255)

                      gggATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAgcgc                                                                >  1:2677798/1‑93 (MQ=255)

                      gggATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAgcgc                                                                <  2:4395134/93‑1 (MQ=255)

                       ggATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCa                                                               >  1:1564523/1‑93 (MQ=255)

                        gATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAg                                                              >  1:1293506/1‑93 (MQ=255)

                        gATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAg                                                              >  1:437728/1‑93 (MQ=255)

                         aTCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAgg                                                             >  2:2957370/1‑93 (MQ=255)

                         aTCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAgg                                                             >  2:2607784/1‑93 (MQ=255)

                          tCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAggg                                                            <  2:4638660/93‑1 (MQ=255)

                            tGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGaa                                                          <  2:4298342/93‑1 (MQ=255)

                             ggCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGaaa                                                         <  1:3282077/93‑1 (MQ=255)

                             ggCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGaaa                                                         <  2:2006606/93‑1 (MQ=255)

                               cTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGc                                                       >  2:3913147/1‑93 (MQ=255)

                               cTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGAGATAAGCGCAGGGAAAGc                                                       >  1:4095704/1‑93 (MQ=255)

                                tGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACTGACGGCGATAAGCGCAGGGAAAGCt                                                      >  2:4452829/1‑93 (MQ=255)

                                tGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGCAAGCt                                                      >  1:3150056/1‑93 (MQ=255)

                                  tAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCa                                                    >  2:714728/1‑93 (MQ=255)

                                  tAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCa                                                    >  2:1229999/1‑93 (MQ=255)

                                  tAGTCCACATTGTCGTAACTGGCTGTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCa                                                    <  2:2038141/93‑1 (MQ=255)

                                    gTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATc                                                  >  2:2268055/1‑93 (MQ=255)

                                     tCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATcc                                                 <  1:4287763/93‑1 (MQ=255)

                                        acaTTGTCGTAACTCGCTTTGAAGTACGGGACTTCGCGTTTTTCTGTGCACATGCTGACGGACGGCGATAAGCGCAGGGAGAGGTCATCCCAt                                              >  2:2553238/1‑93 (MQ=255)

                                         caTTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAtt                                             <  2:4200321/93‑1 (MQ=255)

                                         caTTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAtt                                             >  2:3567524/1‑93 (MQ=255)

                                         caTTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAtt                                             <  1:2418578/93‑1 (MQ=255)

                                         caTTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAtt                                             <  1:4676230/93‑1 (MQ=255)

                                          aTTGTCGTAACTGGCTTTGAAGTACGGGCCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAttt                                            >  2:4525050/1‑93 (MQ=255)

                                           ttGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAtttt                                           >  2:3382996/1‑93 (MQ=255)

                                            tGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCAttttt                                          >  1:1033443/1‑93 (MQ=255)

                                             gTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTc                                         >  1:4053036/1‑93 (MQ=255)

                                              tCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTcc                                        <  1:3016914/93‑1 (MQ=255)

                                              tCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTcc                                        >  2:2796737/1‑93 (MQ=255)

                                                gTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCg                                      >  2:4426691/1‑93 (MQ=255)

                                                 tAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGc                                     >  1:4592097/1‑93 (MQ=255)

                                                  aaCTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCa                                    >  2:178211/1‑93 (MQ=255)

                                                   aCTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAg                                   >  1:4423679/1‑93 (MQ=255)

                                                   aCTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAg                                   >  2:1576663/1‑93 (MQ=255)

                                                   aCTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAg                                   >  1:124807/1‑93 (MQ=255)

                                                     tGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCt                                 >  2:2727913/1‑93 (MQ=255)

                                                       gCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTc                               <  1:144303/93‑1 (MQ=255)

                                                       gCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTc                               <  1:674197/93‑1 (MQ=255)

                                                        cTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAGGCGCAGGGAAAGCTCATCCCATTTTTCCCGAAGCTTCg                              >  1:1737153/1‑93 (MQ=255)

                                                        cTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTca                              >  2:3313301/1‑92 (MQ=255)

                                                        cTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCg                              >  2:1680610/1‑93 (MQ=255)

                                                         tttGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGa                             <  1:1254468/93‑1 (MQ=255)

                                                         tttGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGa                             <  1:458984/93‑1 (MQ=255)

                                                         tttGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGa                             >  2:1305831/1‑93 (MQ=255)

                                                          ttGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGAc                            <  2:2666239/93‑1 (MQ=255)

                                                          ttGAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGAAGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGAc                            >  1:3489228/1‑93 (MQ=255)

                                                            gAAGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACgg                          >  1:3670137/1‑93 (MQ=255)

                                                             aaGTACGCGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACggg                         <  1:2683598/93‑1 (MQ=255)

                                                              aGTACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGc                        >  1:1327507/1‑93 (MQ=255)

                                                                tACGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTGCCCGCAGCTTCGACGGGCTg                      <  2:3059913/93‑1 (MQ=255)

                                                                 aCGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGa                     >  1:2047442/1‑93 (MQ=255)

                                                                 aCGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGa                     >  2:1813509/1‑93 (MQ=255)

                                                                  cGGGTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAg                    >  1:4417483/1‑93 (MQ=255)

                                                                   gggTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGc                   <  2:1501980/93‑1 (MQ=255)

                                                                    ggTCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCa                  >  2:2484413/1‑93 (MQ=255)

                                                                      tCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACg                >  2:4533105/1‑93 (MQ=255)

                                                                      tCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACg                >  1:1025501/1‑93 (MQ=255)

                                                                      tCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACg                >  2:828116/1‑93 (MQ=255)

                                                                      tCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACg                >  1:3709188/1‑93 (MQ=255)

                                                                      tCTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACg                >  2:2117457/1‑93 (MQ=255)

                                                                       cTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGt               >  1:1943993/1‑93 (MQ=255)

                                                                       cTTCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGt               <  1:3488222/93‑1 (MQ=255)

                                                                        ttCGCGTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGtt              >  2:589685/1‑93 (MQ=255)

                                                                          cgcgTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGCTAc            >  2:1016694/1‑93 (MQ=255)

                                                                            cgTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACgg          >  2:4257260/1‑93 (MQ=255)

                                                                             gTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGc         <  2:4173700/93‑1 (MQ=255)

                                                                             gTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGc         <  1:3364288/93‑1 (MQ=255)

                                                                             gTTTTTCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGc         >  2:3927815/1‑93 (MQ=255)

                                                                              tttttCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCa        >  2:2236353/1‑93 (MQ=255)

                                                                              tttttCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCa        >  1:3327523/1‑93 (MQ=255)

                                                                                tttCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCAcc      >  1:3938403/1‑93 (MQ=255)

                                                                                  tCTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACATTACGGCACCAg    <  1:1399265/93‑1 (MQ=255)

                                                                                   cTGTGCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGa   <  1:2573102/93‑1 (MQ=255)

                                                                                    tgtgTACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGaa  <  1:1092363/93‑1 (MQ=255)

                                                                                    tgtgCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGaa  <  2:4281247/93‑1 (MQ=255)

                                                                                    tgtgCACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGaa  <  1:3614338/93‑1 (MQ=255)

TGCTCACCCCCTGAATCCTGCCGGGATCTGGCTGTAGTCCACATTGTCGTAACTGGCTTTGAAGTACGGGTCTTCGCGTTTTTCTGTGTACGTGCTGACGGACGGCGATAAGCGCAGGGAAAGCTCATCCCATTTTTCCCGCAGCTTCGACGGGCTGAGCACGTTACGGCACCAGAA  >  NC_000913/2470895‑2471071

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