BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	667,973	C→T	R96H (CGC→CAC)	ybeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	667,973	0	C	T	98.0%	503.7 / ‑4.5	146	R96H (CGC→CAC)	ybeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/2); new base T (72/71); total (73/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCGAGC  >  NC_000913/667884‑668062

cGGTTTGTACCCAGTCCGGCATTTTCGTTCCCGCGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGAc                                                                                        <  1:2178747/93‑1 (MQ=255)

cGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGCCAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGag                                                                                        >  2:3306655/1‑92 (MQ=255)

 ggTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACg                                                                                       >  2:1374654/1‑93 (MQ=255)

  gTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACgg                                                                                      <  2:2025154/93‑1 (MQ=255)

   tttGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGc                                                                                     <  1:3577288/93‑1 (MQ=255)

   tttGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGc                                                                                     >  2:3192115/1‑93 (MQ=255)

    ttGTACCCAGTCCGGCATTTTCGTTCCCACGGTGACAAGTTGAAGCTTCACGCTTTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGct                                                                                    >  1:3997241/1‑93 (MQ=255)

      gTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGctct                                                                                  <  1:2502838/93‑1 (MQ=255)

       tACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGctctc                                                                                 <  2:4673576/93‑1 (MQ=255)

        aCCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGctctct                                                                                >  2:788121/1‑93 (MQ=255)

          ccAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCCTACAGGTGACGGCTCTCTTc                                                                              >  2:1012362/1‑93 (MQ=255)

          ccAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTc                                                                              <  1:2177106/93‑1 (MQ=255)

               ccGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCa                                                                         >  2:2577226/1‑93 (MQ=255)

                cGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCAt                                                                        <  1:4022481/93‑1 (MQ=255)

                 ggCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCAtc                                                                       >  2:1089050/1‑92 (MQ=255)

                 ggCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATg                                                                       <  2:2494655/93‑1 (MQ=255)

                 ggCATTTTCGTTCCCACGGCGACAAGTTGCAGCCTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATg                                                                       >  2:831933/1‑93 (MQ=255)

                  gCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGa                                                                      >  1:2798737/1‑93 (MQ=255)

                    aTTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACa                                                                    <  2:2914625/93‑1 (MQ=255)

                     ttttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACAt                                                                   <  2:3710502/93‑1 (MQ=255)

                     ttttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACAt                                                                   >  2:703165/1‑93 (MQ=255)

                     ttttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACAt                                                                   <  2:331605/93‑1 (MQ=255)

                     ttttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACAt                                                                   >  1:3259932/1‑93 (MQ=255)

                      tttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATg                                                                  <  1:1146257/93‑1 (MQ=255)

                       ttCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATgg                                                                 <  2:3552728/93‑1 (MQ=255)

                        tCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGa                                                                >  2:1495137/1‑93 (MQ=255)

                         cGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGAc                                                               >  2:3120641/1‑93 (MQ=255)

                         cGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGAc                                                               >  1:4185467/1‑93 (MQ=255)

                         cGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGAc                                                               >  1:3467435/1‑93 (MQ=255)

                           ttCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACaa                                                             >  1:2951674/1‑93 (MQ=255)

                           ttCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACaa                                                             >  1:4112322/1‑93 (MQ=255)

                            tCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAAt                                                            >  2:2216629/1‑93 (MQ=255)

                              ccACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATca                                                          >  1:4307168/1‑93 (MQ=255)

                                aCGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATcaca                                                        <  2:1409540/93‑1 (MQ=255)

                                 cGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACAt                                                       <  1:1694010/93‑1 (MQ=255)

                                  ggCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATc                                                      <  1:1462299/93‑1 (MQ=255)

                                  ggCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATc                                                      <  1:4415424/93‑1 (MQ=255)

                                     gACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGcc                                                   <  1:220646/93‑1 (MQ=255)

                                     gACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGcc                                                   <  2:1113370/93‑1 (MQ=255)

                                       cAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCa                                                 <  1:3431586/93‑1 (MQ=255)

                                       cAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCa                                                 >  1:1790948/1‑93 (MQ=255)

                                        aaGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCaa                                                >  1:1109582/1‑93 (MQ=255)

                                         aGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCaaa                                               <  1:3636949/93‑1 (MQ=255)

                                           ttGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATc                                             <  2:2966580/93‑1 (MQ=255)

                                           ttGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATc                                             <  2:4085222/93‑1 (MQ=255)

                                             gCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCa                                           >  2:3996456/1‑93 (MQ=255)

                                             gCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCa                                           >  2:4722424/1‑93 (MQ=255)

                                              cAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCAc                                          <  2:3672833/93‑1 (MQ=255)

                                               aGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATGCAGGTGACGGCTCTCTTCCTGCATGACATGGACAAGCACATCGCCCAAATCCACg                                         <  1:2357696/93‑1 (MQ=255)

                                               aGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACg                                         >  2:3078478/1‑93 (MQ=255)

                                               aGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCACAAATCCACg                                         >  2:3078477/1‑93 (MQ=255)

                                                gCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGa                                        >  2:2091792/1‑93 (MQ=255)

                                                gCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGa                                        <  2:1683737/93‑1 (MQ=255)

                                                gCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGACCAAATCCACGa                                        >  1:2722036/1‑93 (MQ=255)

                                                 cTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGAc                                       <  2:4386666/93‑1 (MQ=255)

                                                  ttCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACa                                      >  1:2895098/1‑93 (MQ=255)

                                                   tCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACaa                                     <  2:2269606/93‑1 (MQ=255)

                                                    cACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAAt                                    <  1:1240704/93‑1 (MQ=255)

                                                      cGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATcc                                  >  2:1410619/1‑93 (MQ=255)

                                                      cGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATcc                                  <  1:4409940/93‑1 (MQ=255)

                                                        cATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAg                                >  1:426813/1‑93 (MQ=255)

                                                        cATTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAg                                <  2:2981556/93‑1 (MQ=255)

                                                        cATTAACTCCAGAGTTTATCCAGTTCATGCAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAg                                <  1:4219951/93‑1 (MQ=255)

                                                         aTTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAgc                               >  2:2348067/1‑92 (MQ=255)

                                                         aTTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGt                               >  2:1308614/1‑93 (MQ=255)

                                                         aTTAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGt                               >  2:2935528/1‑93 (MQ=255)

                                                          ttAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATTCAGTc                              >  2:3616454/1‑93 (MQ=255)

                                                          ttAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTc                              <  1:2489522/93‑1 (MQ=255)

                                                           tAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCg                             <  2:4313041/93‑1 (MQ=255)

                                                           tAACTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCg                             >  1:2098235/1‑93 (MQ=255)

                                                            aaCTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCgg                            >  2:830509/1‑93 (MQ=255)

                                                            aaCTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCgg                            <  1:3652652/93‑1 (MQ=255)

                                                            aaCTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATTTCATCGACCAAATCCCCGACAATCCAGtgtg                            >  2:3896719/1‑90 (MQ=255)

                                                             aCTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGCCAATCACATCGCCCAAATCCACGACAATCCAGTCGgc                           <  1:1733449/93‑1 (MQ=255)

                                                              cTCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGgcg                          >  1:3045053/1‑93 (MQ=255)

                                                               tCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGgcgc                         <  1:1535800/93‑1 (MQ=255)

                                                               tCCAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGgcgc                         <  1:3846118/93‑1 (MQ=255)

                                                                ccAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCt                        >  1:2299902/1‑93 (MQ=255)

                                                                ccAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCt                        >  2:410964/1‑93 (MQ=255)

                                                                 cAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTg                       >  1:1572972/1‑93 (MQ=255)

                                                                 cAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTg                       >  2:4636767/1‑93 (MQ=255)

                                                                 cAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTg                       >  1:2292295/1‑93 (MQ=255)

                                                                 cAGAGTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGACCAAATCCACGACAATCCAGTCGGCGCTg                       >  1:2292296/1‑93 (MQ=255)

                                                                  agagTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGt                      <  2:4670597/93‑1 (MQ=255)

                                                                   gagTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtt                     >  2:497038/1‑93 (MQ=255)

                                                                    agTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGttt                    <  2:2328588/93‑1 (MQ=255)

                                                                    agTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGttt                    <  1:3925189/93‑1 (MQ=255)

                                                                     gTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtttt                   >  1:4309193/1‑93 (MQ=255)

                                                                     gTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtttt                   >  2:742782/1‑93 (MQ=255)

                                                                     gTTTTTCCAGTTCATACAGGTGACGGCTCTCTTCATGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtttt                   <  1:4583056/93‑1 (MQ=255)

                                                                       ttttCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCa                 <  1:4629586/93‑1 (MQ=255)

                                                                          tCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCt              <  2:1202994/93‑1 (MQ=255)

                                                                          tCCAGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCt              <  2:1202995/93‑1 (MQ=255)

                                                                             aGTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCt           <  1:2776557/93‑1 (MQ=255)

                                                                              gTTCATACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTa          <  2:4592453/93‑1 (MQ=255)

                                                                              gTTCAGACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTa          <  1:2682483/93‑1 (MQ=255)

                                                                                  aTACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGcc      <  1:2910378/93‑1 (MQ=255)

                                                                                  aTACAGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGcc      >  1:701183/1‑93 (MQ=255)

                                                                                      aGGTGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCGAGc  <  1:639256/93‑1 (MQ=255)

                                                                                      aGGTGACGGCTCCCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCGAGc  <  1:1520989/93‑1 (MQ=255)

CGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCGAGC  >  NC_000913/667884‑668062

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