BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	358,162	A→G	E47G (GAG→GGG)	cynT →	carbonic anhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	358,162	0	A	G	84.8%	747.8 / 128.4	328	E47G (GAG→GGG)	cynT	carbonic anhydrase
Reads supporting (aligned to +/- strand): ref base A (27/23); new base G (133/145); total (160/168)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 329 total aligned reads displayed.

CGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTT  >  NC_000913/358072‑358254

cGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGc                                                                                            <  1:6398897/93‑1 (MQ=255)

  aaGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTg                                                                                          <  1:2737199/93‑1 (MQ=255)

   aGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTgg                                                                                         <  1:1368307/93‑1 (MQ=255)

      gggAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGa                                                                                      <  1:674837/93‑1 (MQ=255)

       ggAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGAc                                                                                     >  1:2798719/1‑93 (MQ=255)

           gCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCaa                                                                                 <  2:576611/93‑1 (MQ=255)

           gCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCaa                                                                                 <  2:6697041/93‑1 (MQ=255)

            ccTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAAc                                                                                >  1:6120459/1‑93 (MQ=255)

              ttGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGt                                                                              >  1:7212443/1‑93 (MQ=255)

              ttGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGt                                                                              <  1:7560/93‑1 (MQ=255)

              ntGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACNt                                                                              >  2:5370000/2‑93 (MQ=255)

               tGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTTGTCCCTGGGCTGGTGACGCAACGTg                                                                             <  2:4358847/93‑1 (MQ=255)

               tGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTg                                                                             >  1:6598894/1‑93 (MQ=255)

                  ttAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGc                                                                          <  1:6843960/93‑1 (MQ=255)

                    aaaCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCt                                                                        <  1:7339090/93‑1 (MQ=255)

                     aaCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTg                                                                       <  1:5585604/93‑1 (MQ=255)

                      aCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTgg                                                                      <  1:2458778/93‑1 (MQ=255)

                      aCAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTgg                                                                      >  1:7045590/1‑93 (MQ=255)

                       cAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGc                                                                     >  2:7489548/1‑93 (MQ=255)

                        aGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCg                                                                    >  1:5355459/1‑93 (MQ=255)

                         gCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGa                                                                   <  2:5611147/93‑1 (MQ=255)

                           tGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATc                                                                 <  1:2666162/93‑1 (MQ=255)

                           tGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATc                                                                 >  1:5470190/1‑93 (MQ=255)

                           tGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCCGGCGATc                                                                 <  1:898768/93‑1 (MQ=255)

                            ggCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCt                                                                <  1:1802359/93‑1 (MQ=255)

                               gACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGtt                                                             <  2:6509114/93‑1 (MQ=255)

                               gACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGtt                                                             <  1:6500462/93‑1 (MQ=255)

                                acacAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTc                                                            >  2:4256757/1‑93 (MQ=255)

                                acacAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGCCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTc                                                            >  1:3963834/1‑93 (MQ=255)

                                  acaGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGt                                                          >  1:3942279/1‑93 (MQ=255)

                                    aGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTa                                                        <  2:7452453/93‑1 (MQ=255)

                                     gCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTAt                                                       <  2:488924/93‑1 (MQ=255)

                                     gCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTAt                                                       >  2:1883080/1‑93 (MQ=255)

                                     gCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTAt                                                       >  2:4204142/1‑93 (MQ=255)

                                     gCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTAt                                                       <  2:6178145/93‑1 (MQ=255)

                                       aaaGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTc                                                     <  1:3680697/93‑1 (MQ=255)

                                        aaGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCg                                                    <  2:1827098/93‑1 (MQ=255)

                                          gCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGAGGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCa                                                  >  1:1045115/1‑93 (MQ=255)

                                          gCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCa                                                  >  1:2048126/1‑93 (MQ=255)

                                           cccGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCaa                                                 >  1:1880432/1‑93 (MQ=255)

                                           cccGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCaa                                                 <  2:2119533/93‑1 (MQ=255)

                                            ccGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAAc                                                >  2:1600693/1‑93 (MQ=255)

                                              gcgcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAAcgc                                              >  2:6004339/1‑93 (MQ=255)

                                              gcgcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAAcgc                                              >  2:3474638/1‑93 (MQ=255)

                                                gcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCgg                                            <  2:4465466/93‑1 (MQ=255)

                                                gcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCgg                                            <  1:4567226/93‑1 (MQ=255)

                                                gcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCgg                                            >  2:3812697/1‑93 (MQ=255)

                                                gcACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCgg                                            <  1:6992373/93‑1 (MQ=255)

                                                 cacaCTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCggg                                           >  1:383932/1‑93 (MQ=255)

                                                 cacaCTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCggg                                           >  2:2657347/1‑93 (MQ=255)

                                                  acacTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGc                                          >  1:4820882/1‑93 (MQ=255)

                                                    acTTTTTATCCCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCaa                                        <  1:623237/93‑1 (MQ=255)

                                                       ttttATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAatat                                     <  2:1026445/93‑1 (MQ=255)

                                                       ttttATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAatat                                     >  1:7071574/1‑93 (MQ=255)

                                                       ttttATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAatat                                     <  1:6765459/93‑1 (MQ=255)

                                                        tttATCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATc                                    >  1:1502025/1‑93 (MQ=255)

                                                        tttATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATc                                    <  1:5300399/93‑1 (MQ=255)

                                                          tATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGt                                  >  2:1149535/1‑93 (MQ=255)

                                                           aTCTCCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTc                                 <  2:3326672/93‑1 (MQ=255)

                                                            tctcCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTcc                                >  1:2005260/1‑93 (MQ=255)

                                                              tcCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCt                              >  1:5875022/1‑93 (MQ=255)

                                                              tcCTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCt                              >  1:7423881/1‑93 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            <  2:5561824/93‑1 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            <  1:1648164/93‑1 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            >  2:680464/1‑93 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            >  2:3764030/1‑93 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            <  2:3087805/93‑1 (MQ=255)

                                                                cTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTc                            >  2:5092013/1‑93 (MQ=255)

                                                                  gCTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCt                          >  2:3682098/1‑93 (MQ=255)

                                                                   cTCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTa                         <  2:3556915/93‑1 (MQ=255)

                                                                    tCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTAc                        >  2:2454278/1‑93 (MQ=255)

                                                                    tCCGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTAc                        >  1:1979901/1‑93 (MQ=255)

                                                                    tCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTAc                        >  2:2974454/1‑93 (MQ=255)

                                                                     ccGACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACg                       >  1:5029484/1‑93 (MQ=255)

                                                                      cGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACgg                      >  2:3966733/1‑93 (MQ=255)

                                                                       gACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACggg                     >  1:5155027/1‑93 (MQ=255)

                                                                       gACAGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACggg                     <  1:485640/93‑1 (MQ=255)

                                                                       gACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACggg                     <  2:1392801/93‑1 (MQ=255)

                                                                          aGCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCg                  <  1:4228607/93‑1 (MQ=255)

                                                                           gCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCgg                 <  1:3053363/93‑1 (MQ=255)

                                                                           gCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCgg                 >  2:4623347/1‑93 (MQ=255)

                                                                           gCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCgg                 >  1:2218744/1‑93 (MQ=255)

                                                                           gCCGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCgg                 >  1:3353438/1‑93 (MQ=255)

                                                                             cGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGaa               <  2:5896057/93‑1 (MQ=255)

                                                                             cGTCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGaa               <  1:3199303/93‑1 (MQ=255)

                                                                             cGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGaa               <  1:2836394/93‑1 (MQ=255)

                                                                             cGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGaa               <  1:6062081/93‑1 (MQ=255)

                                                                               tCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAAcc             >  2:7289784/1‑93 (MQ=255)

                                                                               tCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAAcc             >  2:6320223/1‑93 (MQ=255)

                                                                               tCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAAcc             >  2:856239/1‑93 (MQ=255)

                                                                               tCTGGTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAAcc             >  1:5969206/1‑93 (MQ=255)

                                                                                 tGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCg           >  1:4064076/1‑93 (MQ=255)

                                                                                  ggTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCgg          <  2:7192774/93‑1 (MQ=255)

                                                                                   gTCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGt         <  1:1717995/93‑1 (MQ=255)

                                                                                    tCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTg        <  2:2232588/93‑1 (MQ=255)

                                                                                    tCCCTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTg        <  1:248765/93‑1 (MQ=255)

                                                                                      ccTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGc      <  2:5214949/93‑1 (MQ=255)

                                                                                      ccTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGc      >  1:5763701/1‑93 (MQ=255)

                                                                                      ccTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGc      <  1:4471889/93‑1 (MQ=255)

                                                                                      ccTGGGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGc      <  1:1273857/93‑1 (MQ=255)

CGAAGCGGGAAGCCTTGTTTAAACAGCTGGCGACACAGCAAAGCCCGCGCACACTTTTTATCTCCTGCTCCGACAGCCGTCTGGTCCCTGAGCTGGTGACGCAACGTGAGCCTGGCGATCTGTTCGTTATTCGCAACGCGGGCAATATCGTCCCTTCCTACGGGCCGGAACCCGGTGGCGTTT  >  NC_000913/358072‑358254

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