BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,980,297	(C)_8→9	noncoding (79/87 nt)	leuT →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,980,289	1	.	C	88.6%	502.9 / 62.5	185	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (3/18); new base C (58/106); total (61/124)
Fisher's exact test for biased strand distribution p-value = 8.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 228 total aligned reads displayed.

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAC  >  NC_000913/3980197‑3980376

cATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑t                                                                                            >  1:646867/1‑93 (MQ=255)

  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑cc                                                                                         <  1:1525762/93‑1 (MQ=255)

     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGCTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc                                                                                      >  2:5923187/1‑93 (MQ=255)

       aGAAGTTGTGACAATGCGAAGGGGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑T‑cccccac                                                                                    >  2:2683247/1‑91 (MQ=255)

        gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc‑ccc                                                                                  >  1:4461496/1‑93 (MQ=255)

         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc‑cccc                                                                                 <  1:4242489/93‑2 (MQ=255)

         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc‑cccc                                                                                 <  1:3823666/93‑2 (MQ=255)

         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc‑cccc                                                                                 <  2:1657058/93‑2 (MQ=255)

          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑ccccc‑cccct                                                                                <  1:3970932/93‑3 (MQ=255)

          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑G‑ccccc‑ccccc                                                                                >  2:19388/1‑91 (MQ=255)

               tgACAATGCGAAGGTGGCGGAATTGGTAGACGCCCTAGCTTCAGGTGTTAGTTTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGcac                                                                            <  1:2990981/93‑1 (MQ=11)

                 aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGcacca                                                                          <  2:5663585/93‑1 (MQ=16)

                 aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGcacca                                                                          <  2:2539799/93‑1 (MQ=16)

                 aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGcacca                                                                          <  1:5885135/93‑1 (MQ=16)

                  cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGcaccac                                                                         >  1:3223518/1‑93 (MQ=16)

                   aaTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACg                                                                        <  1:3188383/93‑1 (MQ=16)

                    aTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACgc                                                                       >  1:4179076/1‑92 (MQ=12)

                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑c                                                                     >  2:6304719/1‑93 (MQ=255)

                     tGCGAAGGTGGCAGAATTGGCAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGC‑C‑CCCCC‑CCCTCGCACCACGA‑c                                                                     <  2:2610350/93‑1 (MQ=11)

                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑Ct                                                                    <  2:5548762/93‑1 (MQ=255)

                       cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑Ctt                                                                   <  1:5155682/93‑1 (MQ=255)

                       cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑CCCCC‑CCCTCGCACCACGA‑Cttt                                                                  <  1:4129609/93‑1 (MQ=38)

                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑Cttt                                                                  <  2:17314/93‑1 (MQ=255)

                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTa                                                                 <  1:5996325/93‑1 (MQ=255)

                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTa                                                                 >  1:5551416/1‑93 (MQ=255)

                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTa                                                                 >  2:3942845/1‑93 (MQ=255)

                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTa                                                                 >  2:5306171/1‑93 (MQ=255)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGTTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGC‑C‑CCCCC‑CCCTCGCACCACGA‑CTTTaa                                                                <  2:672972/93‑1 (MQ=14)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaa                                                                >  2:2184875/1‑93 (MQ=255)

                           gttGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaa                                                               <  1:1134540/91‑1 (MQ=255)

                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaa                                                               <  1:2598622/93‑1 (MQ=255)

                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaa                                                               >  1:1600716/1‑93 (MQ=255)

                           ggTGACGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaa                                                               <  1:1992455/93‑1 (MQ=255)

                            gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaaa                                                              >  1:6095822/1‑92 (MQ=255)

                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTaaaaaa                                                            >  1:6391715/1‑93 (MQ=255)

                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGaa                                                            <  2:5727435/93‑1 (MQ=255)

                               gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAAt                                                           <  2:841340/93‑1 (MQ=255)

                                cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAAtt                                                          <  2:5426945/93‑1 (MQ=255)

                                cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAAtt                                                          <  1:6222873/93‑1 (MQ=255)

                                 ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTg                                                         <  1:3340932/93‑1 (MQ=255)

                                 ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTg                                                         >  1:1929225/1‑93 (MQ=255)

                                  gAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGa                                                        <  2:4929338/93‑1 (MQ=255)

                                   aaTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGaa                                                       <  2:516237/93‑1 (MQ=255)

                                    aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAAc                                                      <  2:6481775/93‑1 (MQ=255)

                                     ttGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCCCCG‑CCTTAAAGAAATGAAca                                                     >  2:5880672/1‑92 (MQ=11)

                                     ttGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACt                                                     >  2:3845791/1‑93 (MQ=255)

                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGCCC‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACt                                                     <  1:6376883/93‑1 (MQ=17)

                                        gTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaa                                                  <  1:4982203/93‑1 (MQ=255)

                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaaa                                                 <  2:1045426/93‑1 (MQ=255)

                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaaaa                                                <  2:2839402/93‑1 (MQ=255)

                                         tAGACGCGATAGCTTCATGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaaa                                                 <  1:2457836/93‑1 (MQ=18)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGC‑C‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaaaa                                                <  1:4900624/93‑1 (MQ=255)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTaaaaa                                                >  2:784520/1‑93 (MQ=255)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCAACACGA‑CTTTAAAGAATTGAACTaaaaa                                                >  2:1887769/1‑93 (MQ=255)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTGCGGACGTGGGGGTTCAAGC‑C‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAAt                                               <  1:1461779/93‑1 (MQ=18)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAAt                                               <  2:5258655/93‑1 (MQ=255)

                                             cgcgCTAGCTTCAGTTGTTAGTGTCCTTACCGACGTGGTGTTTCAAG‑‑TCCCCCC‑CCCTCGCCCCACGA‑CTTTAAAGAATTGAACTAAAAATTc                                             <  1:666256/93‑1 (MQ=21)

                                              gcgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCa                                            >  2:1640255/1‑93 (MQ=255)

                                              gcgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCa                                            <  1:778915/93‑1 (MQ=255)

                                              gcgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑TACCCCC‑CCCTCACACCACGA‑GATTAAAAAAATGAAATaaaaaatta                                            >  1:3715691/1‑89 (MQ=14)

                                               cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCaa                                           >  2:6419458/1‑93 (MQ=255)

                                               cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCaa                                           <  2:5241295/93‑1 (MQ=255)

                                               cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑TCCCCCC‑CCCTCCCCCCACCACCTTTAAAAAATTAAACAAAAAAttaa                                            >  1:2496824/1‑91 (MQ=9)

                                                gcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCaaa                                          >  2:3107222/1‑93 (MQ=255)

                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCCCACCACCA‑CTTTAAAAAATTTAAATAAAAAATCaaaaa                                        >  2:4717973/1‑93 (MQ=17)

                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                      >  1:868534/1‑93 (MQ=255)

                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  2:3373495/93‑1 (MQ=255)

                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  2:6077932/93‑1 (MQ=255)

                                                      ttgaGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  1:4737460/90‑1 (MQ=255)

                                                      ttCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    >  1:5204209/1‑93 (MQ=255)

                                                      ttCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  1:2233818/93‑1 (MQ=255)

                                                      ttCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  1:3580312/93‑1 (MQ=255)

                                                      ttCAGGTGTGAGTGTCCTGACGGACTTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  2:1777749/93‑1 (MQ=25)

                                                       tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTTGCCCCACGA‑ATTTAAAAAAATGCACTAAAAATTTAAAAAACAga                                   >  1:2750127/1‑92 (MQ=17)

                                                       tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGt                                   <  2:2362911/93‑1 (MQ=255)

                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  1:2700027/93‑1 (MQ=255)

                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  1:1025198/93‑1 (MQ=255)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTCAAAATTCAACAAGCAGTAt                                 >  2:4241510/1‑93 (MQ=21)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAt                                 <  2:4463193/93‑1 (MQ=255)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAt                                 <  1:4867499/93‑1 (MQ=255)

                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  1:1635129/93‑1 (MQ=255)

                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  2:6543815/93‑1 (MQ=255)

                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  1:204024/93‑1 (MQ=255)

                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  1:2712904/93‑1 (MQ=255)

                                                            ggtTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  2:218973/92‑1 (MQ=255)

                                                             gtTACTGTGCTCCCGGCCTTGGGGGTGCAAG‑‑T‑CCCCCGCCCTCGCACCACTT‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                             <  1:1077348/93‑1 (MQ=255)

                                                              ttAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                            <  2:1129588/93‑1 (MQ=255)

                                                              ttAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                            <  1:3405442/93‑1 (MQ=255)

                                                               tAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGc                           <  2:436094/93‑1 (MQ=255)

                                                                   gtCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGt                       <  2:5176898/93‑1 (MQ=255)

                                                                     ccTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                     <  1:3385624/93‑1 (MQ=255)

                                                                      cTTACGGACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgc                    <  2:3712590/93‑1 (MQ=255)

                                                                           ggACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc               <  1:701553/93‑1 (MQ=255)

                                                                            gACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              >  2:1514310/1‑93 (MQ=255)

                                                                            gACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              <  1:6731963/93‑1 (MQ=255)

                                                                            gACGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              <  1:3051053/93‑1 (MQ=255)

                                                                             aCGTGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGGAGCGCAGCtt             >  1:5421793/1‑93 (MQ=255)

                                                                                tGGGGGTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGGATTTCGGCGAGGAGCGCAGCTTggg          >  1:1376957/1‑92 (MQ=39)

                                                                                  ggggTTCAAG‑‑TCCCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg        >  2:4170740/1‑93 (MQ=255)

                                                                                       tCAAG‑‑T‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc  <  2:3536298/93‑1 (MQ=255)

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑CCCCC‑CCCTCGCACCACGA‑CTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAC  >  NC_000913/3980197‑3980376

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