BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,322,882	C→T	T146T (ACG→ACA)	phnB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,322,882	0	C	T	99.6%	862.4 / NA	235	T146T (ACG→ACA)	phnB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (123/111); minor base G (1/0); total (124/111)
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 257 total aligned reads displayed.

CGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTG  >  NC_000913/4322792‑4322971

 gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgtt                                                                                        >  2:1033250/1‑93 (MQ=255)

 gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgtt                                                                                        >  1:5672250/1‑93 (MQ=255)

  aTTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTgttg                                                                                       >  2:5143719/1‑93 (MQ=255)

   tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTTCCGGCGGGGTATTGTgtttg                                                                                      >  2:5327778/1‑91 (MQ=255)

    ttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttggt                                                                                     >  1:3739571/1‑93 (MQ=255)

     tGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttggtt                                                                                    >  1:3961125/1‑93 (MQ=255)

      ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTgttggttg                                                                                   >  1:2967761/1‑93 (MQ=255)

       gTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGgttgt                                                                                  >  1:609836/1‑93 (MQ=255)

        ttaattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGgttgtt                                                                                 <  2:765567/93‑1 (MQ=255)

         taattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGgttgttg                                                                                >  1:4029486/1‑93 (MQ=255)

         taattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTgttgggtggtt                                                                                >  1:3302371/1‑90 (MQ=255)

         taattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTttgttgttg                                                                                >  2:181697/1‑93 (MQ=255)

           attaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGgttgttgtt                                                                              >  1:5807611/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTga                                                                           <  1:5862696/93‑1 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTga                                                                           >  1:5177973/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTga                                                                           >  2:3793731/1‑93 (MQ=255)

              aaaaGTGACTAACAGATGAAGAGATAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCGCCCGGCGGGGTATTGTGGTGTCTGTTGttttg                                                                           >  2:4382779/1‑91 (MQ=255)

               aaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTgac                                                                          <  2:1416450/93‑1 (MQ=255)

               aaaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTgac                                                                          <  1:541580/93‑1 (MQ=255)

                 aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTgacga                                                                        <  1:3228738/93‑1 (MQ=255)

                    gACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACAt                                                                     <  2:4173284/93‑1 (MQ=255)

                       tAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTaa                                                                  >  2:3615894/1‑93 (MQ=255)

                        aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                 <  1:287189/93‑1 (MQ=255)

                        aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAat                                                                 >  1:4224860/1‑93 (MQ=255)

                          cAGATGAAGATTTAACGGGAATTCGATGACAGTGCGGGAGGGCGTGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatca                                                               <  1:3854654/93‑1 (MQ=255)

                          cAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatca                                                               >  1:5936807/1‑93 (MQ=255)

                            gATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAatcatc                                                             <  1:1019267/93‑1 (MQ=255)

                            gATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAatcatc                                                             >  2:1464420/1‑93 (MQ=255)

                            gATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACAATAatcatc                                                             >  1:5536113/1‑93 (MQ=255)

                             aTGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATcc                                                            >  1:1904287/1‑93 (MQ=255)

                             aTGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGGTGGTGGTTGACGACATTAATCATcc                                                            >  2:3016129/1‑93 (MQ=255)

                             aTGAAGAGTGAACGGGAATTCGATGACAGTGCGGGAGGGCGGTCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATcc                                                            <  2:4849428/93‑1 (MQ=255)

                              tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCa                                                           >  2:833737/1‑93 (MQ=255)

                              tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCa                                                           <  2:2112003/93‑1 (MQ=255)

                              tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCa                                                           <  2:425993/93‑1 (MQ=255)

                              tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCa                                                           <  1:4894748/93‑1 (MQ=255)

                               gAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGCCATTAATCATCCAc                                                          >  1:730571/1‑93 (MQ=255)

                                aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACg                                                         <  2:2460599/93‑1 (MQ=255)

                                 agagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACgg                                                        <  1:5375421/93‑1 (MQ=255)

                                  gagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGt                                                       >  2:3779840/1‑93 (MQ=255)

                                  gagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTTTTTTGTTTGTTGTTGTTTTTTGGCATTAATCATTCAccgg                                                       >  2:1097282/1‑90 (MQ=255)

                                    gTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTAc                                                     >  1:249783/1‑93 (MQ=255)

                                    gTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTAc                                                     <  2:192044/93‑1 (MQ=255)

                                      tAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGc                                                   >  2:571006/1‑93 (MQ=255)

                                       aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGcc                                                  <  2:2549967/93‑1 (MQ=255)

                                        aCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCa                                                 >  1:1543594/1‑93 (MQ=255)

                                         cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCaa                                                <  2:5654988/93‑1 (MQ=255)

                                         cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCaa                                                >  2:4632082/1‑93 (MQ=255)

                                          gggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTTATCATCCACGGTACGCCaaa                                               >  1:3368003/1‑93 (MQ=255)

                                           ggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAt                                              >  2:3199336/1‑93 (MQ=255)

                                           ggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAt                                              >  2:5039879/1‑93 (MQ=255)

                                             aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTTTTTGACGACATTAATCATCCACGGTACGCCAAAttt                                            <  2:2796739/93‑1 (MQ=255)

                                             aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAttt                                            >  2:4785218/1‑93 (MQ=255)

                                             aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAttt                                            >  1:355316/1‑93 (MQ=255)

                                               ttCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTAt                                          >  1:4557305/1‑93 (MQ=255)

                                               ttCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTAt                                          <  1:686270/93‑1 (MQ=255)

                                                  gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCgg                                       <  1:1808022/93‑1 (MQ=255)

                                                   aTGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGt                                      >  2:4001256/1‑93 (MQ=255)

                                                    tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTg                                     >  1:1848934/1‑93 (MQ=255)

                                                    tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTg                                     >  1:6164550/1‑93 (MQ=255)

                                                    tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTg                                     >  1:1715725/1‑93 (MQ=255)

                                                      aCAGTGCGGGAGGGCGGGCCTCCCGGCGGGGTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGAc                                   >  2:4500763/1‑93 (MQ=255)

                                                        aGTGCGGGAGGGCGGGCCTCCCGGCGGGATATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCCAATTTTTCGGTGAcgt                                 >  2:2582572/1‑91 (MQ=255)

                                                         gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACttt                                <  2:1744042/93‑1 (MQ=255)

                                                            cGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGcc                             <  1:4268654/93‑1 (MQ=255)

                                                            cGGGAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGcc                             <  1:4477526/93‑1 (MQ=255)

                                                              ggAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaa                           <  2:586939/93‑1 (MQ=255)

                                                               gAGGGCGGGCCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaaa                          <  2:5520949/93‑1 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  1:2398866/93‑1 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  1:1073316/93‑1 (MQ=255)

                                                                   gcgggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCa                      <  2:1575709/93‑1 (MQ=255)

                                                                     gggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATg                    >  1:6183179/1‑93 (MQ=255)

                                                                      ggcCTCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATgc                   <  1:3838985/93‑1 (MQ=255)

                                                                          tCCCGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc               >  1:2739816/1‑93 (MQ=255)

                                                                            ccGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAg             >  1:1260993/1‑93 (MQ=255)

                                                                            ccGGCGGGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAg             >  2:753763/1‑93 (MQ=255)

                                                                             cggcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGa            >  1:3473398/1‑93 (MQ=255)

                                                                               gcgtGTTATGGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa          <  1:5826625/93‑1 (MQ=255)

                                                                               gcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCCTGCGCCCAGaaa          >  1:4754047/1‑93 (MQ=255)

                                                                               gcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa          >  2:1280606/1‑93 (MQ=255)

                                                                               gcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa          <  1:823719/93‑1 (MQ=255)

                                                                               gcggGTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa          >  2:3517303/1‑93 (MQ=255)

                                                                                 gggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGGCTTTGCCAAAGCCATGCGCCCAGAAAGt        >  2:2194487/1‑93 (MQ=255)

                                                                                 gggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt        >  2:59177/1‑93 (MQ=255)

                                                                                 gggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt        >  1:2815272/1‑93 (MQ=255)

                                                                                 gggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt        >  2:5357409/1‑93 (MQ=255)

                                                                                 gggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt        >  2:2060072/1‑93 (MQ=255)

                                                                                  ggTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGtt       >  1:5237966/1‑93 (MQ=255)

                                                                                  ggTTAGTGTGTTGGTTGGCGTGCGTCGACATTAATCAGCCACGTTACGCCTAACTCATCGGTGACGTTGCCAAAGCCATGCGCCCAGAAAGtt       >  2:6296945/1‑93 (MQ=255)

                                                                                   gTTATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGttt      <  1:3020032/93‑1 (MQ=255)

                                                                                    ttATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTc     <  2:2892769/93‑1 (MQ=255)

                                                                                    ttATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTc     <  2:2422303/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  2:1982684/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    >  1:2663024/1‑93 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    <  2:147032/93‑1 (MQ=255)

                                                                                     tATTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTcc    >  2:3239780/1‑93 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTGATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt   >  1:2195992/1‑93 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt   >  1:795868/1‑93 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt   >  2:1189551/1‑93 (MQ=255)

                                                                                      aTTGTGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt   >  2:6058079/1‑93 (MQ=255)

CGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTG  >  NC_000913/4322792‑4322971

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