BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,225	T→A	D126E (GAT→GAA)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,225	0	T	A	85.8%	419.4 / 67.3	155	D126E (GAT→GAA)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (10/12); new base A (55/78); total (65/90)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

Only 100 of 156 total aligned reads displayed.

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTT  >  NC_000913/3765184‑3765317

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                           <  1:6190330/93‑1 (MQ=11)

cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                                           <  2:3600433/93‑1 (MQ=255)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                          <  1:6031325/93‑1 (MQ=12)

 gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                                          >  1:2279092/1‑93 (MQ=11)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                         >  1:7408718/1‑93 (MQ=11)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                         >  1:6053360/1‑93 (MQ=255)

  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                                         <  1:1793720/93‑1 (MQ=255)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                        <  1:5373614/93‑1 (MQ=255)

   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc                                        >  1:1841865/1‑93 (MQ=11)

    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa                                       <  2:2473433/93‑1 (MQ=255)

     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt                                      >  2:6552153/1‑93 (MQ=11)

      tatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa                                     >  2:6471833/1‑93 (MQ=14)

      tatGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa                                     >  1:7123895/1‑93 (MQ=255)

      tatGACCCGCTGCTGGGCCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa                                     >  2:4019234/1‑93 (MQ=14)

       atGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa                                    >  1:2565487/1‑93 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   >  2:7017575/1‑93 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   >  2:3999048/1‑93 (MQ=14)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   <  2:7332181/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   <  2:214392/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   <  2:269937/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   >  1:1388918/1‑93 (MQ=14)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   <  2:5592055/93‑1 (MQ=255)

        tGACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa                                   >  2:6810002/1‑93 (MQ=255)

         gACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGGTGGAATTTTTATCAGTATCCGTTGAATCTGGTCATCAAt                                  >  2:6672703/1‑93 (MQ=11)

         gACCCGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt                                  <  2:5679361/93‑1 (MQ=255)

           cccGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                                <  2:4420312/93‑1 (MQ=255)

           cccGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                                <  2:5993345/93‑1 (MQ=255)

           cccGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt                                <  1:503658/93‑1 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                               >  1:935215/1‑93 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                               <  1:5225013/93‑1 (MQ=255)

            ccGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa                               <  1:7343578/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                              >  1:3105227/1‑93 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                              >  2:4101477/1‑93 (MQ=17)

             cGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                              <  1:2975663/93‑1 (MQ=255)

             cGCTGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg                              >  1:1283297/1‑93 (MQ=17)

              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                             >  2:6823404/1‑93 (MQ=255)

              gctgctGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa                             <  1:6981792/93‑1 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                            >  1:3505997/1‑93 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                            <  1:6125868/93‑1 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                            <  1:117759/93‑1 (MQ=255)

               ctgctgGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt                            <  2:447573/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           <  2:7413496/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           <  2:3031367/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           <  2:268446/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           <  1:5170924/93‑1 (MQ=255)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           >  2:2438878/1‑93 (MQ=17)

                tGCTGGGGCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc                           <  2:611210/93‑1 (MQ=255)

                  ctggggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCg                         <  2:2858003/91‑1 (MQ=255)

                   tggggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGc                        >  1:1897124/2‑93 (MQ=255)

                    ggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                       >  1:4145236/1‑93 (MQ=255)

                    ggggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                       <  1:2839939/93‑1 (MQ=255)

                    ggggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCa                       >  2:3760681/1‑93 (MQ=255)

                     gggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                      >  2:4485671/1‑93 (MQ=255)

                     gggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                      <  1:6102669/93‑1 (MQ=255)

                     gggCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCaa                      >  2:4979699/1‑93 (MQ=255)

                      ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAg                     <  2:6346424/93‑1 (MQ=255)

                       gCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAgg                    <  2:6295542/93‑1 (MQ=255)

                       gCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAgg                    <  2:2296140/93‑1 (MQ=255)

                       gCGGTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAgg                    >  1:3171953/1‑93 (MQ=255)

                         ggTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                  <  1:4323997/93‑1 (MQ=255)

                         ggTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGtt                  >  1:3726873/1‑93 (MQ=255)

                          gTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                 >  1:608481/1‑93 (MQ=255)

                          gTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                 <  1:2072656/93‑1 (MQ=255)

                          gTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                 >  1:1825375/1‑93 (MQ=255)

                          gTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                 <  1:134816/93‑1 (MQ=255)

                          gTATATCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGttt                 <  2:3769481/93‑1 (MQ=255)

                            atatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg               <  2:2081231/93‑1 (MQ=255)

                            atatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg               >  1:7024445/1‑93 (MQ=255)

                            atatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg               <  2:4680053/93‑1 (MQ=255)

                            atatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTgg               >  2:5219799/1‑93 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              <  2:3867517/93‑1 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              >  2:4756223/1‑93 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              >  1:2554249/1‑93 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              >  1:2476395/1‑93 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              >  1:1461373/1‑93 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              <  2:5937212/93‑1 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              <  1:679143/93‑1 (MQ=255)

                             tatCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGt              >  1:4649515/1‑93 (MQ=255)

                              atCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGtt             <  2:2988380/93‑1 (MQ=255)

                               tCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTg            >  1:3110989/1‑93 (MQ=255)

                               tCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTg            >  2:2179458/1‑93 (MQ=255)

                               tCACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTg            <  1:3038417/93‑1 (MQ=255)

                                cACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGa           <  1:3718897/93‑1 (MQ=255)

                                cACTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGa           <  1:6711083/93‑1 (MQ=255)

                                 aCTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGat          >  2:5401869/1‑93 (MQ=255)

                                  cTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata         <  1:3174974/93‑1 (MQ=255)

                                  cTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata         >  1:6351321/1‑93 (MQ=255)

                                  cTCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGata         <  2:4976357/93‑1 (MQ=255)

                                   tCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat        >  2:6180309/1‑93 (MQ=255)

                                   tCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat        <  2:1076605/93‑1 (MQ=255)

                                   tCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat        <  2:3622777/93‑1 (MQ=255)

                                   tCAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGatat        >  1:100016/1‑93 (MQ=255)

                                    cAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGAtata       <  1:161083/93‑1 (MQ=255)

                                    cAGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGAtata       <  2:3392529/93‑1 (MQ=255)

                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaa      <  2:3510273/93‑1 (MQ=255)

                                     aGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaa      >  1:565540/1‑93 (MQ=255)

                                     aGGAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaa      <  1:1592049/93‑1 (MQ=255)

                                      ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa     <  2:2523145/93‑1 (MQ=255)

                                      ggAACCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATaaa     <  2:5206352/93‑1 (MQ=255)

                                         tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAAttt  >  1:5698086/1‑93 (MQ=255)

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTT  >  NC_000913/3765184‑3765317

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