BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,954,264	A→T	D881E (GAT→GAA)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,954,264	0	A	T	98.9%	620.5 / ‑0.7	178	D881E (GAT→GAA)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (79/97); total (80/98)
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.50e-01

Only 100 of 187 total aligned reads displayed.

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGC  >  NC_000913/2954174‑2954351

cTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTc                                                                                       >  1:2487939/1‑93 (MQ=255)

cTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGCCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTc                                                                                       <  2:1758301/93‑1 (MQ=255)

 tGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGTGCCACGATTTTATCACGAGAATGTAATCGCATATTGCCGCGTTCg                                                                                      >  1:3817454/1‑93 (MQ=255)

 tGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCg                                                                                      >  1:3686468/1‑93 (MQ=255)

 tGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCg                                                                                      >  2:1957579/1‑93 (MQ=255)

 tGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCg                                                                                      >  2:4349305/1‑93 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                     <  1:425606/93‑1 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                     >  1:3258578/1‑93 (MQ=255)

  gATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGa                                                                                     >  2:4088047/1‑93 (MQ=255)

    tGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGaaa                                                                                   <  1:1536618/93‑1 (MQ=255)

     ggAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAAt                                                                                  >  2:3650288/1‑93 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                 >  1:1133821/1‑93 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                 <  1:3262868/93‑1 (MQ=255)

      gaagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATc                                                                                 <  1:1273586/93‑1 (MQ=255)

       aagaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCt                                                                                >  1:3542236/1‑93 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                               <  1:2394338/93‑1 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                               <  2:1387889/93‑1 (MQ=255)

        agaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCtt                                                                               <  2:3635235/93‑1 (MQ=255)

         gaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCttt                                                                              >  2:3207139/1‑93 (MQ=255)

         gaaATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCttt                                                                              <  2:1597067/93‑1 (MQ=255)

          aaaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCGACGATTTTAACACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTa                                                                             <  2:1015115/93‑1 (MQ=255)

          aaaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTa                                                                             <  1:3989222/93‑1 (MQ=255)

          aaaTCAGCAAGTTGTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTa                                                                             <  1:597146/93‑1 (MQ=255)

           aaTCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTAc                                                                            >  1:1744796/1‑93 (MQ=255)

             tCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACtt                                                                          >  2:864150/1‑93 (MQ=255)

             tCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACtt                                                                          >  2:1280774/1‑93 (MQ=255)

              cAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTGCTTa                                                                         <  2:2535785/93‑1 (MQ=255)

              cAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTa                                                                         >  1:2759142/1‑93 (MQ=255)

              cAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTa                                                                         <  1:3285877/93‑1 (MQ=255)

               aGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTac                                                                        >  1:4082915/1‑92 (MQ=255)

               aGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTaa                                                                        >  1:1036532/1‑93 (MQ=255)

                 cAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACt                                                                      >  2:1433136/1‑93 (MQ=255)

                 cAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACt                                                                      >  1:3705186/1‑93 (MQ=255)

                 cAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACt                                                                      <  2:913091/93‑1 (MQ=255)

                 cAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACt                                                                      <  2:1504664/93‑1 (MQ=255)

                  aaGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACtt                                                                     >  2:3523585/1‑93 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGGCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                    <  2:4023438/93‑1 (MQ=255)

                   aGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTc                                                                    <  1:764635/93‑1 (MQ=255)

                    gTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCg                                                                   <  2:724899/93‑1 (MQ=255)

                     tttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGa                                                                  >  2:976076/1‑93 (MQ=255)

                      ttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGAt                                                                 <  2:2606368/93‑1 (MQ=255)

                      ttttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGAt                                                                 <  1:3470615/93‑1 (MQ=255)

                       tttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATg                                                                <  2:3984699/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                               <  2:477402/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                               <  1:410322/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                               <  2:389642/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                               <  2:2414043/93‑1 (MQ=255)

                        ttGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGc                                                               <  2:213552/93‑1 (MQ=255)

                         tGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGct                                                              <  2:3382895/93‑1 (MQ=255)

                          gcggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGctt                                                             >  1:807046/1‑93 (MQ=255)

                            ggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                           >  1:1124389/1‑93 (MQ=255)

                            ggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                           <  1:1387107/93‑1 (MQ=255)

                            ggcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttct                                                           <  2:4133339/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                          >  1:4441515/1‑93 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                          <  2:810361/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGcttctt                                                          <  2:2925864/93‑1 (MQ=255)

                             gcgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGcttctt                                                          <  1:3783422/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                         <  2:428268/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                         <  1:4109245/93‑1 (MQ=255)

                              cgTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCttcttc                                                         <  2:2644571/93‑1 (MQ=255)

                               gTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGGAATCTTTGCTTAACTTCGATGCTTCTCCg                                                        >  2:2246245/1‑93 (MQ=255)

                                tCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGc                                                       <  2:623148/93‑1 (MQ=255)

                                     agTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGc                                                  >  1:2925320/1‑93 (MQ=255)

                                      gTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCg                                                 >  1:1869451/1‑93 (MQ=255)

                                       ttttATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGt                                                >  2:3220876/1‑93 (MQ=255)

                                        tttATCTTGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGtt                                               <  2:3316607/93‑1 (MQ=255)

                                          tATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTg                                             <  1:1005616/93‑1 (MQ=255)

                                           ttCTGGGGCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGc                                            <  1:1346631/92‑1 (MQ=255)

                                           aTCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGc                                            <  2:2847852/93‑1 (MQ=255)

                                           aTCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGc                                            <  2:2994811/93‑1 (MQ=255)

                                           aTCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGc                                            <  1:265450/93‑1 (MQ=255)

                                            tCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCg                                           <  2:761288/93‑1 (MQ=255)

                                               gggCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTg                                        >  2:583673/1‑93 (MQ=255)

                                                 gCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGcc                                      <  1:2016345/93‑1 (MQ=255)

                                                 gCCACGATTTTATCACGCGAAGCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGcc                                      <  1:3721054/93‑1 (MQ=255)

                                                  ccACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCt                                     <  1:1574791/93‑1 (MQ=255)

                                                   cACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTg                                    <  1:2799998/93‑1 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                   >  1:4715514/1‑93 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                   <  2:665372/93‑1 (MQ=255)

                                                    aCGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgc                                   >  1:3086654/1‑93 (MQ=255)

                                                      gATTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcag                                 >  1:1956512/1‑93 (MQ=255)

                                                       aTTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcagc                                >  1:3043141/1‑93 (MQ=255)

                                                       aTTTTATCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTgcagc                                >  2:3375677/1‑93 (MQ=255)

                                                          ttATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATc                             >  2:3296838/1‑93 (MQ=255)

                                                             tCACGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTgg                          >  2:1413475/1‑93 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                        <  1:4502181/93‑1 (MQ=255)

                                                               aCGCGAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGt                        >  2:4506481/1‑93 (MQ=255)

                                                                   gAATCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAAtt                    >  1:4061459/1‑93 (MQ=255)

                                                                      tCGAAGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTAcc                 >  1:4492049/1‑93 (MQ=255)

                                                                         aaGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGcc              <  1:4414745/93‑1 (MQ=255)

                                                                         aaGCGCATATTGCCGCGTTAGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGcc              <  1:2220941/93‑1 (MQ=255)

                                                                          aGCTCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCct             <  2:2073403/93‑1 (MQ=255)

                                                                          aGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCct             <  2:469170/93‑1 (MQ=255)

                                                                          aGCGCATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCct             <  2:2073402/93‑1 (MQ=255)

                                                                           gcgcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctg            >  2:1606905/1‑93 (MQ=255)

                                                                            cgcATATTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCTTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCctgc           <  1:2079062/93‑1 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa       >  2:1159833/1‑93 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa       <  2:206006/93‑1 (MQ=255)

                                                                                tatTGCCGCGTTCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGa       >  1:2842586/1‑93 (MQ=255)

                                                                                   tGCCGCGTTCGAAGTCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGAttt    <  2:1208557/93‑1 (MQ=255)

CTGATGGAAGAAATCAGCAAGTTTTTGCGGCGTCAGCAGTTTTATCTGGGCCACGATTTTATCACGCGAATCGAAGCGCATATTGCCGCGATCGAAATCTTTACTTAACTTCGATGCTTCTTCGCCGAGCGTTTGCGGTGCCTGCAGCATCTGGGTAATTACCGCCTGCTGGATTTGC  >  NC_000913/2954174‑2954351

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