BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,958,877	2 bp→CG	coding (586‑587/2022 nt)	rep →	rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,958,877	0	G	C	100.0%	1132.4 / NA	355	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (183/172); total (183/172)
*	NC_000913	3,958,878	0	C	G	100.0%	1138.6 / NA	356	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (183/173); total (183/173)

Only 100 of 362 total aligned reads displayed.

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964

||

   tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCAACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                 >  2:3788241/1‑93 (MQ=255)

    aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                <  1:4097074/93‑1 (MQ=255)

    aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                <  1:5934378/93‑1 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCTATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  1:4193083/93‑1 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  1:2162542/93‑1 (MQ=255)

      gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                              <  2:2016427/93‑1 (MQ=255)

      gggCTGTATGATGCACACCGGAATGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                              <  1:4788036/93‑1 (MQ=255)

       ggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGaagaag                                                                             <  1:6921164/93‑1 (MQ=255)

        gCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGt                                                                            <  2:5286108/93‑1 (MQ=255)

         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  2:3154447/93‑1 (MQ=255)

         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  2:7309364/93‑1 (MQ=255)

         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  2:5187409/93‑1 (MQ=255)

             atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                       >  2:1080048/1‑93 (MQ=255)

             atgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCa                                                                       <  1:728774/93‑1 (MQ=255)

               gatgCCCACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  2:4914135/93‑1 (MQ=255)

               gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  1:6687069/93‑1 (MQ=255)

                atgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgc                                                                    <  1:1353156/93‑1 (MQ=255)

                 tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                   >  1:5766996/1‑93 (MQ=255)

                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                  >  2:6873609/1‑93 (MQ=255)

                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                  <  2:3314233/93‑1 (MQ=255)

                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                >  1:4484696/1‑93 (MQ=255)

                    acacCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTg                                                                <  2:6357491/93‑1 (MQ=255)

                        cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                            >  2:5466728/1‑93 (MQ=255)

                        cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                            >  1:510712/1‑93 (MQ=255)

                        cTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAg                                                            >  2:6299818/1‑93 (MQ=255)

                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           >  1:7061629/1‑93 (MQ=255)

                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           >  1:2811068/1‑93 (MQ=255)

                           aaaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAATAAGTCCGCAAGCGCTGGCAGAAc                                                         >  2:4346081/1‑93 (MQ=255)

                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                        <  2:6371965/93‑1 (MQ=255)

                               ccTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaaa                                                     >  1:6646226/1‑93 (MQ=255)

                                cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                                    >  1:792591/1‑93 (MQ=255)

                                cTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAt                                                    >  1:6155460/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   >  2:715126/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   >  2:1586959/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   >  1:7211984/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   <  2:7861432/93‑1 (MQ=255)

                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  <  1:6733065/93‑1 (MQ=255)

                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  <  1:2616942/93‑1 (MQ=255)

                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  >  2:4517384/1‑93 (MQ=255)

                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  <  1:4694288/93‑1 (MQ=255)

                                   tAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCg                                                 >  2:385793/1‑93 (MQ=255)

                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                <  2:604270/93‑1 (MQ=255)

                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                >  2:5567450/1‑93 (MQ=255)

                                       gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                             <  2:890028/93‑1 (MQ=255)

                                       gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                             <  1:5457181/93‑1 (MQ=255)

                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                            >  1:5095454/1‑93 (MQ=255)

                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                            <  2:1852093/93‑1 (MQ=255)

                                         tctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATct                                           <  2:4026949/93‑1 (MQ=255)

                                         tctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATct                                           <  1:2491972/93‑1 (MQ=255)

                                          ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                                          >  2:4571913/1‑93 (MQ=255)

                                           tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                         >  2:1117492/1‑93 (MQ=255)

                                            cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                        <  2:6480747/93‑1 (MQ=255)

                                            cGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgct                                        <  2:2937120/93‑1 (MQ=255)

                                             gACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgctg                                       >  1:2790058/1‑93 (MQ=255)

                                              aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                                      >  2:2181635/1‑93 (MQ=255)

                                              aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                                      >  2:2617129/1‑93 (MQ=255)

                                              aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                                      <  1:2958567/93‑1 (MQ=255)

                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                     >  1:1579969/1‑93 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  2:6611706/1‑93 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    <  2:5854278/93‑1 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    <  1:533405/93‑1 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  2:2347478/1‑93 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  2:3244857/1‑93 (MQ=255)

                                                 tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                   >  1:3951001/1‑93 (MQ=255)

                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                  <  1:4988244/93‑1 (MQ=255)

                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                  <  2:7601302/93‑1 (MQ=255)

                                                     tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCCCTATCTGCTGGTGGATGa                               >  2:1374382/1‑93 (MQ=255)

                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                              >  1:3489775/1‑93 (MQ=255)

                                                        tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                            >  1:7355402/1‑93 (MQ=255)

                                                         cTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTAt                           >  2:227764/1‑93 (MQ=255)

                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                          >  1:3270248/1‑93 (MQ=255)

                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCAGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTAGTGGATGAGTATc                          >  1:3118108/1‑93 (MQ=255)

                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                         <  2:3551867/93‑1 (MQ=255)

                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  2:3479351/1‑93 (MQ=255)

                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       <  2:6986158/93‑1 (MQ=255)

                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  1:6378146/1‑93 (MQ=255)

                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       >  1:3671412/1‑93 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAggc                      >  1:359068/1‑92 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      <  1:2443126/93‑1 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      <  1:7669414/93‑1 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      >  1:7587526/1‑93 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      <  1:400068/93‑1 (MQ=255)

                                                              tttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      <  2:5231806/93‑1 (MQ=255)

                                                                 tttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                   <  1:2218050/92‑1 (MQ=255)

                                                                 aTTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGCCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                   <  2:1510133/93‑1 (MQ=255)

                                                                     ccGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaac               <  2:820829/93‑1 (MQ=255)

                                                                     ccGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaac               >  1:151644/1‑93 (MQ=255)

                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac             <  2:5704242/93‑1 (MQ=255)

                                                                        aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc            >  2:7490567/1‑93 (MQ=255)

                                                                          gTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag          <  1:6271150/93‑1 (MQ=255)

                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  2:2877991/1‑93 (MQ=255)

                                                                              ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag      >  2:497607/1‑93 (MQ=255)

                                                                              ctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagccag      >  1:1194000/1‑93 (MQ=255)

                                                                               tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt     >  1:5914321/1‑93 (MQ=255)

                                                                               tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt     >  2:615427/1‑93 (MQ=255)

                                                                               tgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGt     >  2:5074849/1‑93 (MQ=255)

                                                                                gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa    <  2:173629/93‑1 (MQ=255)

                                                                                gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa    >  2:1138595/1‑93 (MQ=255)

                                                                                 ctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCTGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAt   <  2:6573402/93‑1 (MQ=255)

                                                                                 ctgcAACGCAATGAAGAAGCCCGCAAGCGCGGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTAt   <  1:182953/93‑1 (MQ=255)

||

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964

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