BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,210,945	C→T	S86S (AGC→AGT)	rpoD →	RNA polymerase sigma(70)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,210,945	0	C	T	99.3%	1068.3 / ‑6.0	306	S86S (AGC→AGT)	rpoD	RNA polymerase sigma(70)
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (149/154); total (151/155)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

GGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGA  >  NC_000913/3210855‑3211034

ggAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAgtgt                                                                                         <  1:6680871/93‑1 (MQ=255)

ggAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAgtgt                                                                                         <  2:5714902/93‑1 (MQ=255)

 gaagaaGCTCCGTATGCTGATGATCTGAGTCTGTCTGAAAACACCGCGGACTAAGATGCTGCCGAAGCCGCCGCGCAGGTTCTTTCCAGtgtg                                                                                        <  2:7769058/93‑1 (MQ=255)

 gaagaaGCGCCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGtgtg                                                                                        >  2:17591/1‑93 (MQ=255)

 gaagaaGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGtgtg                                                                                        >  2:2039585/1‑93 (MQ=255)

   agaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGa                                                                                      >  2:5013304/1‑93 (MQ=255)

    gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGGAGGGGCTTTCCATCGTGGat                                                                                     >  2:2714346/1‑92 (MQ=255)

    gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGaa                                                                                     >  1:1038644/1‑93 (MQ=255)

    gaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGaa                                                                                     <  2:3246914/93‑1 (MQ=255)

     aagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAAt                                                                                    <  2:6023088/93‑1 (MQ=255)

      agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   <  1:5378501/93‑1 (MQ=255)

      agCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATc                                                                                   <  2:2368536/93‑1 (MQ=255)

       gCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCt                                                                                  >  2:1502036/1‑93 (MQ=255)

         aCCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGa                                                                                >  1:5774567/1‑93 (MQ=255)

          ccGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaa                                                                               >  2:3963742/1‑93 (MQ=255)

          ccGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaa                                                                               >  2:739909/1‑93 (MQ=255)

          ccGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaa                                                                               <  2:5425908/93‑1 (MQ=255)

           cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              >  2:5110374/1‑93 (MQ=255)

           cGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGaaa                                                                              >  1:4939501/1‑93 (MQ=255)

            ggATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAAt                                                                             <  1:2903109/93‑1 (MQ=255)

            ggATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAAt                                                                             <  2:6540325/93‑1 (MQ=255)

            ggATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAAt                                                                             >  1:488959/1‑93 (MQ=255)

             gATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATc                                                                            <  2:2750473/93‑1 (MQ=255)

             gATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATc                                                                            >  2:7636489/1‑93 (MQ=255)

             gATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATc                                                                            >  2:4460532/1‑93 (MQ=255)

             gATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATc                                                                            <  1:5990610/93‑1 (MQ=255)

              aTGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCg                                                                           <  1:3080253/93‑1 (MQ=255)

               tGCCGATGATCTGATGCTGGCTGAAAACACCGCGGCCGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCgg                                                                          >  2:4189106/1‑93 (MQ=255)

               tGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCgg                                                                          <  1:5207912/93‑1 (MQ=255)

               tGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCgg                                                                          >  2:3405441/1‑93 (MQ=255)

                gCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCggg                                                                         >  2:6679285/1‑93 (MQ=255)

                gCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCggg                                                                         <  1:2097841/93‑1 (MQ=255)

                   gatgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGgcgc                                                                      >  1:7375382/1‑93 (MQ=255)

                   gatgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGgcgc                                                                      >  1:1979800/1‑93 (MQ=255)

                    atgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCa                                                                     <  2:1270750/93‑1 (MQ=255)

                    atgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCa                                                                     <  1:6960895/93‑1 (MQ=255)

                    atgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCa                                                                     >  1:6343818/1‑93 (MQ=255)

                    atgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCa                                                                     <  2:6425069/93‑1 (MQ=255)

                     tgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCAc                                                                    <  1:822601/93‑1 (MQ=255)

                     tgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCAc                                                                    <  1:267092/93‑1 (MQ=255)

                     tgatCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCAc                                                                    <  2:701587/93‑1 (MQ=255)

                       atCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGa                                                                  >  2:2515624/1‑93 (MQ=255)

                         cTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACt                                                                <  2:1951344/93‑1 (MQ=255)

                         cTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACt                                                                <  1:2992321/93‑1 (MQ=255)

                          tGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTg                                                               >  1:2792190/1‑93 (MQ=255)

                           gATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGa                                                              >  2:3942709/1‑93 (MQ=255)

                             tGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGAcc                                                            <  2:2792189/93‑1 (MQ=255)

                              gctggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGAccc                                                           <  1:7547553/93‑1 (MQ=255)

                               ctggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCg                                                          <  1:5087822/93‑1 (MQ=255)

                                tggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCgg                                                         >  2:5270885/1‑93 (MQ=255)

                                tggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCgg                                                         >  1:3754326/1‑93 (MQ=255)

                                tggctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCgg                                                         >  1:1251910/1‑93 (MQ=255)

                                  gctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTa                                                       <  1:6435386/93‑1 (MQ=255)

                                  gctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGtc                                                       >  2:4890156/1‑92 (MQ=255)

                                    tgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACg                                                     >  1:4843631/1‑93 (MQ=255)

                                      aaaaCACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCa                                                   <  2:1003291/93‑1 (MQ=255)

                                          cacCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTCCGCATGTa                                               >  1:5494746/1‑93 (MQ=255)

                                           acCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTAc                                              >  2:3047231/1‑93 (MQ=255)

                                           acCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTAc                                              >  1:4007172/1‑93 (MQ=255)

                                            ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             <  2:7257158/93‑1 (MQ=255)

                                            ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             >  2:1626286/1‑93 (MQ=255)

                                            ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             <  2:5599597/93‑1 (MQ=255)

                                            ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACa                                             <  1:325733/93‑1 (MQ=255)

                                             cgcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACAt                                            <  2:7846780/93‑1 (MQ=255)

                                              gcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATg                                           <  1:1125720/93‑1 (MQ=255)

                                              gcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATg                                           <  2:4808742/93‑1 (MQ=255)

                                              gcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATg                                           <  1:7183149/93‑1 (MQ=255)

                                               cgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGc                                          <  1:7815665/93‑1 (MQ=255)

                                               cgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGc                                          <  1:3498561/93‑1 (MQ=255)

                                                 gACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGt                                        >  2:5515177/1‑93 (MQ=255)

                                                     aaGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGaaa                                    >  1:1733550/1‑93 (MQ=255)

                                                      aGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                   <  1:4275711/93‑1 (MQ=255)

                                                      aGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAAt                                   >  1:4342338/1‑93 (MQ=255)

                                                       gATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATg                                  >  1:5319845/1‑93 (MQ=255)

                                                         tgctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATggg                                >  1:2479543/1‑93 (MQ=255)

                                                          gctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGc                               <  1:4754662/93‑1 (MQ=255)

                                                             gcCGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCAcc                            <  2:7605974/93‑1 (MQ=255)

                                                               cGAAGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGt                          <  2:1110904/93‑1 (MQ=255)

                                                                  aGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGa                       <  2:905611/93‑1 (MQ=255)

                                                                  aGCCGCCGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGa                       <  1:6281246/93‑1 (MQ=255)

                                                                   gccgccGCGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGaa                      >  2:4692485/1‑93 (MQ=255)

                                                                    ccgccgCGCAGGTGCTTTCCAGTGTGGCATCTGAAATCGGGCGACCGCCTGACCCGGTCCGCATGTACATGCGGGCAATGGGCACCGTTgcag                     >  1:2274766/1‑90 (MQ=255)

                                                                     cgccgcGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACt                    <  2:7088906/93‑1 (MQ=255)

                                                                       ccgcGCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGt                  <  2:753607/93‑1 (MQ=255)

                                                                        cgcgCAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGtt                 <  1:2195449/93‑1 (MQ=255)

                                                                        cgcgCAGCTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGtt                 <  1:4495133/93‑1 (MQ=255)

                                                                          cgcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGa               >  1:6803033/1‑93 (MQ=255)

                                                                           gcAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAc              >  1:4183840/1‑93 (MQ=255)

                                                                            cAGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAcc             <  1:1624476/93‑1 (MQ=255)

                                                                             aGGTGCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGAccc            >  1:750390/1‑93 (MQ=255)

                                                                                 gCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGa        >  1:6191132/1‑93 (MQ=255)

                                                                                 gCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGa        >  1:3607514/1‑93 (MQ=255)

                                                                                 gCTTTCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGa        >  2:6323099/1‑93 (MQ=255)

                                                                                   tttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg      >  1:6550395/1‑93 (MQ=255)

                                                                                   tttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg      >  2:1876265/1‑93 (MQ=255)

                                                                                   tttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAg      >  1:1427560/1‑93 (MQ=255)

                                                                                    ttCCAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAgg     <  2:5863183/93‑1 (MQ=255)

                                                                                     tCCAGAGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGc    <  1:62327/93‑1 (MQ=255)

                                                                                      ccAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCg   <  1:3905109/93‑1 (MQ=255)

                                                                                      ccAGTGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCg   <  1:5685625/93‑1 (MQ=255)

GGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGA  >  NC_000913/3210855‑3211034

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