BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,186,123	C→T	R353C (CGC→TGC)	pepT →	putative peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,186,123	0	C	T	99.6%	1020.7 / ‑6.7	278	R353C (CGC→TGC)	pepT	putative peptidase T
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (136/141); total (137/141)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

CGCGAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAT  >  NC_000913/1186033‑1186212

cgcgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCtgc                                                                                         <  1:7886183/93‑1 (MQ=255)

 gcgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCtgcg                                                                                        >  1:3151346/1‑93 (MQ=255)

  cgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCgg                                                                                       <  2:2606840/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  1:6825206/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  2:1794322/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCggt                                                                                      >  2:9210898/1‑93 (MQ=255)

     gaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCGCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtgg                                                                                    <  1:895560/93‑1 (MQ=255)

      aaaGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtggt                                                                                   <  1:5149975/93‑1 (MQ=255)

      aaaGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtggt                                                                                   >  1:4853346/1‑93 (MQ=255)

         gTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTAcc                                                                                >  2:9617656/1‑93 (MQ=255)

         gTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTAcc                                                                                >  1:3905446/1‑93 (MQ=255)

           ggTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGa                                                                              <  2:6577163/93‑1 (MQ=255)

           ggTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGa                                                                              <  2:5454843/93‑1 (MQ=255)

           ggTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGa                                                                              <  1:9455153/93‑1 (MQ=255)

             ttGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACg                                                                            >  2:9813734/1‑93 (MQ=255)

              tGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACgg                                                                           >  1:9940626/1‑93 (MQ=255)

              tGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACgg                                                                           >  1:2657226/1‑93 (MQ=255)

              tGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACgg                                                                           <  2:9565014/93‑1 (MQ=255)

                    tCCGCATATTCTCGATATCGCCCAGCAGGCGATGGGCGATTGCGATATTGAACCGGAACTGACACCGATCTGCGGTGGTACCGACGGCGCGCa                                                                     >  1:5300109/1‑93 (MQ=255)

                    tCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCa                                                                     >  2:9283496/1‑93 (MQ=255)

                    tCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCa                                                                     >  1:1292784/1‑93 (MQ=255)

                      cGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGt                                                                   <  2:4714722/93‑1 (MQ=255)

                       gCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGtt                                                                  >  2:2976471/1‑93 (MQ=255)

                       gCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGtt                                                                  >  1:7437859/1‑93 (MQ=255)

                       gCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGtt                                                                  <  2:1035787/93‑1 (MQ=255)

                         atatTCTTGATATCGCCCAGCGGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCTCGCAGTTGt                                                                <  2:7905160/93‑1 (MQ=255)

                          tatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTc                                                               >  1:4713119/1‑93 (MQ=255)

                           atTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCg                                                              >  1:1840199/1‑93 (MQ=255)

                           atTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCg                                                              >  2:2214076/1‑93 (MQ=255)

                           atTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCg                                                              >  1:2559349/1‑93 (MQ=255)

                            ttCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGt                                                             >  2:9079666/1‑93 (MQ=255)

                             tctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGtt                                                            <  1:3489722/93‑1 (MQ=255)

                             tctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGtt                                                            <  2:9437495/93‑1 (MQ=255)

                             tctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGtt                                                            <  1:5867305/93‑1 (MQ=255)

                             tctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGtt                                                            <  1:4480638/93‑1 (MQ=255)

                               ccGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTa                                                          <  1:3423428/92‑1 (MQ=255)

                                 gATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATg                                                        <  2:8234662/93‑1 (MQ=255)

                                  atatCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATgg                                                       >  1:4802303/1‑93 (MQ=255)

                                   tatCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATggg                                                      >  2:6235962/1‑93 (MQ=255)

                                    atCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGa                                                     <  1:7302250/93‑1 (MQ=255)

                                      cGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGAtt                                                   <  2:1182815/93‑1 (MQ=255)

                                       gCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTa                                                  <  2:1446329/93‑1 (MQ=255)

                                       gCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGGGGTACCGACGGCGCGCAGGTGTCGTTTATGGGATTa                                                  >  2:5717418/1‑93 (MQ=255)

                                        cccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTAc                                                 >  2:7508700/1‑93 (MQ=255)

                                        cccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTAc                                                 >  1:9768719/1‑93 (MQ=255)

                                        cccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGGACCGACGGCGCGCAGTTGTCGTTTATGGGATTAc                                                 >  1:5478407/1‑93 (MQ=255)

                                          cagcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCg                                               <  1:3017826/93‑1 (MQ=255)

                                          cagcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCg                                               <  1:5705846/93‑1 (MQ=255)

                                          cagcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCg                                               <  2:8581919/93‑1 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             >  2:8167969/1‑93 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             >  1:3584111/1‑93 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             <  2:5283795/93‑1 (MQ=255)

                                             cagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGc                                            >  2:8800093/1‑93 (MQ=255)

                                               ggCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGccc                                          >  1:1130684/1‑93 (MQ=255)

                                               ggCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGccc                                          >  1:2453181/1‑93 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAccc                                     >  2:2775711/1‑93 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     <  2:9737271/93‑1 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     >  2:4211208/1‑93 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     <  2:3078846/93‑1 (MQ=255)

                                                      cgcgATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTg                                   <  1:8603254/93‑1 (MQ=255)

                                                       gcgATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTTACGTGCCACAACCTGt                                  >  2:3650979/1‑93 (MQ=255)

                                                       gcgATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGt                                  >  2:8376399/1‑93 (MQ=255)

                                                         gATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTc                                <  2:6670851/93‑1 (MQ=255)

                                                          aTTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCa                               >  2:2651846/1‑93 (MQ=255)

                                                            tGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACt                             >  2:2465418/1‑93 (MQ=255)

                                                             gCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTg                            <  1:6084435/93‑1 (MQ=255)

                                                               gATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGc                          <  2:1901855/93‑1 (MQ=255)

                                                               gATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGc                          >  2:4384803/1‑93 (MQ=255)

                                                               gATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGc                          >  1:7598024/1‑93 (MQ=255)

                                                                atatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCg                         >  1:3712375/1‑93 (MQ=255)

                                                                 tatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCgg                        >  2:4137915/1‑93 (MQ=255)

                                                                 tatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCgg                        >  1:4192260/1‑93 (MQ=255)

                                                                 tatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCgg                        >  1:6395323/1‑93 (MQ=255)

                                                                 tatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCgg                        <  1:7129053/93‑1 (MQ=255)

                                                                  atTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGt                       <  2:3167470/93‑1 (MQ=255)

                                                                   ttGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTTTGGGATTACCGTGCCCGGACCTGTGCACTGGAGGtt                      >  2:4831650/1‑93 (MQ=255)

                                                                   ttGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGtt                      <  2:9113577/93‑1 (MQ=255)

                                                                     gAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTAc                    <  2:1567840/93‑1 (MQ=255)

                                                                      aaCCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACa                   <  2:3385465/93‑1 (MQ=255)

                                                                       aCCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACaa                  >  1:7973245/1‑93 (MQ=255)

                                                                       aCCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACaa                  >  1:9195692/1‑93 (MQ=255)

                                                                        ccGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAAc                 <  1:8330676/93‑1 (MQ=255)

                                                                        ccGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAAc                 >  1:9075829/1‑93 (MQ=255)

                                                                        ccGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAAc                 >  2:8019658/1‑93 (MQ=255)

                                                                          ggAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTa               >  1:6947855/1‑93 (MQ=255)

                                                                          ggAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTa               <  1:1048163/93‑1 (MQ=255)

                                                                          ggAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTa               <  2:2936658/93‑1 (MQ=255)

                                                                             aCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCa            >  1:1371837/1‑93 (MQ=255)

                                                                             aCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCa            >  1:1008300/1‑93 (MQ=255)

                                                                               tGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATg          <  2:786373/93‑1 (MQ=255)

                                                                                gAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATgg         >  2:9768551/1‑93 (MQ=255)

                                                                                 aaaCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGt        <  1:7676858/93‑1 (MQ=255)

                                                                                 aaaCCGATCTGCGGTGGGACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGt        <  2:7055218/93‑1 (MQ=255)

                                                                                  aaCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCTCTGGCGGTTACAACTATCATGGTa       <  1:5771023/93‑1 (MQ=255)

                                                                                  aaCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTa       <  2:5861726/93‑1 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      <  1:9295775/93‑1 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      >  1:384844/1‑93 (MQ=255)

                                                                                     cGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGc    >  2:4272176/1‑93 (MQ=255)

                                                                                       aTCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTCCCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAt  <  1:2147038/93‑1 (MQ=255)

                                                                                       aTCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAt  >  1:5818059/1‑93 (MQ=255)

CGCGAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAT  >  NC_000913/1186033‑1186212

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