BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,133,368	G→T	E652* (GAG→TAG)	katG →	catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,133,368	0	G	T	98.4%	939.5 / 4.4	256	E652* (GAG→TAG)	katG	catalase
Reads supporting (aligned to +/- strand): ref base G (3/1); new base T (136/116); total (139/117)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 261 total aligned reads displayed.

CAACTTCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACGAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGAT  >  NC_000913/4133277‑4133457

 aaCTTCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAg                                                                                         >  1:5624270/1‑93 (MQ=255)

   cTTCGATGGCAGCAAAAACGGCGTCTTCACTGACTGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTg                                                                                       >  2:2658772/1‑93 (MQ=255)

   cTTCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTg                                                                                       >  1:574544/1‑93 (MQ=255)

   cTTCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTg                                                                                       >  1:2497388/1‑93 (MQ=255)

    ttCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTgg                                                                                      <  2:5708378/93‑1 (MQ=255)

     tCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGGTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTGGTGGa                                                                                     >  2:2072032/1‑93 (MQ=255)

        aTGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAg                                                                                  <  1:1478342/93‑1 (MQ=255)

           gcagcaAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGcga                                                                               <  2:2380522/93‑1 (MQ=255)

            cagcaAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGcgac                                                                              <  1:5692120/93‑1 (MQ=255)

            cagcaAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGcgac                                                                              <  2:296103/93‑1 (MQ=255)

              gcaAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGcgaccg                                                                            >  1:710180/1‑93 (MQ=255)

               caAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGcgaccga                                                                           <  2:2750272/93‑1 (MQ=255)

                 aaaaCGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACcgacg                                                                         >  1:4783324/1‑93 (MQ=255)

                 aaaaCGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACcgacg                                                                         >  1:636233/1‑93 (MQ=255)

                 aaaaCGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACcgacg                                                                         >  2:2970975/1‑93 (MQ=255)

                   aaCGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGaa                                                                       <  2:1435999/93‑1 (MQ=255)

                   aaCGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGaa                                                                       >  2:2867598/1‑93 (MQ=255)

                     cGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATc                                                                     <  2:1857795/93‑1 (MQ=255)

                      ggCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCg                                                                    >  2:5329706/1‑93 (MQ=255)

                         gTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGCTACTAACGGAAAGCGACAGACGCATCGaaa                                                                 >  1:1586173/1‑93 (MQ=255)

                          tCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAag                                                                >  2:3622276/1‑93 (MQ=255)

                          tCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAag                                                                >  2:1328187/1‑93 (MQ=255)

                          tCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAag                                                                >  1:3223574/1‑93 (MQ=255)

                           cTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAaga                                                               >  2:892893/1‑93 (MQ=255)

                             tCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGc                                                             <  2:3375867/93‑1 (MQ=255)

                             tCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGc                                                             <  1:4440278/93‑1 (MQ=255)

                              cACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCt                                                            >  1:5137230/1‑93 (MQ=255)

                               aCTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTg                                                           >  1:2453270/1‑93 (MQ=255)

                               aCTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTg                                                           >  1:4343856/1‑93 (MQ=255)

                                cTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGt                                                          <  2:2515469/93‑1 (MQ=255)

                                cTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGt                                                          <  1:6068130/93‑1 (MQ=255)

                                 tGACGGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGtt                                                         <  2:3131253/93‑1 (MQ=255)

                                  gACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTc                                                        <  1:5414168/93‑1 (MQ=255)

                                  gACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTc                                                        >  1:5224713/1‑93 (MQ=255)

                                   aCCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCg                                                       <  1:2920366/93‑1 (MQ=255)

                                    ccGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGa                                                      >  1:496117/1‑93 (MQ=255)

                                     cgcgTTGTCGTATTGATCAATGTCGTCTTCGTGAACTTGCTGGATATGCGTTACTAGGGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGaa                                                     <  1:3811444/93‑1 (MQ=255)

                                     cgcgTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGTCCGACGAATCGAAAGAGCTGTTCGaa                                                     <  2:4019212/93‑1 (MQ=255)

                                     cgcgTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGaa                                                     >  1:1781934/1‑93 (MQ=255)

                                      gcgTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAg                                                    <  1:3978871/93‑1 (MQ=255)

                                      gcgTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAg                                                    >  1:1212500/1‑93 (MQ=255)

                                      gcgTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATACGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAg                                                    <  2:4215847/93‑1 (MQ=255)

                                          tGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccg                                                >  2:350414/1‑93 (MQ=255)

                                          tGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccg                                                >  2:165104/1‑93 (MQ=255)

                                          tGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccg                                                <  2:598676/93‑1 (MQ=255)

                                          tGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccg                                                >  2:1073783/1‑93 (MQ=255)

                                          tGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccg                                                <  1:4001292/93‑1 (MQ=255)

                                            gCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtg                                              <  1:5811054/93‑1 (MQ=255)

                                              gTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgac                                            <  2:5877479/93‑1 (MQ=255)

                                               tATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgacc                                           >  2:730780/1‑93 (MQ=255)

                                               tATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgacc                                           >  2:2447893/1‑93 (MQ=255)

                                                 ttGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgaccgt                                         >  2:6128441/1‑93 (MQ=255)

                                                 ttGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgaccgt                                         <  1:1951338/93‑1 (MQ=255)

                                                  tGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGccgtgaccgtg                                        >  1:1002448/1‑93 (MQ=255)

                                                    aGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGaa                                      >  2:5538548/1‑93 (MQ=255)

                                                     gCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAATGAGCTGTTCGAAGGCCGTGACCGTGaaa                                     >  2:4349216/1‑93 (MQ=255)

                                                     gCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGaaa                                     <  1:3193562/93‑1 (MQ=255)

                                                       caTGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAAcc                                   <  2:4767683/92‑1 (MQ=255)

                                                       aaTGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAAcg                                   >  1:6080504/1‑92 (MQ=255)

                                                       aaTGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAAcc                                   >  2:5404573/1‑93 (MQ=255)

                                                        aTGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCg                                  >  1:1025432/1‑93 (MQ=255)

                                                        aTGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCg                                  >  2:3865927/1‑93 (MQ=255)

                                                         tGACTTCTTCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCgg                                 <  1:2505853/93‑1 (MQ=255)

                                                            cttcttCGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGa                              >  2:562518/1‑93 (MQ=255)

                                                             ttcttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGaa                             >  1:3580395/1‑93 (MQ=255)

                                                             ttcttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGaa                             >  2:1781071/1‑93 (MQ=255)

                                                             ttcttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGaa                             >  1:2842199/1‑93 (MQ=255)

                                                              tcttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAg                            >  1:467368/1‑93 (MQ=255)

                                                              tcttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAg                            >  1:3698318/1‑93 (MQ=255)

                                                               cttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGt                           <  2:5010907/93‑1 (MQ=255)

                                                               cttcGTGAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGt                           <  1:2565554/93‑1 (MQ=255)

                                                                ttcGTGAACTTGCTGGATATGCGTTACTTGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTg                          >  2:5777060/1‑93 (MQ=255)

                                                                     gAACTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAAtt                     <  1:210707/93‑1 (MQ=255)

                                                                       agttgttgGATATGCGTTAATAGTGGAAAGCGAACGACGAATCGAAAGAGCTGTTCGAAGGCAGTGACCGTGAAACCGGCGAAGTGAAATTTa                   <  2:2563032/91‑1 (MQ=255)

                                                                       aCTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAATTGGCGAAGTGAAATTTa                   >  2:5499364/1‑93 (MQ=255)

                                                                       aCTTGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTa                   <  1:4527879/93‑1 (MQ=255)

                                                                         ttGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACg                 >  2:2033618/1‑93 (MQ=255)

                                                                          tGCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACgg                >  2:5652497/1‑93 (MQ=255)

                                                                           gcgggATATGCGTTACGAGTGGAGAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACAGGCGAAGTGAAATTTACGGc               <  1:4219508/90‑1 (MQ=255)

                                                                           gCTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGc               >  2:2201266/1‑93 (MQ=255)

                                                                            cTGGATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGcc              <  1:5578389/93‑1 (MQ=255)

                                                                              ggATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAg            >  1:1383464/1‑93 (MQ=255)

                                                                               gATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAgg           <  1:3071868/93‑2 (MQ=255)

                                                                               gATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGc           <  1:1243857/93‑1 (MQ=255)

                                                                               gATATGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGc           <  2:5089094/93‑1 (MQ=255)

                                                                                atatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGcc          <  2:5595799/93‑1 (MQ=255)

                                                                                atatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGcc          <  1:426022/93‑1 (MQ=255)

                                                                                atatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGcc          <  2:3206518/93‑1 (MQ=255)

                                                                                atatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGcc          <  2:4702174/93‑1 (MQ=255)

                                                                                atatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGcc          <  1:343087/93‑1 (MQ=255)

                                                                                 tatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCg         >  2:5138400/1‑93 (MQ=255)

                                                                                 tatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCg         <  2:4503446/93‑1 (MQ=255)

                                                                                 tatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCg         >  1:2189620/1‑93 (MQ=255)

                                                                                 tatGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGCGAAATTTACGGCCAGCCg         >  1:1537980/1‑93 (MQ=255)

                                                                                  atGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGt        >  2:5205617/1‑93 (MQ=255)

                                                                                   tGCGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGGGAAATTTACGGCCAGCCGTg       >  2:5290074/1‑93 (MQ=255)

                                                                                     cGTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCg     <  2:4315441/93‑1 (MQ=255)

                                                                                      gTTACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCgg    >  1:3784331/1‑93 (MQ=255)

                                                                                       ttACTAGTTGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGa   >  1:364652/1‑93 (MQ=255)

                                                                                       ttACTAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGa   <  1:2423185/93‑1 (MQ=255)

CAACTTCGATGGCAGCAAAAACGGCGTCTTCACTGACCGCGTTGGCGTATTGAGCAATGACTTCTTCGTGAACTTGCTGGATATGCGTTACGAGTGGAAAGCGACCGACGAATCGAAAGAGCTGTTCGAAGGCCGTGACCGTGAAACCGGCGAAGTGAAATTTACGGCCAGCCGTGCGGAT  >  NC_000913/4133277‑4133457

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