BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,829,419	C→T	R446C (CGC→TGC)	yicH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,829,419	0	C	T	99.4%	546.0 / ‑5.5	159	R446C (CGC→TGC)	yicH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (77/81); total (77/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 173 total aligned reads displayed.

ACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC  >  NC_000913/3829329‑3829508

aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGc                                                                                         >  2:2182736/1‑93 (MQ=255)

  gCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGccc                                                                                       <  1:4163176/93‑1 (MQ=255)

   cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCg                                                                                      >  2:3914775/1‑93 (MQ=255)

    tagTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGt                                                                                     <  2:2815007/91‑1 (MQ=255)

     ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTc                                                                                    >  1:1272160/1‑93 (MQ=255)

     ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTc                                                                                    <  2:3719808/93‑1 (MQ=255)

       ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTTCCCGttgt                                                                                  >  2:3569090/1‑90 (MQ=255)

        tACCGATCATAAATGGGGGGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCt                                                                                 <  2:1928186/93‑1 (MQ=255)

        tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCt                                                                                 <  2:730623/93‑1 (MQ=255)

         ccccaTCATAAATGGGGCGTCTGGCGTGGCTCGCCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATCTTCGTTGCCCGCCACTg                                                                                <  2:3976460/89‑1 (MQ=255)

         aCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTTGTTGCCCGTCGCTg                                                                                >  2:1012746/1‑93 (MQ=255)

         aCCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTg                                                                                <  1:3500824/93‑1 (MQ=255)

          ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTgg                                                                               <  1:1477681/93‑1 (MQ=255)

          ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTgg                                                                               >  2:1980720/1‑93 (MQ=255)

          ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTgg                                                                               <  1:2705292/93‑1 (MQ=255)

           cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGgc                                                                              >  2:2905583/1‑93 (MQ=255)

            gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTTTTGATGTTCGTTGCCCGGCGCTGgcc                                                                             >  1:289847/1‑92 (MQ=255)

            gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGgcg                                                                             >  1:1735066/1‑93 (MQ=255)

            gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGGCGCTGgcg                                                                             >  2:1370752/1‑93 (MQ=255)

             aTCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCTTTGCCCGTCGCTGgcgc                                                                            >  1:4248861/1‑93 (MQ=255)

              tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCt                                                                           >  2:450203/1‑93 (MQ=255)

              tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCt                                                                           <  2:1456604/93‑1 (MQ=255)

               cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTg                                                                          <  2:438274/93‑1 (MQ=255)

                aTAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGa                                                                         <  1:907454/93‑1 (MQ=255)

                 tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAc                                                                        >  2:2579956/1‑93 (MQ=255)

                 tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGAc                                                                        >  2:3527061/1‑93 (MQ=255)

                    aTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGGCGCCGGCGAGACATTCAATCGTGTTGATGTTTGTTGATCGTCGCTGGCGCTGACCGc                                                                     >  2:3133192/1‑93 (MQ=255)

                    aTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGc                                                                     >  2:3877763/1‑93 (MQ=255)

                         gCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCcaca                                                                >  2:187775/1‑93 (MQ=255)

                           gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                              <  1:2136337/93‑1 (MQ=255)

                           gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                              <  1:4414354/93‑1 (MQ=255)

                            tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagca                                                             <  1:4408099/93‑1 (MQ=255)

                            tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAAcagca                                                             >  1:582371/1‑93 (MQ=255)

                              tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                                           >  2:15464/1‑93 (MQ=255)

                              tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCCCGACATTCAATCGTGTTGATGTTCGTTGACCGTCGCTGGCGCTGACCGCCAACagcagc                                                           >  2:4685030/1‑93 (MQ=255)

                               ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAgcagca                                                          <  2:688004/93‑1 (MQ=255)

                                gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCAc                                                         <  1:3901298/93‑1 (MQ=255)

                                  gTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                       >  2:2348730/1‑93 (MQ=255)

                                       tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                  <  2:2491675/93‑1 (MQ=255)

                                       tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                  >  2:4570334/1‑93 (MQ=255)

                                       tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                                                  >  2:1231592/1‑93 (MQ=255)

                                        cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                                                 >  1:2430300/1‑93 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                               <  2:3440569/93‑1 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                               >  1:713950/1‑93 (MQ=255)

                                          gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                                               >  1:894771/1‑93 (MQ=255)

                                            gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                                             >  2:4730053/1‑93 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAgg                                            >  1:3461192/1‑92 (MQ=255)

                                             aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                                            >  2:2265998/1‑93 (MQ=255)

                                                  gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                                       >  2:1254261/1‑93 (MQ=255)

                                                  gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                                       >  1:3444479/1‑93 (MQ=255)

                                                    aTGCCGCCGCCGCGACATTCAATCGTGTTTATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                     <  2:3725580/93‑1 (MQ=255)

                                                    aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                     <  2:3810869/93‑1 (MQ=255)

                                                    aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                                     <  1:1809772/93‑1 (MQ=255)

                                                     tGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                                    >  1:3587401/1‑93 (MQ=255)

                                                     tGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGATGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                                    <  1:1994843/93‑1 (MQ=255)

                                                      gccgccgccGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                   <  1:3130780/93‑1 (MQ=255)

                                                      gccgccgccGCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                   <  1:1843540/93‑1 (MQ=255)

                                                       ccgccgccgCGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                  >  1:4279524/1‑93 (MQ=255)

                                                          ccgccgcGACATTCACTCGTGTTGATTTTCGTTGCCCGTCGCGGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAt                               <  2:301744/93‑1 (MQ=255)

                                                          ccgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAt                               >  2:3913505/1‑93 (MQ=255)

                                                           cgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAtt                              <  2:1445245/93‑1 (MQ=255)

                                                           cgccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAtt                              <  1:4154704/93‑1 (MQ=255)

                                                            gccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt                             >  2:1047043/1‑93 (MQ=255)

                                                             ccgcGACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa                            <  2:1890552/93‑1 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGGACGGTGAATATCAGCGAACTGAGTGCATTTACt                          >  1:1203586/1‑93 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGCATATCAGCGAACTGAGCGCCTTTACt                          >  2:840065/1‑93 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                          <  1:2812920/93‑1 (MQ=255)

                                                               gcgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt                          >  1:2534819/1‑93 (MQ=255)

                                                                cgACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg                         <  2:4304663/93‑1 (MQ=255)

                                                                 gACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGa                        >  1:1844726/1‑93 (MQ=255)

                                                                 gACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGa                        >  2:2109055/1‑93 (MQ=255)

                                                                 gACATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCCAACTGAGAGTAATTACtca                        >  1:2593176/1‑91 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  2:2870218/1‑93 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  1:1054168/1‑93 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  1:3688419/1‑93 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  2:4321573/1‑93 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  1:2310705/1‑93 (MQ=255)

                                                                  aCATTCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaa                       >  1:2678245/1‑93 (MQ=255)

                                                                     ttCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGaaaaa                    >  2:3666656/1‑93 (MQ=255)

                                                                     ttCAATCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCCAACTGAGTGCATTTACTGaaaaa                    >  2:3245141/1‑93 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 >  2:943591/1‑93 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 >  2:2873629/1‑93 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 <  2:4504826/93‑1 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 <  2:499062/93‑1 (MQ=255)

                                                                        aaTCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                 <  2:588330/93‑1 (MQ=255)

                                                                          tCGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAt               <  1:3163987/93‑1 (MQ=255)

                                                                           cGTGTTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAtt              <  1:519222/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCGCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc             <  1:2786964/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc             <  1:2608133/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc             <  2:162221/93‑1 (MQ=255)

                                                                            gtgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc             <  2:2105581/93‑1 (MQ=255)

                                                                             tgtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCt            <  2:4128744/93‑1 (MQ=255)

                                                                              gtTGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTg           >  1:787633/1‑93 (MQ=255)

                                                                               ttGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg          <  1:2661212/93‑1 (MQ=255)

                                                                               ttGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg          >  1:2733086/1‑93 (MQ=255)

                                                                               ttGATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg          >  2:798701/1‑93 (MQ=255)

                                                                                tGATCTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGa         <  1:2540856/93‑1 (MQ=255)

                                                                                 gATGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGaa        <  1:4740739/93‑1 (MQ=255)

                                                                                  aTGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAg       >  1:1481516/1‑93 (MQ=255)

                                                                                   tGTTCGTTGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGc      >  1:2034629/1‑93 (MQ=255)

ACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC  >  NC_000913/3829329‑3829508

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