BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,580	T→A	I389N (ATC→AAC)	katG →	catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,580	0	T	A	87.5%	861.1 / 110.3	288	I389N (ATC→AAC)	katG	catalase
Reads supporting (aligned to +/- strand): ref base T (25/11); new base A (139/113); total (164/124)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-02

Only 100 of 303 total aligned reads displayed.

ATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAA  >  NC_000913/4132489‑4132671

 tCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGTGAAGAact                                                                                           <  2:1075489/93‑1 (MQ=255)

  cccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactc                                                                                          <  1:57689/93‑1 (MQ=255)

  cccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactc                                                                                          <  1:6549137/93‑1 (MQ=255)

  cccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactc                                                                                          <  2:4144001/93‑1 (MQ=255)

   ccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactct                                                                                         <  1:935829/93‑1 (MQ=255)

    cGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAtctctc                                                                                        <  1:4210893/93‑1 (MQ=255)

    cGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAtctctc                                                                                        >  2:3920865/1‑93 (MQ=255)

     ggATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcg                                                                                       >  2:4920289/1‑93 (MQ=255)

      gATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgt                                                                                      >  2:7106739/1‑93 (MQ=255)

      gATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgt                                                                                      >  2:2161907/1‑93 (MQ=255)

       aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                     <  1:632227/93‑1 (MQ=255)

       aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                     >  2:5149667/1‑93 (MQ=255)

       aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                     >  1:3253392/1‑93 (MQ=255)

          cGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGtt                                                                                  >  2:6743096/1‑93 (MQ=255)

           gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGttt                                                                                 >  1:5872908/1‑93 (MQ=255)

           gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGttt                                                                                 >  1:166584/1‑93 (MQ=255)

            tttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTc                                                                                >  2:4465231/1‑93 (MQ=255)

            tttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTc                                                                                >  1:71751/1‑93 (MQ=255)

            tttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTc                                                                                <  1:2065487/93‑1 (MQ=255)

             ttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTcc                                                                               >  2:5955597/1‑93 (MQ=255)

             ttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTcc                                                                               <  2:1183886/93‑1 (MQ=255)

                aTCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCa                                                                            >  1:7479601/1‑93 (MQ=255)

                 tCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCaa                                                                           <  1:6416006/93‑1 (MQ=255)

                 tCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCaa                                                                           >  1:3430914/1‑93 (MQ=255)

                  ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAAc                                                                          <  1:4660474/93‑1 (MQ=255)

                  ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAAc                                                                          >  1:1656691/1‑93 (MQ=255)

                  ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAAc                                                                          <  2:3862921/93‑1 (MQ=255)

                    gTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGa                                                                        >  2:1455790/1‑93 (MQ=255)

                     tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGAt                                                                       <  1:7418244/93‑1 (MQ=255)

                     tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGAt                                                                       >  2:7489531/1‑93 (MQ=255)

                      cGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATc                                                                      <  2:7415526/93‑1 (MQ=255)

                       nnagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                     <  2:1574351/91‑1 (MQ=255)

                       gaagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                     >  1:137284/1‑93 (MQ=255)

                         agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGc                                                                   >  1:1276153/1‑93 (MQ=255)

                          gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                  >  2:2416016/1‑93 (MQ=255)

                          gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                  <  2:2103729/93‑1 (MQ=255)

                          gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                  >  2:497736/1‑93 (MQ=255)

                          gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                  <  2:6804725/93‑1 (MQ=255)

                          gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                  <  1:5222723/93‑1 (MQ=255)

                           aaaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAg                                                                 <  2:247414/93‑1 (MQ=255)

                             aCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGc                                                               <  2:5094690/93‑1 (MQ=255)

                             aCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGc                                                               >  1:7038476/1‑93 (MQ=255)

                               gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGt                                                             >  1:6720417/1‑93 (MQ=255)

                                tAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGtt                                                            >  2:5214731/1‑93 (MQ=255)

                                     cGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACg                                                       >  2:3821716/1‑93 (MQ=255)

                                      gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGa                                                      >  1:4796607/1‑93 (MQ=255)

                                       aCAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGaa                                                     >  2:1512192/1‑93 (MQ=255)

                                        cAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAg                                                    >  1:2573715/1‑93 (MQ=255)

                                         aaTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGc                                                   >  1:4578969/1‑93 (MQ=255)

                                          aTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGcc                                                  <  2:329740/93‑1 (MQ=255)

                                           tGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCt                                                 >  1:434112/1‑93 (MQ=255)

                                           tGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCt                                                 <  2:6580632/93‑1 (MQ=255)

                                            gCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGTGTTCAACGAAGCCtt                                                >  2:2549569/1‑93 (MQ=255)

                                            gCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCtt                                                <  1:2116697/93‑1 (MQ=255)

                                              tGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTg                                              <  2:5649400/93‑1 (MQ=255)

                                              tGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTg                                              >  2:5377829/1‑93 (MQ=255)

                                               ggTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGc                                             <  1:6014726/93‑1 (MQ=255)

                                                  gaccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCg                                          <  1:5339920/93‑1 (MQ=255)

                                                    ccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTg                                        >  2:3708225/1‑93 (MQ=255)

                                                       accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                     >  1:3607442/1‑93 (MQ=255)

                                                       accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCCCAGGCGTCCACGGAGGCCTTTGCCCGTGCCt                                     >  2:5998464/1‑93 (MQ=255)

                                                        ccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTg                                    <  1:2277345/93‑1 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  >  2:2880357/1‑93 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  <  2:5478150/93‑1 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  >  2:4288963/1‑93 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  >  2:828500/1‑93 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  <  2:1620292/93‑1 (MQ=255)

                                                          tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                  >  1:4135824/1‑93 (MQ=255)

                                                           gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                                 >  2:6485612/1‑93 (MQ=255)

                                                              gCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCaa                              >  1:2980582/1‑93 (MQ=255)

                                                               cTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCaaa                             <  1:1605968/93‑1 (MQ=255)

                                                                   gTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                         <  1:6500203/93‑1 (MQ=255)

                                                                   gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                         <  2:192509/93‑1 (MQ=255)

                                                                   gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                         <  1:5630589/93‑1 (MQ=255)

                                                                    ttttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAc                        <  2:3145511/93‑1 (MQ=255)

                                                                    ttttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCCGGTTCAAACTGAc                        >  2:5352324/1‑93 (MQ=255)

                                                                     tttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACg                       >  2:5413919/1‑93 (MQ=255)

                                                                      tttaTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGc                      <  2:2337347/90‑1 (MQ=255)

                                                                       ttaTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCCCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGca                     <  1:379996/91‑1 (MQ=255)

                                                                       tGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGca                     >  2:1421134/1‑93 (MQ=255)

                                                                        gATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGTGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcac                    >  1:6279309/1‑93 (MQ=255)

                                                                         aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGTGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                   <  1:5028390/93‑1 (MQ=255)

                                                                         aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                   >  1:5432391/1‑93 (MQ=255)

                                                                          tCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGGTCAAACCTACGccaag                  >  1:2374593/1‑89 (MQ=255)

                                                                          tCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGTGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAg                  >  2:1756869/1‑93 (MQ=255)

                                                                          tCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAg                  <  2:4700371/93‑1 (MQ=255)

                                                                           ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg                 >  2:2780750/1‑93 (MQ=255)

                                                                           ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg                 >  1:3325356/1‑93 (MQ=255)

                                                                           ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg                 <  1:5141586/93‑1 (MQ=255)

                                                                            cTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAggg                <  1:98088/93‑1 (MQ=255)

                                                                            cTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAggg                >  1:2532334/1‑93 (MQ=255)

                                                                             tGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGa               >  2:6127547/1‑93 (MQ=255)

                                                                               aGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGata             <  1:1013623/93‑1 (MQ=255)

                                                                                gTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGatat            >  1:2429654/1‑93 (MQ=255)

                                                                                   cGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATggg         <  1:1809694/93‑1 (MQ=255)

                                                                                    gagaAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgggg        >  1:2884825/1‑93 (MQ=255)

                                                                                    gagaAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgggg        >  2:338890/1‑93 (MQ=255)

                                                                                     agaagaACTCTCGTCGTTTTCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGc       >  2:6867482/1‑93 (MQ=255)

                                                                                      gaagaaCTCTCGTCGTTTCCTCAACGATCCGCAGGTGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGcc      >  1:4881459/1‑93 (MQ=255)

                                                                                      gaagaaCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGcc      <  1:1224674/93‑1 (MQ=255)

ATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAA  >  NC_000913/4132489‑4132671

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