BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,173,376	A→G	noncoding (44/76 nt)	thrT →	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,173,376	0	A	G	99.7%	1068.6 / NA	345	noncoding (44/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (160/184); total (160/185)
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 355 total aligned reads displayed.

CCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAT  >  NC_000913/4173286‑4173465

 cAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGt                                                                                        <  1:1132428/93‑1 (MQ=255)

 cAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGt                                                                                        >  1:3739376/1‑93 (MQ=255)

 cAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGt                                                                                        >  1:1839913/1‑93 (MQ=255)

    gCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCgg                                                                                     <  2:4142535/93‑1 (MQ=255)

     cTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGc                                                                                    >  2:484391/1‑93 (MQ=255)

     cTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGc                                                                                    <  1:2207294/93‑1 (MQ=255)

      tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCa                                                                                   <  1:5523438/93‑1 (MQ=255)

       gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAg                                                                                  >  2:7702514/1‑93 (MQ=255)

         tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGtt                                                                                >  1:3374969/1‑93 (MQ=255)

         tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGtt                                                                                >  2:4254835/1‑93 (MQ=255)

         tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGtt                                                                                >  1:134157/1‑93 (MQ=255)

           atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCg                                                                              >  1:2484855/1‑93 (MQ=255)

           atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCg                                                                              >  2:6352627/1‑93 (MQ=255)

            tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGTTCGGCAGTTCGa                                                                             >  2:3397685/1‑93 (MQ=255)

            tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGa                                                                             >  2:4566401/1‑93 (MQ=255)

            tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGa                                                                             >  1:6141653/1‑93 (MQ=255)

             gCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGaa                                                                            <  2:894967/93‑1 (MQ=255)

              cGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAAt                                                                           <  2:7227576/93‑1 (MQ=255)

               gggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATc                                                                          <  1:3027738/93‑1 (MQ=255)

               ggcttCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATc                                                                          <  2:2878676/90‑1 (MQ=255)

                ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCt                                                                         >  1:4166719/1‑93 (MQ=255)

                ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCt                                                                         >  1:7914199/1‑93 (MQ=255)

                 gTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTg                                                                        <  1:6968150/93‑1 (MQ=255)

                       ttCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATc                                                                  <  1:2281353/93‑1 (MQ=255)

                         cccGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAg                                                                <  1:3458827/93‑1 (MQ=255)

                          ccGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGc                                                               >  1:6822453/1‑93 (MQ=255)

                             cTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcacc                                                            >  1:3898384/1‑93 (MQ=255)

                              tGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcacca                                                           >  2:1108832/1‑93 (MQ=255)

                               gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcaccac                                                          <  2:5413173/93‑1 (MQ=255)

                               gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcaccac                                                          >  2:3579513/1‑93 (MQ=255)

                               gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcaccac                                                          <  2:7850421/93‑1 (MQ=255)

                                 ccGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActt                                                        <  2:6068199/93‑1 (MQ=255)

                                 ccGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActt                                                        <  2:1676211/93‑1 (MQ=255)

                                 ccGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGCGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActt                                                        >  2:2565249/1‑93 (MQ=255)

                                   gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtct                                                      <  1:995758/93‑1 (MQ=255)

                                   gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtct                                                      <  1:2719540/93‑1 (MQ=255)

                                    cactaaGATTTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGTTGGGGTGGGCAGTTCGAATCTGCCTATCAGCACCActtctt                                                     <  1:7414007/89‑1 (MQ=18)

                                    cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtctt                                                     <  1:633673/93‑1 (MQ=255)

                                     tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCttt                                                    >  2:6892976/1‑93 (MQ=255)

                                      ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   >  2:7898074/1‑93 (MQ=255)

                                      ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   >  1:3990176/1‑93 (MQ=255)

                                      ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   >  1:7624404/1‑93 (MQ=255)

                                       cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtc                                                  >  2:6173854/1‑93 (MQ=255)

                                       cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtc                                                  <  2:7755362/93‑1 (MQ=255)

                                        aaGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtct                                                 >  2:214937/1‑93 (MQ=255)

                                          gATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTctcc                                               <  2:1359262/93‑1 (MQ=255)

                                          gATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTctcc                                               <  1:4584971/93‑1 (MQ=255)

                                            tgtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTctcctc                                             >  1:4776623/1‑93 (MQ=255)

                                            agtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTctcctc                                             <  1:2384689/92‑1 (MQ=255)

                                              tgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTccc                                           >  2:3719206/1‑93 (MQ=255)

                                               gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCt                                          <  1:7779425/93‑1 (MQ=255)

                                               gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCt                                          <  1:5096528/93‑1 (MQ=255)

                                               gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTCTCTCCTCCCt                                          <  1:2805140/93‑1 (MQ=255)

                                                 tGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                        >  2:1414610/1‑93 (MQ=255)

                                                 tGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                        >  1:181243/1‑93 (MQ=255)

                                                 tGATATAGCTCAGATGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                        >  1:5228320/1‑93 (MQ=255)

                                                  gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtt                                       >  2:7124890/1‑93 (MQ=255)

                                                   atatAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttt                                      >  2:3819880/1‑93 (MQ=255)

                                                    tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttt                                     >  2:1981489/1‑93 (MQ=255)

                                                      taGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttt                                   >  2:3216498/1‑93 (MQ=255)

                                                      taGCTCAGTTGGTAGAGCGCACCCTTGGGAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttt                                   >  1:4246817/1‑93 (MQ=255)

                                                       aGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCCCCCTGTTTTTTc                                  <  1:2982091/93‑1 (MQ=255)

                                                        gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTcc                                 <  1:1471949/93‑1 (MQ=255)

                                                        gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTcc                                 <  1:6569380/93‑1 (MQ=255)

                                                        gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTcc                                 >  1:4847489/1‑93 (MQ=255)

                                                        gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTcc                                 <  2:7381997/93‑1 (MQ=255)

                                                          tCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCtt                               >  1:5686302/1‑93 (MQ=255)

                                                           cAGTTGGTAGAGCGCGCCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  1:428939/93‑1 (MQ=255)

                                                           cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  2:3269620/93‑1 (MQ=255)

                                                           cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  1:6430233/93‑1 (MQ=255)

                                                           cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  2:6646952/93‑1 (MQ=255)

                                                           cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  2:6522307/93‑1 (MQ=255)

                                                           cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  2:2971611/93‑1 (MQ=255)

                                                            aGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCt                             >  1:1961098/1‑93 (MQ=255)

                                                             gTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                            <  1:4450640/93‑1 (MQ=255)

                                                              ttGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGt                           <  2:678161/93‑1 (MQ=255)

                                                               tGGTAGAGCGCCCCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTACTCCTCCCAGTTTTTTCCTTCTGtt                          <  2:7046177/93‑1 (MQ=25)

                                                                ggTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCACCCCCCTGTTTTTTCCTTCTGttt                         >  1:1642017/1‑93 (MQ=255)

                                                                 gTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTa                        <  2:1839913/93‑1 (MQ=255)

                                                                    gagCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                     <  1:1033319/93‑1 (MQ=255)

                                                                     agCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGc                    <  1:7173475/93‑1 (MQ=255)

                                                                      gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   <  2:5968836/93‑1 (MQ=255)

                                                                        gcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                 >  2:5764857/1‑93 (MQ=255)

                                                                        gcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                 >  2:5320197/1‑93 (MQ=255)

                                                                        gcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                 <  1:4919509/93‑1 (MQ=255)

                                                                           cccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaa              <  1:6021129/93‑1 (MQ=255)

                                                                           cccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaa              <  2:2484854/93‑1 (MQ=255)

                                                                           cccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaa              >  2:2400080/1‑93 (MQ=255)

                                                                              ttGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacag           >  2:6270659/1‑92 (MQ=255)

                                                                              ttGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacag           >  2:394693/1‑92 (MQ=255)

                                                                              ttGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa           <  2:450288/93‑1 (MQ=255)

                                                                              ttGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa           <  1:7708357/93‑1 (MQ=255)

                                                                                 gTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc        <  2:5909186/93‑1 (MQ=255)

                                                                                  tAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCg       >  2:1268007/1‑93 (MQ=255)

                                                                                  tAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCg       <  2:4058006/93‑1 (MQ=255)

                                                                                   aaGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg      >  2:5163386/1‑93 (MQ=255)

                                                                                    aGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCggg     <  2:3166826/93‑1 (MQ=255)

                                                                                     gggTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAGGTCGGGc    >  2:4837710/1‑93 (MQ=255)

                                                                                      ggTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa   <  1:878976/93‑1 (MQ=255)

                                                                                      ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa   <  1:1411049/93‑1 (MQ=255)

CCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAT  >  NC_000913/4173286‑4173465

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