BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,250,165	C→T	R74R (CGC→CGT)	yhaH →	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,250,165	0	C	T	88.5%	565.5 / 41.1	191	R74R (CGC→CGT)	yhaH	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (6/16); new base T (45/124); total (51/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 196 total aligned reads displayed.

TTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAG  >  NC_000913/3250073‑3250255

ttAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCcgccgc                                                                                            >  1:2306553/1‑93 (MQ=255)

  aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCt                                                                                          <  2:6143968/93‑1 (MQ=255)

     cTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCa                                                                                       <  2:6112054/93‑1 (MQ=255)

       ggCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACg                                                                                     >  1:3130486/1‑93 (MQ=255)

         cAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCATGac                                                                                   >  2:659016/1‑93 (MQ=255)

         cAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGac                                                                                   >  2:4826490/1‑93 (MQ=255)

          aaCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGaca                                                                                  <  1:1713017/93‑1 (MQ=255)

          aaCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGaca                                                                                  >  1:1935255/1‑93 (MQ=255)

           aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGacac                                                                                 >  1:1629215/1‑93 (MQ=255)

           aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGacac                                                                                 <  1:4729051/93‑1 (MQ=255)

           aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGacac                                                                                 <  1:4893655/93‑1 (MQ=255)

            cgcgCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCTGCCGTCTGCACGACacc                                                                                >  2:4263324/1‑93 (MQ=255)

              cgcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccga                                                                              <  2:1819613/93‑1 (MQ=255)

               gcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgac                                                                             <  1:5987927/93‑1 (MQ=255)

               gcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTCTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgac                                                                             <  1:6298725/93‑1 (MQ=255)

                  ggtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGc                                                                          <  2:7388564/93‑1 (MQ=255)

                  ggtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGc                                                                          >  2:6952583/1‑93 (MQ=255)

                    tggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTc                                                                        >  2:3676430/1‑93 (MQ=255)

                     ggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCg                                                                       <  2:191673/93‑1 (MQ=255)

                       cGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGc                                                                     <  2:7347272/93‑1 (MQ=255)

                       cGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGc                                                                     >  1:7235017/1‑93 (MQ=255)

                         aaGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCgt                                                                   <  1:5815382/93‑1 (MQ=255)

                           ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                 >  2:5003857/1‑93 (MQ=255)

                           ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                 <  1:3714288/93‑1 (MQ=255)

                           ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                 <  1:1311091/93‑1 (MQ=255)

                           ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                 >  1:4131863/1‑93 (MQ=255)

                           ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                 <  1:1156666/93‑1 (MQ=255)

                            gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                                <  1:5028755/93‑1 (MQ=255)

                            gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                                <  1:1108605/93‑1 (MQ=255)

                            gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                                <  1:4365447/93‑1 (MQ=255)

                            gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                                <  1:2541921/93‑1 (MQ=255)

                            gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                                >  2:6194903/1‑93 (MQ=255)

                             cATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggtg                                                               <  1:6462237/93‑1 (MQ=255)

                              aTCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGtggtgg                                                              >  2:1103608/1‑93 (MQ=255)

                               tCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTggg                                                             >  1:1213953/1‑93 (MQ=255)

                                ccTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGc                                                            <  2:4201978/93‑1 (MQ=255)

                                 cTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCa                                                           <  2:5646643/93‑1 (MQ=255)

                                  tGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCAc                                                          <  1:4106830/93‑1 (MQ=255)

                                   gacgacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACt                                                         <  1:187526/93‑1 (MQ=255)

                                     cgacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGc                                                       >  1:3848621/1‑93 (MQ=255)

                                      gacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCt                                                      >  2:1445565/1‑93 (MQ=255)

                                      gacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCt                                                      >  2:2306773/1‑93 (MQ=255)

                                       acAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTa                                                     <  2:4582025/93‑1 (MQ=255)

                                       acAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTa                                                     <  2:11913/93‑1 (MQ=255)

                                        cAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAt                                                    <  2:3628369/93‑1 (MQ=255)

                                         aaTCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAtt                                                   >  1:763473/1‑93 (MQ=255)

                                         aaTCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAtt                                                   <  1:678541/93‑1 (MQ=255)

                                          aTCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTc                                                  <  1:948465/93‑1 (MQ=255)

                                           tCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCt                                                 >  1:3573972/1‑93 (MQ=255)

                                            cTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCtt                                                <  1:2552039/93‑1 (MQ=255)

                                             tACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTa                                               >  2:3457012/1‑93 (MQ=255)

                                              aCGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTaa                                              >  1:2085060/1‑93 (MQ=255)

                                              aCGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTaa                                              <  2:7413477/93‑1 (MQ=255)

                                               cGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAAt                                             <  2:4915792/93‑1 (MQ=255)

                                                  tATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATccc                                          >  2:553949/1‑93 (MQ=255)

                                                  tATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATccc                                          <  1:513193/93‑1 (MQ=255)

                                                   aTTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCg                                         >  2:703808/1‑93 (MQ=255)

                                                    ttCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGt                                        <  2:5327537/93‑1 (MQ=255)

                                                       tGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTa                                     <  1:3521447/93‑1 (MQ=255)

                                                       tGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTa                                     >  1:5449609/1‑93 (MQ=255)

                                                       tGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTa                                     <  2:4671444/93‑1 (MQ=255)

                                                        nntgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTAt                                    <  1:2868125/91‑1 (MQ=255)

                                                        ggTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTAt                                    <  1:4051949/93‑1 (MQ=255)

                                                         gtgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATc                                   >  1:2206766/1‑93 (MQ=255)

                                                          tgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCg                                  <  2:2428136/93‑1 (MQ=255)

                                                          tgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCg                                  >  2:6539759/1‑93 (MQ=255)

                                                          tgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCg                                  <  2:3189925/93‑1 (MQ=255)

                                                            tttttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggc                                <  1:87135/93‑1 (MQ=255)

                                                                 aCCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggctggct                           >  2:2076220/1‑93 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         <  2:2483913/93‑1 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         <  2:3974584/93‑1 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         >  2:6236874/1‑93 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         <  1:5627263/93‑1 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         <  2:1519187/93‑1 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGGttt                         >  1:1387023/1‑92 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         <  1:3899503/93‑1 (MQ=255)

                                                                   cGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                         >  1:5827559/1‑93 (MQ=255)

                                                                    gtggtgGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttat                        <  2:5899539/93‑1 (MQ=255)

                                                                     tggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatt                       <  2:720265/93‑1 (MQ=255)

                                                                      ggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatta                      <  2:896783/93‑1 (MQ=255)

                                                                      ggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatta                      <  2:4068884/93‑1 (MQ=255)

                                                                       gtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTtattat                     <  1:4262374/93‑1 (MQ=255)

                                                                       gtggGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTtattat                     <  1:65798/93‑1 (MQ=255)

                                                                        tggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCGGGCTTATTatc                    <  2:4658466/93‑1 (MQ=255)

                                                                         gggCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatca                   <  1:3323838/93‑1 (MQ=255)

                                                                           gCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatcatc                 <  1:267775/93‑1 (MQ=255)

                                                                           gCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatcatc                 <  2:2343956/93‑1 (MQ=255)

                                                                             ggTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGt               >  2:4390398/1‑93 (MQ=255)

                                                                               ttCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCt             <  2:2396002/93‑1 (MQ=255)

                                                                                 cAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCttt           <  2:2557114/93‑1 (MQ=255)

                                                                                  aGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTa          <  2:5955343/93‑1 (MQ=255)

                                                                                  aGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTa          <  1:4154790/93‑1 (MQ=255)

                                                                                   gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGGTTATTTTTATCGTCTTTaa         >  1:2607435/1‑93 (MQ=255)

                                                                                   gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa         <  1:5610986/93‑1 (MQ=255)

                                                                                   gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa         <  1:4079455/93‑1 (MQ=255)

                                                                                   gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa         <  2:1279528/93‑1 (MQ=255)

                                                                                   gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa         <  1:5884251/93‑1 (MQ=255)

                                                                                        gccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAAATGcc    >  1:3505134/1‑93 (MQ=255)

                                                                                         ccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCa   <  1:7433691/93‑1 (MQ=255)

                                                                                          cgccTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAg  <  2:2875901/93‑1 (MQ=255)

TTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCAG  >  NC_000913/3250073‑3250255

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