BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,204,258	T→C	W182R (TGG→CGG)	yehM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,204,258	0	T	C	99.2%	368.2 / ‑4.7	122	W182R (TGG→CGG)	yehM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (57/64); total (57/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 147 total aligned reads displayed.

CCCTCTGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCTGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGATGTGGCGC  >  NC_000913/2204168‑2204347

cccTCTGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcgg                                                                                         >  2:3207350/1‑93 (MQ=255)

 ccTCTGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcgga                                                                                        >  1:4231857/1‑93 (MQ=255)

  ctctGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggat                                                                                       <  1:3687649/93‑1 (MQ=255)

  ctctGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggat                                                                                       >  2:1456380/1‑93 (MQ=255)

   tctGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggatt                                                                                      >  2:1754644/1‑93 (MQ=255)

    ctGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggattg                                                                                     >  1:1755832/1‑93 (MQ=255)

      gcgcTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggattgga                                                                                   <  2:711622/93‑1 (MQ=255)

      gcgcTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggattgga                                                                                   <  2:1072458/93‑1 (MQ=255)

       cgcTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCcggattggat                                                                                  <  1:2599257/93‑1 (MQ=255)

        gcTGGCGCCCTATTTTGCCCAGCTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATg                                                                                 <  2:2725314/93‑1 (MQ=255)

        gcTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATg                                                                                 <  1:129534/93‑1 (MQ=255)

        gcTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATg                                                                                 <  2:1637530/93‑1 (MQ=255)

         cTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATgg                                                                                <  1:2828178/93‑1 (MQ=255)

         cTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATgg                                                                                >  1:152654/1‑93 (MQ=255)

          tGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATggg                                                                               <  2:825760/93‑1 (MQ=255)

           ggCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGc                                                                              <  2:2536260/93‑1 (MQ=255)

            gcgcACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCg                                                                             >  1:2124348/1‑93 (MQ=255)

            gcgcACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCg                                                                             >  1:2124349/1‑93 (MQ=255)

            gcgcACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCg                                                                             <  2:4204586/93‑1 (MQ=255)

             cgcACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGa                                                                            <  2:1718791/93‑1 (MQ=255)

              gcACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGAt                                                                           <  2:723924/93‑1 (MQ=255)

                ccTATTTTGCCCAACGGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgc                                                                         <  2:2232275/92‑1 (MQ=255)

                aCTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgc                                                                         <  2:542209/93‑1 (MQ=255)

                aCTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgc                                                                         >  1:781692/1‑93 (MQ=255)

                 cTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgca                                                                        >  2:896539/1‑93 (MQ=255)

                  tATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgcag                                                                       >  2:766880/1‑93 (MQ=255)

                  tATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgcag                                                                       <  2:2945248/93‑1 (MQ=255)

                  tATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgcag                                                                       >  1:3348118/1‑93 (MQ=255)

                   aTTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgcagc                                                                      <  1:1195799/93‑1 (MQ=255)

                   aTTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATgcagc                                                                      >  1:639772/1‑93 (MQ=255)

                    ttttGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGCTGGGAGATgcagca                                                                     >  2:1436002/1‑93 (MQ=255)

                     tttGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCTCCGGATTGGATGGGCGATGCAGCaa                                                                    >  2:4081652/1‑93 (MQ=255)

                       tGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCaaaa                                                                  <  2:3484874/93‑1 (MQ=255)

                        gCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaat                                                                 >  1:2057448/1‑93 (MQ=255)

                          ccAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataa                                                               <  1:1043837/93‑1 (MQ=255)

                          ccAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataa                                                               <  2:3760874/93‑1 (MQ=255)

                          ccAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataa                                                               <  1:853852/93‑1 (MQ=255)

                           cAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataat                                                              >  2:2287147/1‑93 (MQ=255)

                           cAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataat                                                              >  2:2169659/1‑93 (MQ=255)

                           cAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAaataat                                                              <  1:4101342/93‑1 (MQ=255)

                            aaCTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTTTGGCCCGCCGGATTGGATGGTCGATGCAGCAAAATAATg                                                             <  2:785607/93‑1 (MQ=255)

                            aaCTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATg                                                             <  1:1462825/93‑1 (MQ=255)

                            aaCTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATg                                                             <  2:1367908/93‑1 (MQ=255)

                             aCTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATgg                                                            <  1:3628642/93‑1 (MQ=255)

                               tGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCCCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCg                                                          <  1:2937171/93‑1 (MQ=255)

                                 cGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGGGATGCAGCAAAAAAATGGGGag                                                        >  2:3630899/1‑92 (MQ=255)

                                 cGGGGAGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGGTGCAGCAAAATAATGGCGAc                                                        >  2:3375160/1‑93 (MQ=255)

                                      ggccGCTTCCGGCGATGCGATTATTCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACgtgtg                                                   <  2:2925991/90‑2 (MQ=255)

                                       gACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGttt                                                  >  2:4082956/1‑92 (MQ=255)

                                       gACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTa                                                  >  2:2234263/1‑93 (MQ=255)

                                       gACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTa                                                  <  1:1964488/93‑1 (MQ=255)

                                        aCGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCTGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAg                                                 <  1:2943844/93‑1 (MQ=255)

                                        aCGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAg                                                 >  2:23272/1‑93 (MQ=255)

                                         cGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGt                                                <  2:2648132/93‑1 (MQ=255)

                                          gCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGtt                                               >  2:4473542/1‑93 (MQ=255)

                                          gCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGtt                                               >  2:2915492/1‑93 (MQ=255)

                                           ccTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTg                                              >  2:1380527/1‑93 (MQ=255)

                                            cTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGt                                             >  1:527922/1‑93 (MQ=255)

                                            cTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGt                                             >  2:3716840/1‑93 (MQ=255)

                                             tCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTc                                            >  1:2791762/1‑93 (MQ=255)

                                             tCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTc                                            <  1:4198101/93‑1 (MQ=255)

                                              ccGGCGATGCTCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCt                                           <  2:1319807/93‑1 (MQ=255)

                                              ccGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCt                                           <  2:1319806/93‑1 (MQ=255)

                                                ggCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGc                                         >  1:1804492/1‑93 (MQ=255)

                                                 gCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCg                                        >  2:4077972/1‑93 (MQ=255)

                                                  cGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCgg                                       <  1:3044327/93‑1 (MQ=255)

                                                   gATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGt                                      >  1:4410385/1‑93 (MQ=255)

                                                   gATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGt                                      >  2:213839/1‑93 (MQ=255)

                                                     tGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtgg                                    >  2:2777992/1‑93 (MQ=255)

                                                      gcgcTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggc                                   >  2:3383700/1‑93 (MQ=255)

                                                      gcgcTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGAATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggc                                   <  2:1653176/93‑1 (MQ=255)

                                                        gcTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggctg                                 >  1:1350984/1‑93 (MQ=255)

                                                         cTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggctgg                                <  1:946907/93‑1 (MQ=255)

                                                         cTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggctgg                                <  2:310085/93‑1 (MQ=255)

                                                         cTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGCCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggctgc                                >  1:982294/1‑92 (MQ=255)

                                                          tCAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGtggctggc                               <  2:4555449/93‑1 (MQ=255)

                                                           cAATCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCa                              >  1:3797259/1‑93 (MQ=255)

                                                            aaTCGTCAGCGCGAAGCCTTTATGGACCGCCGGATTGGAAGGGCGATTCAGCAAAATAATGGCGACGAGTTTGTTGTCTGCGTATGCTGGCAc                             >  2:4609553/1‑93 (MQ=255)

                                                             aTCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACg                            >  1:1170138/1‑93 (MQ=255)

                                                              tCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGc                           >  2:1610090/1‑93 (MQ=255)

                                                              tCGTCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGc                           <  2:2676359/93‑1 (MQ=255)

                                                                 tCAGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTcc                        <  1:2153045/93‑1 (MQ=255)

                                                                   aGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCgg                      <  1:3790770/93‑1 (MQ=255)

                                                                   aGCGCGAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCgc                      >  1:1871111/1‑92 (MQ=255)

                                                                     cgcgAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCa                    <  2:4282922/93‑1 (MQ=255)

                                                                        gAAGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTg                 <  1:2593262/93‑1 (MQ=255)

                                                                         aaGCCTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTgg                >  2:1949765/1‑93 (MQ=255)

                                                                             cTTTATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCaaa            >  2:3403375/1‑93 (MQ=255)

                                                                              tttATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAg           >  2:496611/1‑93 (MQ=255)

                                                                              tttATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAg           >  1:2025357/1‑93 (MQ=255)

                                                                               ttATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGa          <  2:1930036/93‑1 (MQ=255)

                                                                               ttATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGa          <  1:1689077/93‑1 (MQ=255)

                                                                               ttATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGa          >  1:4181341/1‑93 (MQ=255)

                                                                                tATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGAt         >  1:1640715/1‑93 (MQ=255)

                                                                                tATGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGAt         <  1:4101139/93‑1 (MQ=255)

                                                                                  tGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGAtgt       <  2:1915395/93‑1 (MQ=255)

                                                                                  tGGCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGAtgt       <  2:2057448/93‑1 (MQ=255)

                                                                                   ggCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGAtgtg      <  2:3763361/93‑1 (MQ=255)

                                                                                    gCCCGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCACCGGCACTGGCAAAGATGTgg     <  2:2986050/93‑1 (MQ=255)

                                                                                       cGCCGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGATGTGgcgc  >  2:4436682/1‑93 (MQ=255)

CCCTCTGCGCTGGCGCACTATTTTGCCCAACTGCGGGGCGACGCCTCCGGCGATGCGCTCAATCGTCAGCGCGAAGCCTTTATGGCCCGCTGGATTGGATGGGCGATGCAGCAAAATAATGGCGACGTGTTAGTTGTCTGCGGTGGCTGGCACGCTCCGGCACTGGCAAAGATGTGGCGC  >  NC_000913/2204168‑2204347

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