BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,461,874	T→C	P104P (CCT→CCC)	b1399 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,461,874	0	T	C	99.0%	961.9 / ‑2.9	305	P104P (CCT→CCC)	b1399	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (145/157); total (145/160)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 320 total aligned reads displayed.

CATCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGG  >  NC_000913/1461784‑1461963

   cGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATg                                                                                      <  1:2792744/93‑1 (MQ=255)

    gggCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATgg                                                                                     >  1:610150/1‑93 (MQ=255)

    gggCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATgg                                                                                     <  1:7880809/93‑1 (MQ=255)

     ggCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATggg                                                                                    >  2:3053051/1‑93 (MQ=255)

      gCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGa                                                                                   >  2:5348860/1‑93 (MQ=255)

         aCGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATgg                                                                                >  2:285538/1‑93 (MQ=255)

         aCGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATgg                                                                                >  2:7990690/1‑93 (MQ=255)

         aCGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATgg                                                                                >  2:2003434/1‑93 (MQ=255)

           ggaGTTTCTATAGCCTCAGTGATAAAGGCTTGCGTCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTa                                                                              <  1:5353922/91‑1 (MQ=255)

           gCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCCTGGGGTGGTa                                                                              >  1:6024877/1‑93 (MQ=255)

           gCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTa                                                                              >  2:8800701/1‑93 (MQ=255)

           gCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTa                                                                              >  2:102192/1‑93 (MQ=255)

           gCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGAAGAGCAACCCGCATGGGATGGTa                                                                              >  1:8067285/1‑93 (MQ=255)

             aGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTaaa                                                                            >  1:9776799/1‑93 (MQ=255)

             aGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTaaa                                                                            <  1:6776217/93‑1 (MQ=255)

             aGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTaaa                                                                            <  1:4949599/93‑1 (MQ=255)

              gTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTaaag                                                                           >  1:9584472/1‑92 (MQ=255)

               tttCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGGCACCCGCATGGGATAGTGAATg                                                                          >  2:3739126/1‑93 (MQ=255)

                ttCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATgg                                                                         <  1:557238/93‑1 (MQ=255)

                 tCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGc                                                                        <  1:4139823/93‑1 (MQ=255)

                   tataGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTc                                                                      >  1:3319574/1‑93 (MQ=255)

                   tataGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTc                                                                      >  1:4799272/1‑93 (MQ=255)

                      aGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTg                                                                   <  1:8909065/93‑1 (MQ=255)

                       gCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGt                                                                  >  1:2574243/1‑93 (MQ=255)

                          tCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGc                                                               >  1:3295449/1‑93 (MQ=255)

                          tCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGc                                                               <  2:6135327/93‑1 (MQ=255)

                           cAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGct                                                              >  1:716141/1‑93 (MQ=255)

                            aGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGctc                                                             >  2:1805523/1‑93 (MQ=255)

                             gTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGctct                                                            >  2:6533392/1‑93 (MQ=255)

                               gATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCg                                                          <  1:5245332/93‑1 (MQ=255)

                               gATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCg                                                          <  1:4072606/93‑1 (MQ=255)

                                 tAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGa                                                        >  2:163101/1‑93 (MQ=255)

                                 tAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGa                                                        <  2:8459899/93‑1 (MQ=255)

                                  aaaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGaa                                                       <  1:868939/93‑1 (MQ=255)

                                  aaaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGaa                                                       >  1:3543197/1‑93 (MQ=255)

                                   aaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAg                                                      <  2:1423627/93‑1 (MQ=255)

                                        ttGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTa                                                 <  1:777118/93‑1 (MQ=255)

                                         tGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAg                                                >  1:9071208/1‑93 (MQ=255)

                                            gcCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATa                                             <  2:9834329/93‑1 (MQ=255)

                                             ccTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATaa                                            >  1:5189947/1‑93 (MQ=255)

                                              cTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATaaa                                           <  2:4670340/93‑1 (MQ=255)

                                               tGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAt                                          >  1:4363333/1‑93 (MQ=255)

                                               tGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAt                                          >  1:2136915/1‑93 (MQ=255)

                                                gACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAta                                         >  2:8650971/1‑92 (MQ=255)

                                                gACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATc                                         >  1:3083657/1‑93 (MQ=255)

                                                 aCGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCa                                        <  1:4634358/93‑1 (MQ=255)

                                                  cgcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCaa                                       <  1:2657847/93‑1 (MQ=255)

                                                  cgcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCaa                                       <  2:6354577/93‑1 (MQ=255)

                                                   gcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAAc                                      <  2:974131/93‑1 (MQ=255)

                                                     gACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgc                                    <  2:8318541/93‑1 (MQ=255)

                                                     gACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgc                                    >  1:5884961/1‑93 (MQ=255)

                                                     caCGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgc                                    <  2:8508020/92‑1 (MQ=255)

                                                      aCGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgct                                   >  1:318269/1‑93 (MQ=255)

                                                       cGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctg                                  >  1:4574175/1‑93 (MQ=255)

                                                        gggCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctgg                                 >  1:5385728/1‑93 (MQ=255)

                                                        gggCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctgg                                 >  2:5516021/1‑93 (MQ=255)

                                                         ggCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctggc                                >  1:3737528/1‑93 (MQ=255)

                                                         ggCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctggc                                >  2:7607850/1‑93 (MQ=255)

                                                         ggCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctggc                                >  2:3202434/1‑93 (MQ=255)

                                                          gCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctggct                               <  2:518400/93‑1 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             <  1:7561870/93‑1 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             >  1:7372317/1‑93 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             <  2:4030722/93‑1 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            >  2:362972/1‑93 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            <  1:7063874/93‑1 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGATTAGATAAATCAACGCTGGCTGAt                            >  2:6660348/1‑93 (MQ=255)

                                                              aaaGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATg                           >  2:5759736/1‑93 (MQ=255)

                                                                aGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTc                         <  2:2248665/93‑1 (MQ=255)

                                                                aGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTc                         <  2:3318772/93‑1 (MQ=255)

                                                                 gTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCa                        >  1:5741557/1‑93 (MQ=255)

                                                                  tAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaa                       >  2:7341391/1‑93 (MQ=255)

                                                                   aaaaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaa                      >  1:7212519/1‑93 (MQ=255)

                                                                    aaaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaaa                     <  1:324750/93‑1 (MQ=255)

                                                                    aaaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaaa                     >  1:2038839/1‑93 (MQ=255)

                                                                    aaaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGATGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaaa                     <  1:10043837/93‑1 (MQ=255)

                                                                       tttATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGc                  >  2:2823326/1‑93 (MQ=255)

                                                                       tttATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGc                  <  2:1267860/93‑1 (MQ=255)

                                                                           tCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGtt              <  1:2396024/93‑1 (MQ=255)

                                                                           tCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGtt              <  1:3890799/93‑1 (MQ=255)

                                                                            cgcgCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTg             >  1:2319340/1‑93 (MQ=255)

                                                                            cgcgCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTg             <  2:3903322/93‑1 (MQ=255)

                                                                            cgcgCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTg             >  1:6302752/1‑93 (MQ=255)

                                                                             gcgcAGTGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGa            >  2:7207740/1‑93 (MQ=255)

                                                                             gcgcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGa            <  2:7043891/93‑1 (MQ=255)

                                                                             gcgcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGa            <  1:7973448/93‑1 (MQ=255)

                                                                              cgcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGAt           <  2:6244637/93‑1 (MQ=255)

                                                                               tcAGAGCAACCCGCCTGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATc          >  2:3525225/2‑93 (MQ=255)

                                                                               gcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTTAAAAAGCAGTTGATc          >  1:8278586/1‑93 (MQ=255)

                                                                               gcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATc          <  2:9357578/93‑1 (MQ=255)

                                                                                cAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCt         <  1:1887836/93‑1 (MQ=255)

                                                                                 agagCAACCCGCATGGGGTGGGAAAAGGCTCCTGTTGCTCTCGGAAGGGTTTGATAAATTAACGCTGGCTGATGTTAAAAAAAAGTTGATCTg        >  2:9611264/1‑93 (MQ=255)

                                                                                 agagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTg        >  2:5646339/1‑93 (MQ=255)

                                                                                 agagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTg        <  2:9100821/93‑1 (MQ=255)

                                                                                  gagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCtcg       >  1:8262188/1‑91 (MQ=255)

                                                                                  gagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTgg       >  2:4313776/1‑93 (MQ=255)

                                                                                   agCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGc      >  2:4864842/1‑93 (MQ=255)

                                                                                     cAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCaa    <  2:4484383/93‑1 (MQ=255)

                                                                                      aaCCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAg   >  1:149642/1‑93 (MQ=255)

                                                                                       aCCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAgg  >  2:106758/1‑93 (MQ=255)

                                                                                       aCCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAgg  >  1:9309470/1‑93 (MQ=255)

CATCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGG  >  NC_000913/1461784‑1461963

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